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Abstract 

Unmodeled  Dynamics  in  Robust  Nonlinear  Control 

by 

Murat  Arcak 

Since  it  is  common  to  employ  reduced  models  for  control  design,  robustness  to 
unmodeled  dynamics  is  a  crucial  design  criterion.  Recent  advances  in  nonlinear  con¬ 
trol  theory  have  led  to  a  number  of  recursive  design  procedures  for  which  applications 
and  extensions  are  being  reported  at  an  increasing  rate.  However,  the  robustness  of 
these  designs  in  the  presence  of  unmodeled  dynamics  has  received  very  little  attention. 

The  purpose  of  this  dissertation  is  to  develop  systematic  redesign  procedures 
that  render  nonlinear  control  laws  robust  against  unmodeled  dynamics.  We  consider 
classes  of  unmodeled  dynamics  characterized  by  their  structural  properties  such  as 
input-to- state  stability ,  passivity ,  minimum  phaseness,  relative  degree ,  and  discuss 
their  destabilizing  effects  on  closed- loop  stability.  Using  recently  developed  nonlinear 
feedback  tools  such  as  nonlinear  small- gain  theorems  and  feedback  passivation ,  we 
develop  redesign  methods  for  each  class  of  unmodeled  dynamics  considered. 

Part  One  of  the  dissertation  presents  robust  redesigns  under  the  assumption 
that  the  full  state  of  the  plant  is  available  for  measurement.  Our  redesigns  start  with 
nominal  control  laws  such  as  those  designed  by  backstepping  and  forwarding ,  and 
robustify  them  to  achieve  global  asymptotic  stability  in  the  presence  of  unmodeled 
dynamics. 

Part  Two  addresses  output-feedback  design  issues  and  presents  a  new  nonlinear 
observer  design.  Compared  to  other  areas  of  nonlinear  control  theory,  progress  in 
nonlinear  output-feedback  design  has  been  slower  due  to  the  absence  of  constructive 
observer  design  methods.  For  systems  with  monotonic  nonlinearities,  we  introduce 
a  new  global  observer  design  which  results  in  a  nonlinear  observer  error  system  rep¬ 
resented  as  the  feedback  interconnection  of  a  linear  system  and  a  time-varying  mul¬ 
tivariable  sector  nonlinearity.  Using  efficient  numerical  methods  available  for  linear 
matrix  inequalities,  observer  gain  matrices  are  computed  to  satisfy  the  circle  criterion 
and,  hence,  to  drive  the  observer  error  to  zero. 

Due  to  the  absence  of  a  separation  principle  for  nonlinear  systems,  the  avail¬ 
ability  of  an  observer  does  not  mean  that  it  can  be  used  for  feedback  control.  We 
discuss  how  the  new  observer  can  be  incorporated  in  output-feedback  design,  and  pro¬ 
pose  a  small-gain  method  for  output-feedback  control  design  with  robustness  against 
unmodeled  dynamics.  The  design  is  illustrated  on  the  jet  engine  compressor  example. 
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Chapter  1 
Introduction 


Nonlinear  control  theory  has  undergone  a  period  of  significant  progress  in  the  last 
decade.  The  emergence  of  new  analysis  tools  such  as  input- to- state  stability ,  nonlinear 
small-gain  theorems,  and  the  idea  of  rendering  a  system  passive  by  feedback  have  lead 
to  systematic  design  procedures  such  as  backstepping  and  forwarding.  Although  such 
design  methods  are  specialized  for  certain  classes  of  systems,  their  extensions  and 
combined  use  enlarge  their  range  of  applicability.  Nonlinear  control  methods  have 
already  been  successfully  applied  to  the  control  of  electrical  motors,  diesel  engines, 
ships,  jet  engine  compressors,  and  are  promising  for  emerging  technologies  such  as 
microelectromechanical  systems. 

The  increasing  demand  for  nonlinear  control  makes  it  necessary  to  improve  the 
practicality  of  design  methods  by  studying  the  effects  of  system  uncertainty,  distur¬ 
bances,  incomplete  and  noisy  state  measurements.  While  efforts  in  this  direction  have 
been  successful  in  several  specific  problems  including  robustness  against  disturbances 
and  parametric  uncertainties,  progress  in  other  areas  has  been  slower. 

The  main  purpose  of  this  dissertation  is  to  study  the  robustness  of  nonlinear 
control  methods  in  the  presence  of  unmodeled  dynamics.  Because  it  is  common  to 
employ  low  order  models  for  control  design,  robustness  to  unmodeled  dynamics  is  a 
crucial  design  criterion.  In  general,  a  control  design  based  on  a  nominal  model  fails 
to  achieve  stabilization  in  the  presence  of  unmodeled  dynamics.  We  study  classes  of 
nonlinear  control  laws  such  as  those  designed  by  backstepping  and  forwarding,  and 
develop  systematic  redesign  procedures  for  their  robustification  against  unmodeled 
dynamics. 

A  common  class  of  unmodeled  dynamics  are  those  that  appear  at  the  plant 
input.  Among  this  class,  two  types  of  unmodeled  dynamics  affect  closed-loop  stabil¬ 
ity  properties  in  fundamentally  different  ways  and,  hence,  call  for  different  redesign 
strategies.  The  first  type  is  relative  degree  zero  and  minimum  phase  unmodeled  dy- 
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namics.  The  second  type  consists  of  unmodeled  dynamics  that  fail  to  meet  at  least 
one  of  the  relative  degree  zero  and  minimum  phase  conditions  and,  hence,  exhibit 
phase-lag  at  high  frequencies.  Although  such  “phase-lag”  unmodeled  dynamics  are 
common  in  actuators,  most  results  in  the  literature  only  address  relative  degree  zero 
and  minimum  phase  unmodeled  dynamics.  In  this  work  we  take  a  more  ambitious 
path  and  study  both  types. 

The  second  part  of  the  dissertation  addresses  problems  in  which  only  the  plant 
output  is  measured.  Progress  in  nonlinear  output-feedback  control  is  hindered  by  two 
obstacles.  First,  nonlinear  observers  are  available  only  for  very  restrictive  classes  of 
systems.  Next,  the  availability  of  an  observer  does  not  imply  that  it  can  be  used  for 
output-feedback  control,  because  the  separation  principle  does  not  hold.  In  this  work 
we  introduce  new  tools  for  both  observer  design  and  observer-based  control  design. 

For  systems  with  monotonic  nonlinearities,  we  introduce  a  new  global  observer 
design  which  results  in  a  nonlinear  observer  error  system  represented  as  the  feedback 
interconnection  of  a  linear  system  and  a  time-varying  multivariable  sector  nonlin¬ 
earity.  Using  linear  matrix  inequality  (LMI)  software,  observer  gain  matrices  are 
computed  to  satisfy  the  circle  criterion  and,  hence,  to  drive  the  observer  error  to  zero. 
We  discuss  how  the  new  observer  can  be  incorporated  in  output-feedback  design, 
and  propose  a  small-gain  method  for  output-feedback  control  design  with  robustness 
against  unmodeled  dynamics.  The  design  is  illustrated  on  the  jet  engine  compressor 
example. 

As  in  most  nonlinear  designs,  our  results  are  applicable  to  classes  of  systems 
characterized  by  their  structural  properties  and  types  of  nonlinearities.  In  contrast 
to  a  single  design  methodology  that  encompasses  all  nonlinear  systems  of  interest, 
our  approach  offers  the  advantage  of  exploiting  structural  properties,  and  avoiding 
conservative  results.  In  this  work  we  emphasize  global  designs,  that  is,  designs  for 
the  entire  domain  in  which  the  system  model  is  valid. 

In  Section  1.1  below,  we  review  preliminary  concepts  that  will  be  used  through¬ 
out  the  dissertation.  Section  1.2  presents  a  preview  of  the  main  topics  discussed  in 
the  dissertation.  Section  1.3  contains  a  list  of  notation  and  acronyms  used  in  the 
dissertation. 
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1.1  Nonlinear  Feedback  Concepts 


1.1.1  ISS  Small-Gain  Theorem 

For  systems  with  disturbances,  a  classical  total  stability  concept  is  due  to  Malkin 
[60],  Krasovskii  [49],  and  Hahn  [23].  Sontag  [86]  replaced  this  concept  with  a  more 
useful  concept  of  input-to- state  stability  (ISS).  The  system 

x  =  f(x,w),  /( 0,0)  =0,  (1.1) 

where  w  is  a  disturbance  input,  is  ISS  if  there  exist  a  class-/C£  function  /?(•,  •)  and  a 
class -K.  function  y(-)  such  that 

\x(t)\  <  max  j/3(|x(0)|,f)  ,7  ^sup  Jw(r)|^  |  .  (1.2) 

When  the  effect  of  the  initial  condition  ft  vanishes  as  t  — >  00,  the  remaining  term  y(-) 
is  a  class-/C  ISS- gain  of  the  system  (1.1). 

Sontag  and  Wang  [88]  showed  that  the  ISS  property  is  equivalent  to  the  existence 
of  a  positive  definite  and  proper  ISS-Lyapunov  function  V ( x )  such  that 

M  >  p(\w\)  =>  LfV(x,w)  < —cr(\x\) ,  (1.3) 

where  p(-)  and  er(-)  are  class -1C  functions.  With  this  /)(•),  the  ISS-gain  y(-)  in  (1.2)  is 
the  composition  j(-)  =  af1  o  a2  o  p(-),  where 

cri(\x\)  <  V (x)  <  a2(\x\) .  (1.4) 

Further  characterizations  of  the  ISS  property  are  presented  by  Sontag  and  Wang  [89]. 

The  concept  of  ISS  gain  led  to  a  new  version  of  the  ISS  small-gain  theorem 
by  Jiang  et  al.  [39],  which  includes  the  effect  of  initial  conditions  and  represents  an 
extension  of  an  earlier  theorem  by  Mareels  and  Hill  [61].  It  is  now  illustrated  on  the 
interconnected  subsystems 

x\  =  /i(ah,x2),  (1-5) 

X2  =  f-2(x2,x1). 

If  the  x  1  -subsystem  with  x2  as  its  input  has  ISS-gain  7i(-),  and  the  ^-subsystem  with 
x\  as  its  input  has  ISS-gain  72(-)>  then  the  interconnection  is  globally  asymptotically 
stable  (GAS)  if,  for  all  s  >  0, 

7i°72(s)<s,  (1.6) 

which  is  the  ISS  small-gain  condition.  Significant  extensions  and  design  applications 
of  the  ISS  small-gain  theorem  are  given  by  Teel  [93]. 
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1.1.2  Nonlinear  Relative  Degree  and  Zero  Dynamics 

The  development  of  nonlinear  geometric  methods  was  a  remarkable  achievement 
of  the  1980’s,  presented  in  the  books  by  Isidori  [25],  Nijmeijer  [67],  Marino  [63]  and 
in  the  numerous  papers  referenced  therein.  Here,  we  describe  two  geometric  con¬ 
cepts  that  will  be  used  throughout  the  dissertation:  nonlinear  relative  degree  and 
zero  dynamics.  These  indispensable  tools  bring  into  focus  the  common  input-output 
structure  of  linear  and  nonlinear  systems. 

For  the  nonlinear  system 


x  =  f{x)  +  g{x)u, 

y  =  h{x)  +  j(x)u ,  x  £  JRn  ,u,y  G  M ,  (1.7) 


the  relative  degree  at  a  point  x *  is  zero  if  j(x*)  /  0,  it  is  one  if  j(x*)  is  identically 
zero  on  a  neighborhood  of  x*  and  Lgh  :=  j£g(x)  /  0  at  x*.  This  is  so  because 


dh  . 

y  =  —x  =  Lfh  +  Lghu , 


(1.8) 


so  that,  if  Lgh  is  nonzero,  then  the  input  u[t)  appears  in  the  expression  for  the  first 
derivative  y(t)  of  the  output  y(t).  If  Lgh  is  zero,  we  can  differentiate  y  once  more  and 
check  whether  u  appears  in  the  expression  for  y(t),  etc. 

When  the  system  (1.7)  has  relative  degree  one,  its  input-output,  linearization  is 
performed  with  the  feedback  transformation 


u  =  (Lgh)  1(v  —  Lfh)  =$■  y  =  v, 


(1.9) 


which  cancels  the  nonlinearities  in  the  y-equation  and  converts  it  into  y  =  v.  Selecting 
new  state  coordinates  in  which  y  is  one  of  the  states,  the  remaining  n  —  1  equations 
with  y(t)  =  0  and  v(t)  =0  constitute  the  zero  dynamics,  that  is,  nonlinear  dynamics 
which  remain  when  the  output  is  kept  at  zero.  If  the  relative  degree  is  two,  then 
the  linear  part  of  the  system  is  y  =  v,  the  chain  of  two  integrators.  In  this  case  the 
zero  dynamics  are  described  by  the  remaining  n  —  2  equations  y{t)  =  y{t)  =  0  and 
v(t)  =  0. 

In  minimum  phase  systems  the  zero  dynamics  are  asymptotically  stable.  In 
weakly  minimum  phase  systems  the  zero  dynamics  are  stable,  but  not  asymptotically 
stable. 


1.1.3  Passivity  and  Feedback  Passivation 

Passivity,  as  a  feedback  concept,  was  first  used  by  Popov  [69,  70]  in  his  frequency 
domain  solution  to  the  absolute  stability  problem.  Popov’s  contribution  led  to  various 
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versions  of  the  circle  criterion  by  Sandberg  [81],  Narendra  and  Goldwyn  [66],  Zames 
[104],  and  to  the  Positive  Real  (PR)  Lemma  by  Yakubovich  [102]  and  Kalman  [40], 
For  a  minimal  realization  (A,  B,  C)  of  H (s),  the  PR  Lemma  shows  that  ReH(ju)  >  0 
is  equivalent  to  the  existence  of  a  P  =  PT  >  0  such  that 

AtP  +  PA<0  (1.10) 

BtP  =  C.  (1.11) 

Thus,  H (s)  being  passive  (positive  real)  means  that  P  satisfies  not  only  the  Lyapunov 
inequality  (1.10),  but  also  the  input-output  constraint  (1.11)  which  restricts  the  rel¬ 
ative  degree  of  H(s)  to  be  less  than  two  and  its  zeros  to  be  stable.  H(s)  is  called 
strictly  positive  real  (SPR)  if  (1.10)  holds  with  strict  inequality. 

For  a  general  nonlinear  system  H  with  state  x,  input  u  and  output  y,  Willems 
[101]  introduced  a  storage  function  S(x),  5(0)  =  0,  and  a  supply  rate  w(u,y),  and 
defined  H  as  dissipative  if  S(x(t))  <  w(u(t),y(t)).  Passivity  is  the  special  case  when 
w  =  uTy.  If  H (s)  is  passive,  then  its  storage  function  is  obtained  from  the  PR  Lemma, 
that  is,  S(x)  =  7£XT Px,  where  P  is  as  in  (l.lO)-(l.ll). 

The  use  of  passivity  as  a  nonlinear  stabilization  tool  was  made  possible  by  the 
Passivity  Theorem  of  Popov  [71]  and  Zames  [104],  which  states  that  the  feedback 
interconnection  of  two  nonlinear  passive  blocks  Hi  and  H%  in  Figure  1.1  is  stable. 
This  theorem  encompasses  the  circle  criterion  where  Hi  =  Hi(s)  is  positive  real  and 
H2  is  a  time- varying  nonlinearity  y?  =  (p(t,  U2),  which  is  passive  because  of  its  sector 
property  U2(j)(t,  u?)  >  0. 


While  passivity  was  a  common  tool  in  adaptive  control,  its  first  use  for  nonlinear 
control  design  was  motivated  by  a  difficulty  encountered  in  feedback  stabilization  of 
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linear-nonlinear  cascade  systems 

x  =  f(x,£) 

£  =  A£  +  Bu,  (1.12) 

resulting  from  input-output  linearization.  The  difficulty  was  that  GAS  of  the  subsys¬ 
tem  x  =  f(x,  0)  is  not  sufficient  to  achieve  GAS  of  the  whole  cascade  with  ^-feedback 
u  =  K£  alone.  Thus,  we  need  a  feedback  from  both  £  and  x,  that  is, 

u  =  +  v(x,g) .  (1.13) 

Such  a  control  law  was  designed  by  Byrnes  and  Isidori  [9]  for  the  special  case  of  (1.12) 
with  £  =  Bu,  where  B  is  a  square  nonsingular  matrix.  Kokotovic  and  Sussmann  [47] 
extended  this  design  to  feedback  passivation  where  the  main  idea  is  to  make  the 
cascade  (1.12)  appear  as  the  feedback  interconnection  of  the  blocks  H\  and  H-2  in 
Figure  1.1.  The  final  result  in  Figure  1.2  is  arrived  at  in  several  steps.  First,  an 
output  rj  of  the  linear  block  H \  is  selected  to  be  the  input  of  the  nonlinear  block  H-2 , 
that  is,  the  x-subsystem  of  (1.12)  is  rewritten  as 

x  =  /(x,0)  +  g(x,C)rj,  (1.14) 

where  several  choices  of  7]  =  C£  may  be  available.  An  output  y  is  then  chosen  to  render 
(1.14)  passive  from  y  to  y.  If  a  Lyapunov  function  V(x)  is  known  for  x  =  /(x,  0)  so 
that  LfV  <  0  for  all  i/O,  then  y  =  LgVT  renders  (1.14)  passive  because 

V  =  LfV  +  LgVr1<LgVr1  =  yTi1.  (1.15) 

Finally,  if  the  linear  block  H±  is  made  PR  by  feedback  /i£,  the  passivity  theorem  will 
be  satisfied  by  closing  the  loop  with  —y  as  in  Figure  1.2. 

This  means  the  nonlinear  feedback  term  in  the  control  law  (1.13)  is  v(x,£)  = 
—y  =  —LgVT.  What  remains  to  be  done  is  to  find  K  and  P  >  0  to  satisfy  the 
Feedback  PR  Lemma 


(A  +  BK)tP  +  P(A  +  BK )  <  0 , 

BtP  =  C.  (1.16) 

Kokotovic  and  Sussmann  [47]  showed  that  an  FPR  solution  exists  if  and  only  if  the 
minimal  representation  (A,B,C)  is  relative  degree  one  and  weakly  minimum  phase. 
Saberi  et  al.  [79]  showed  that  the  weak  minimum  phase  property  of  ( A,B,C )  is 
necessary  unless  some  other  restriction  is  imposed  on  the  nonlinear  part.  An  analysis 
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Figure  1.2:  Feedback  passivation  design. 

by  Sepulchre  [83] ,  extending  earlier  results  of  Sussmann  and  Kokotovic  [91] ,  revealed 
that  the  obstacle  is  the  peaking  phenomenon.  Higher  relative  degree  systems  are 
prone  to  destabilizing  transients  caused  by  fast  or  slow  peaking.  For  nonminimum 
phase  systems  global  stabilization  can  be  achieved  only  with  further  restrictions  on 
the  cross-term  g(x,  £)  as  discussed  by  Braslavsky  and  Middleton  [8],  and  Sepulchre 
and  Arcak  [84],  where  these  restrictions  are  characterized  by  a  relationship  between 
the  locations  of  the  nonminimum  phase  zeros  and  the  growth  of  g(x,£)  in  x  and  £. 

1.2  Preview  of  the  Dissertation 

1.2.1  Robustness  against  Unmodeled  Dynamics 

Control  design  methods  like  feedback  linearization  [25],  backstepping  [53],  for¬ 
warding  [93,  85]  achieve  GAS  for  classes  of  nonlinear  systems.  However,  for  the 
system 


x  =  f(x)  +  g(x)p(£,x,u)  (1.17) 

£  =  q(£,x,u), 

where  the  ^-subsystem  with  the  output  v  =  p(x,  £,  u)  represents  unmodeled  dynamics, 
a  GAS  control  law  u  =  a(x)  designed  for  the  nominal  model  x  =  f(x)  +  g(x)u,  will 
in  general  fail  to  achieve  GAS  of  the  actual  system.  Therefore,  it  is  necessary  to 
redesign  control  laws  to  make  them  robust  against  unmodeled  dynamics.  Several 
redesigns  have  been  based  on  small-gain  and  passivity  properties  of  the  unmodeled 
dynamics. 
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As  an  illustration  of  small-gain  redesigns  proposed  by  Jiang  et  al.  [73,  36,  39,  35, 
34],  Krstic  et  al.  [54,  55],  Praly  and  Wang  [76],  we  let  p(£,x,u)  =  u  +  £,  q(£,x,u)  = 
q(£,x),  and  assume  that  the  unmodeled  dynamics  are  ISS  with  x  considered  as  the 
input,  that  is, 

\m\  <  max  j/?i(|£(0)|,i),7i  f  sup  Jx(r)|^  |  .  (1.18) 

The  nominal  control  law  a(x)  was  designed  for  V(x),  such  that  Lf+gaV  :=  (f(x)  + 

g(x)a(x))  <  0,  Vx  /  0.  We  substitute  p(£,x,u)  =  u  +  £  in  (1.17) 

x  =  f(x)  +  g(x)[u  +  £] ,  (1-19) 

and  redesign  the  control  law  to  assign  an  ISS-gain  from  ^  to  x,  that  is, 

#2(|s(0)I>*)>72  (  sup  |£(r)|>)l  .  (1.20) 

Vo <T<t  J  J 

If  J2O)  is  selected  such  that  that  71  o  72(5)  <  s,  Vs  /  0,  then  the  ISS  small-gain 
theorem  of  Teel  et  al.  [39,  93]  guarantees  GAS  of  the  closed-loop  system.  The  redesign 
of  the  control  law  is  completed  by  a  continuous  approximation  of  the  control  law 

u  =  a(x)  -  sgn(LgV(x))p(\x\) ,  (1.21) 

where  p(-)  is  determined  from  the  desired  ISS-gain  72O)  and  the  Lyapunov  function 
V  (x) .  The  resulting  feedback  system  can  tolerate  all  unmodeled  dynamics  that  satisfy 
the  ISS  property  (1.18),  which  represents  its  ISS-gain  margin. 

Small-gain  redesigns  for  adaptive  control  were  proposed  by  Jiang  and  Praly 
[37,  38],  Jiang  and  Hill  [33],  and  Jiang  [32], 

An  alternative  redesign  by  passivation  does  not  require  that  unmodeled  dynam¬ 
ics  have  bounded  ISS-gain.  Instead,  the  class  of  unmodeled  dynamics  is  restricted 
by  a  passivity  requirement  on  the  ^-subsystem  in  (1.17)  with  u  as  the  input  and 
v  =  p(£,  x ,  u)  as  the  output. 

The  passivation  redesigns  of  Jankovic  et  al.  [31],  extended  by  Hamzi  and  Praly 
[24],  are  based  on  V(x)  as  a  control  Lyapunov  function  (CLF)  for  the  nominal  system 
x  =  f(x)  +  g(x)u.  For  example,  if  V(x)  has  the  property 

LfV(x)  <\LgV(x)\2 ,  Vx/0,  (1.22) 

then  the  control  law 


u  =  —kLgV(x) ,  k  >  1 


(1.23) 


guarantees  GAS  not  only  for  the  nominal  system,  but  also  for  all  stable  unmodeled 
dynamics  which  satisfy  the  dissipativity  condition 

S(£)<vu—yu2,  (1-24) 

k 

where  S(£)  is  a  storage  function.  Thus,  (1.24)  represents  a  ‘stability  margin’  for  the 
control  law  (1-23)  because  it  defines  a  class  of  admissible  unmodeled  dynamics.  This 
stability  margin  is  guaranteed  if,  for  example,  the  control  law  in  (1.23)  is  optimal  with 
the  control  penalty  matrix  R(x)  =  I  because,  then,  the  value  function  V (x)  satisfies 
(1.22).  For  the  case  when  V{x)  fails  to  satisfy  (1.22),  Jankovic  et  al.  [31]  construct  a 
new  V (x)  which  recovers  the  same  stability  margin. 

The  first  result  in  this  dissertation,  presented  in  Chapter  2,  is  a  passivation 
redesign  of  backstepping  control  laws  for 

X  =  $(x)  +  r(x)«  (!-25) 

£  =  ?(£,«)  (1-26) 

v  =  p(£,u), 


which  robustifies  them  against  unmodeled  dynamics  described  by  the  ^-subsystem 
(1.26).  To  preserve  GAS  in  the  presence  of  unmodeled  dynamics  satisfying  (1.24),  we 
want  the  CLF  V(x)  to  be  as  hi  (1-22).  It  was  shown  by  Hamzi  and  Praly  [24]  that 
(1.22)  is  equivalent  to  the  existence  of  another  CLF  V(x)  such  that 


L*V(X)  „  , 

hm  sup  .  .  9  <  l, 

x-0  (LtV(x))  2 

for  some  l  >  0.  Then,  using  a  positive  scalar  function  $(•)  such  that 


lim 


't— XX) 


l  6^ds  =  +°°’  9{y' W)  >  (LrVM)2  ’  VX  *  °  ’ 


/o  (LrV(X))2 

a  redesigned  V(x)  which  satisfies  (1.22)  is 


(1.27) 


(1.28) 


fV{x) 

V(x)=  /  8(s)ds. 


(1.29) 


Jo 

We  derive  conditions  under  which  a  CLF  Vn(x),  constructed  after  n  steps  of  backstep¬ 
ping,  will  satisfy  the  above  condition.  As  an  illustration,  such  a  CLF  is  constructed 
for  the  system 


X  =  X3  +  x 

X  =  V 

£  =  +  (1-30) 

v  =  +  u, 


9 


with  the  unmodeled  dynamics  subsystem  (1.30). 

After  two  steps  of  backstepping  we  obtain  the  CLF 

v2  =  \x2  +  \y\  y  =  x  +  X  +  X3,  (1.31) 

and  using  6(V 2)  =  3  +  18V^2 ,  we  redesign  the  control  law  to  be 

u  =  -3  y-\{X2+y2fy.  (1.32) 

This  control  law  achieves  GAS  for  all  unmodeled  dynamics  satisfying  the  dissipation 
inequality  (1.24)  with  k  =  1,  including  (1.30). 

For  linear  unmodeled  dynamics  we  present  a  version  of  this  redesign  which 
simplifies  its  applications  to  high  order  systems. 

Another  direction  for  achieving  robustness  against  a  wider  class  of  unmodeled 
dynamics  is  the  dynamic  nonlinear  damping  redesign  presented  in  Chapter  3.  This 
redesign  replaces  passivity  and  small-gain  restrictions  by  the  less  restrictive  assump¬ 
tion  that  the  unmodeled  dynamics  subsystem  is  relative  degree  zero  and  minimum 
phase.  For  the  system  (1.17),  a  nominal  control  law  ot{x)  is  redesigned  to  be 

u  =  a(x)  —  k(1  +  \m(t)\  +  \a(x)\)LgV(x),  k  >  0,  (1.33) 

where  the  signal  \m(t)\  constitutes  an  upper  bound  for  the  state  of  the  unmodeled 
dynamics.  We  prove  that  this  redesign  guarantees  GAS  for  sufficiently  large  k. 

To  illustrate  the  dynamic  nonlinear  damping  redesign,  we  consider  the  system 

x  =  x2  +  A  (s)u,  (1-34) 

where  the  poles  A i  of  the  unmodeled  dynamics  A(s)  satisfy  Re{Xj}  <  —5  <  0.  For 
the  nominal  control  law  a(x)  =  —x  —  x2,  the  redesigned  control  law  is 

u  =  —  x  —  x2  —  k(1  +  \m\  +  \a\)x 

rh  =  —5m  +  |it|  , 

which  ensures  that  |m(f)|  is  an  upper  bound  for  the  state  £(t)  of  A(s). 

The  relative  degree  zero  and  minimum  phase  restrictions  of  the  dynamic  nor¬ 
malization  redesign  are  tight,  as  illustrated  by  the  simulation  results  in  Figure  1.3. 
While  the  closed-loop  system  is  GAS  with  the  linear  minimum  phase  unmodeled  dy¬ 
namics  Ai(s),  the  nonminimum  phase  A2(s)  causes  the  closed-loop  solutions  to  grow 
unbounded. 
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Figure  1.3:  Closed-loop  solutions  for  minimum  phase  Ai(s)  (left),  and  nonminimum 
phase  A2 (s)  (right). 

1.2.2  Phase-Lag  Unmodeled  Dynamics 

The  redesigns  discussed  so  far  require  that  the  unmodeled  dynamics  be  relative 
degree  zero  and  minimum  phase  and,  therefore,  exclude  several  practically  important 
classes  of  unmodeled  dynamics  such  as  those  due  to  actuators.  The  high-gain  nature 
of  these  redesigns  results  in  reduced  robustness  to  phase-lag  unmodeled  dynamics. 

The  problem  of  robustness  against  phase-lag  unmodeled  dynamics  has  received 
insufficient  attention  in  the  nonlinear  control  literature.  A  singular  perturbation 
result  by  Sepulchre  et  al.  [85,  Theorem  3.18]  shows  that  asymptotic  stability  can 
be  preserved  with  large  regions  of  attraction  if  the  unmodeled  dynamics  are  much 
faster  than  the  nominal  closed-loop  system.  For  unmodeled  dynamics  with  relative 
degree  greater  than  zero,  Praly  and  Jiang  [74]  have  designed  a  semiglobal  control 
law  that  incorporates  a  high-gain  observer  for  the  unmodeled  dynamics.  For  the 
single  integrator  x  =  v,  Zhang  and  Ioannou  [105]  have  considered  nonminimum  phase 
unmodeled  dynamics  and  achieved  global  asymptotic  stability  using  low-gain  control 
laws. 

In  Chapter  4,  we  achieve  global  asymptotic  stability  for  a  broader  class  of  sys¬ 
tems  with  phase-lag  unmodeled  dynamics.  We  study  systems  in  feedforward  form, 
and  redesign  nested  saturation  control  laws  of  Teel  [92,  93]. 

We  illustrate  our  result  with  a  design  for  the  system 

X'2  =  x\  +  xf  +  v2 

x\  =  v  +  v2  ,  (1.35) 

and  test  its  robustness  against  the  nonminimum  phase  unmodeled  dynamics 

v(s)  =  2  _l  +_l  i  u(s ) '  (L36) 

sz  +  s  +  1 
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For  the  nominal  system,  that  is  for  (1.35)  with  v  =  u,  the  nested  saturation  design  of 
Teel  [92,  93]  is  applicable.  With  two  saturation  functions  4>i(s)  =  sgn(s)  min{|s|,  Aj} 
where  Ai  =  0.5,  A2  =  0.05,  and  y\  =  X{,  y 2  =  x\  +  X2,  the  nested  saturation  control 
law  is 


u  =  -0i (2/1  +  <^2(2/2))  •  (1-37) 

This  control  law  is  not  robust  against  the  unmodeled  dynamics  (1.36):  the  solutions 
of  the  closed-loop  system  (1.35),  (1.36),  (1.37)  grow  unbounded  as  shown  in  Figure 
1.4.  We  redesign  the  control  law  (1.37)  to  be 

u  =  ~4>i(kixi  +  4>2{k2X2))  ■  (1.38) 

With  parameters  Ai  =  0.5,  k\  =  0.4,  A2  =  0.05,  &2  =  0.02,  designed  according  to 
a  small-gain  procedure  described  in  Chapter  4,  the  control  law  (1.38)  renders  the 
closed-loop  system  GAS,  as  illustrated  in  Figure  1.5. 


Xl 

Figure  1.5:  Robust  redesign. 


The  achieved  robustness  property  is  due  to  the  low-gain  design  (1.38),  which  is 
sufficient  for  stabilization  of  systems  in  feedforward  form  such  as  (1.35). 


1.2.3  Output-Feedback  Control 

Compared  with  advances  in  other  areas  of  nonlinear  control  theory,  progress  in 
nonlinear  output-feedback  design  has  been  slower.  First,  nonlinear  observers  are  avail¬ 
able  only  for  very  restrictive  classes  of  systems.  Second,  unlike  linear  systems  where 
the  separation  principle  allows  output-feedback  problems  to  be  solved  by  combin¬ 
ing  state- feedback  controllers  with  observers,  the  availability  of  a  nonlinear  observer 
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does  not  imply  that  it  can  be  used  for  output-feedback  design  because  the  separation 
principle  does  not  hold. 

Global  nonlinear  observer  designs  in  the  literature  severely  restrict  classes  of 
systems  and  nonlinearities.  Early  efforts  by  Thau  [98],  Kou  et  al.  [48]  and  Banks 
[3]  restricted  the  state-dependent  nonlinearities  to  be  globally  Lipschitz.  Under  this 
restriction,  quadratic  Lyapunov  functions  have  been  used  for  observer  design,  with 
various  extensions  by  Tsinias  [99],  Yaz  [103],  Boyd  et  al.  [7,  Section  7.6],  Raghavan 
and  Hedrick  [77],  Eker  and  Astrom  [14],  and  Rajamani  [78]. 

For  systems  in  which  the  nonlinearities  appear  as  functions  of  the  measured 
output,  the  observer  design  is  linear  because  the  nonlinearity  is  canceled  by  an  “output 
injection”  term.  This  class  of  systems  has  been  characterized  by  Krener  and  Isidori 
[50],  Bestle  and  Zeitz  [6],  Besangon  [5],  and  other  authors.  Output-injection  observers 
have  been  incorporated  in  observer-based  control  designs  by  Marino  and  Tornei  [62, 
63],  Kanellakopoulos  et  al.  [41,  53],  Praly  and  Jiang  [73],  and,  for  stochastic  nonlinear 
systems,  by  Deng  and  Krstic  [13,  51],  and  Arslan  and  Ba§ar  [1], 

A  broader  class  of  systems  is  characterized  by  linear  dependence  on  unmeasured 
states.  For  this  class,  dynamic  output-feedback  designs  have  been  proposed  by  Praly 
[72],  Pomet  et  al.  [68],  Marino  and  Tornei  [63],  Freeman  and  Kokotovic  [19],  and 
Praly  and  Kanellakopoulos  [75]. 

The  possibility  to  dominate  the  state  dependent  nonlinearities  by  linear  high- 
gain  has  recently  been  explored  by  Khalil  and  coworkers  [15,  46].  While  —  kx  cannot 
dominate  x 3  globally,  it  can  do  so  for  as  large  |x|  as  desired,  provided  the  gain  k 
is  sufficiently  large.  To  achieve  this  kind  of  “semiglobal”  convergence,  one  must 
avoid  the  destabilizing  effect  of  the  peaking  phenomenon ,  analyzed  by  Sussmann  and 
Kokotovic  [91] .  In  Khalil’s  high-gain  observer,  this  is  achieved  with  saturation  of  the 
observer  signals  before  they  are  fed  to  the  controller.  The  high-gain  observer  has 
been  employed  in  semiglobal  output-feedback  designs  by  Teel  and  Praly  [96],  Lin  and 
Saberi  [57],  Lin  and  Qian  [56],  Praly  and  Jiang  [74],  and  Isidori  et  al.  [28].  The  same 
approach  has  been  employed  in  adaptive  control  by  Jankovic  [29],  Khalil  [43],  and  in 
the  nonlinear  servomechanism  problem  by  Khalil  et  al.  [42,  59,  45],  and  Isidori  [26]. 

The  idea  of  using  a  high-gain  observer  with  saturation  has  led  to  the  semiglobal 
separation  theorem  of  Teel  and  Praly  [95] ,  which  states  that  global  stabilizability  by 
state-feedback  and  uniform  observability  in  the  sense  of  Gauthier  and  Bornard  [21] 
imply  semiglobal  stabilizability  by  dynamic  output  feedback.  Several  extensions  and 
interpretations  of  this  result  have  been  presented  by  Atassi  and  Khalil  [2],  Isidori  [27, 
Section  12.3],  and  Battilotti  [4], 

Global  high-gain  observers  have  been  designed  by  Gauthier  et  al.  [20]  under 
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a  global  Lipschitz  condition  -  a  common  restriction  in  most  global  designs.  In  the 
absence  of  such  a  restriction,  global  stabilization  by  output  feedback  may  not  be 
possible,  as  shown  by  the  counterexamples  of  Mazenc  et  al.  [64], 

1.2.4  A  New  Nonlinear  Observer 

In  Chapter  5  we  present  an  advance  in  observer  design  for  systems  with  state- 
dependent  nonlinearities.  In  the  new  observer,  global  convergence  is  achieved  without 
high-gain.  This  advance  is  made  under  two  restrictions  which  allow  the  observer  error 
system  to  satisfy  the  well  known  multivariable  circle  criterion.  First,  a  linear  matrix 
inequality  (LMI)  is  to  be  feasible,  which  implies  a  positive  real  property  for  the  linear 
part  of  the  observer  error  system.  The  second  restriction  is  that  the  nonlinearities  be 
nondecreasing  functions  of  linear  combinations  of  unmeasured  states.  This  restriction 
ensures  that  the  vector  time- varying  nonlinearity  in  the  observer  error  system  satisfies 
the  sector  condition  of  the  circle  criterion. 

We  illustrate  the  nonlinear  observer  design  on  the  well  known  van  der  Pol  oscil¬ 
lator 


xi  =  x2 ,  (1.39) 

1  3 

X2  =  -X'l  +  X2  ~  -X2  • 
o 

Our  problem  is  to  estimate  X2  when  only  x\  is  measured.  The  system  nonlinearity 
x3  which  violates  the  global  Lipschitz  assumption  depends  on  the  unmeasured  state. 
The  necessary  output  injection  conditions  of  Krener  and  Isidori  [50]  are  not  satisfied 
and,  hence,  none  of  the  existing  global  observer  design  methods  is  applicable. 

Our  idea  is  to  add  a  nonlinear  injection  term  w3  in  the  observer 

x\  =  x2  +  h(xi  —  x\) ,  (1-40) 

1  , 

X2  =  -X\  +  X2  +  h{x\  -  X\ )  — -w  , 

o 

so  that  the  observer  error  e  :=  x  —  x  satisfies 

ei  =  hei+e2,  (1.41) 

e2  =  hei  +  e2  —  ~{x\  -  w3) . 

The  task  of  the  nonlinear  injection  term  w3  is  to  counteract  x\  and  achieve  conver¬ 
gence  of  e(t)  to  zero.  For  w,  we  are  free  to  construct  any  function  from  the  available 
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signals  aq,  aq,  xq.  For  this  construction,  we  use  the  Lyapunov  function  V  =  eT Pe. 
Its  derivative  along  the  observer  error  system  (1.41)  is 


"  h 

1 

'hi' 

T  \  2  T 

0 

+ 

P  e - eTP 

. h 

1 

_  h  1  _ 

3 

1 

V  =  e1  |  P  |  |  +  |  |  P  |  e  -  -e1  P  \  .  |  (x%  -  w3) .  (1.42) 

To  render  its  quadratic  part  negative  definite,  we  select  P,  1 1,  I2  such  that 


P 


h  1 
h  1 


+ 


h 

h 


T 

P  +  uI  <  0, 


(1.43) 


for  some  positive  constant  u.  Then,  we  design  w  to  render  the  remainder  of  (1.42) 
nonpositive.  Denoting 


p  k 
k  m 


(1.44) 


we  see  that  w  must  guarantee  (. ke\  +  me 2)(x|  —  w 3)  >  0 .  To  this  end  we  employ  the 
inequality 


(x'2  —  w)(x\  —  w3)  >  0  VX2,  vj  G  M , 


(1.45) 


which  holds  because  x3  is  a  nondecreasing  function  of  x'2-  In  view  of  (1.45),  we 
construct  w  to  satisfy 


ke  1  +  ?ne2  =  X2  —  w, 


(1.46) 


so  that 


w  =  x 2  —  me2  —  ke  1  =  (1  —  ?n)x2  +  mi2  +  fe(xi  —  xq) .  (1-47) 

Since  X2  is  not  available,  we  select  m  =  1  and  obtain 

w  =  x 2  +  /c(xi  —  xi) ,  (1-48) 

which,  substituted  in  the  observer  (1.40),  yields  V  <  —u\e\2.  This  guarantees  that 
e(t)  — >  0  exponentially,  provided  (1.43)  can  be  satisfied  with  P  >  0,  constrained  by 
m  =  1.  We  express  this  constraint  as 


0 

V 

k  ' 

0 

’  k  ’ 

1 

k 

m 

1 

1 

(1.49) 
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and  combine  it  with  (1.43)  in  the  following  matrix  inequality: 


P 


h  1 

h  1 

(p 


+ 

0 

1 


h 

1 

1 

P  +  izI 

p 

0 

k 

h 

1 

1 

1 

’  k  " 

\T 

1 

) 

0 

<  0. 


(1.50) 


This  inequality  is  linear  in  P,  P  [l\  ,  v  and  k.  Its  feasibility  can  be  determined 

numerically  using  the  efficient  LMI  methods  (see  e.g.  Boyd  et  al.  [7]).  If  we  set  k  =  0 
in  (1.48),  that  is,  if  we  let  w  =  x 2,  then  the  LMI  is  not  feasible.  With  k  7^  0,  the  LMI 
is  feasible  and  a  solution  is 


P  = 


6 

-2 


*  =  2, 


L  = 


’  h  ' 

-6 

_  h  _ 

-16 

k  =  -2. 


(1.51) 


This  shows  that  it  is  crucial  for  the  nonlinear  injection  w  to  contain  the  output  error 
k(x  1  —  xi).  The  resulting  nonlinear  observer 


X\  =  X2  —  6(xi  —  Xl) 

X2  =  -Xl  +  x2  -  16(xi  -  Xl)  -  -(X2  -  2(xi  -  Xl))3  ,  (1.52) 

o 

guarantees  global  convergence  x(t)  — *  x(t). 

The  constraint  m  =  1  imposed  on  P  reveals  that  the  nondecreasing  property 
of  the  nonlinearity  x\  is  not  the  only  precondition  for  a  successful  design.  The  other 
required  property  is  revealed  by  rewriting  the  observer  error  system  (1.41)  as 


e  = 


h  1 
h  1 


e  + 


0 

1 


0, 


■&  :  = 


and  considering 


(1.53) 


z  :=  [k  l]e  =  X2  —  w  (1-54) 

as  the  output  of  the  linear  block.  Then,  (1.53)-(1.54)  is  represented  by  the  block- 
diagram  in  Figure  1.6.  The  key  observation  now  is  that  the  nondecreasing  property 
(1.45)  implies  that  <p(t,z)  :=  |(x|  —  w 3)  is  a  sector  nonlinearity:  z<p(t,z)  >  0.  This 
provides  a  link  with  the  well  known  circle  criterion  (see  e.g.  Khalil  [44]),  which 
guarantees  e(t)  — >  0  if  the  linear  block  is  SPR.  Indeed,  (1-43)  and  (1.49)  constitute 
the  required  SPR  condition. 
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In  Chapter  5  we  generalize  this  observer  design  to  the  control  system 

x  =  Ax  +  G'y(Hx)  +  g(y,  u)  (1.55) 

V  =  Cx, 

in  which  the  state  dependent  nonlinearity  7  {Hx)  is  a  vector  where  each  entry  is 
a  nondecreasing  function  of  a  linear  combination  of  the  states.  We  also  develop  a 
reduced-order  variant  of  this  observer,  which  provides  estimates  only  for  the  unmea¬ 
sured  states.  Both  full-order  and  reduced-order  designs  are  constructive  in  the  sense 
that  the  issues  of  existence  and  the  evaluation  of  observer  matrices  are  resolved  by 
efficient  LMI  computations. 

A  further  advantage  of  this  observer  is  its  robustness  against  inexact  modeling 
of  nonlinearities.  We  give  bounds  within  which  the  observer  error  gradually  increases 
with  an  increase  in  the  modeling  error. 

In  Chapter  6  we  present  structural  conditions  that  characterize  the  feasibility  of 
the  new  observer  design.  We  derive  necessary  and  sufficient  conditions  for  feasibility 
of  the  full-order  observer  design,  and  prove  that  they  are  equivalent  to  the  feasibility 
conditions  of  the  reduced-order  observer  design. 

1.2.5  Robust  Output-Feedback  Design 

Due  to  the  absence  of  a  separation  principle,  the  availability  of  a  nonlinear 
observer  does  not  imply  that  it  can  be  used  for  output-feedback  design.  It  was 
shown  in  Krstic  et  al.  [53,  Section  7.1]  that,  even  with  exponentially  converging 
state  estimates,  a  certainty-equivalence  control  law  may  lead  to  finite  escape  times. 
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Thus,  nonlinear  observers  can  only  be  used  in  conjunction  with  control  laws  that 
guarantee  boundedness  of  the  states  in  the  presence  of  bounded  observer  errors.  For 
classes  of  systems,  such  control  laws  have  been  designed  by  Freeman  and  Kokotovic 
[17,  18],  Krstic  et  al.  [53],  and  Marino  and  Tomei  [63]. 

In  Chapter  7,  the  new  observer  is  incorporated  in  output-feedback  design  in 
conjunction  with  the  observer  backstepping  design  of  Krstic  et  al.  [53].  In  contrast 
to  full  state-feedback  backstepping,  this  design  employs  the  observer  states  as  the 
virtual  control  variables,  and  introduces  nonlinear  damping  terms  to  counteract  the 
destabilizing  effect  of  the  observer  error. 

Finally,  we  use  our  observer-based  control  design  tools  to  develop  a  robust 
output-feedback  scheme  for  systems  with  unmodeled  dynamics  characterized  by  an 
ISS  property.  The  resulting  closed-loop  system  is  depicted  in  Figure  1.7,  where  the 
unmodeled  dynamics  and  the  observer  error  system  blocks  have  finite  ISS  gains.  The 
remaining  task  for  the  control  design  is  to  assign  the  gains  for  the  plant  with  respect 
to  both  observer  error  and  unmodeled  dynamics  inputs  so  as  to  render  the  gains 
of  both  inner  and  outer  loops  less  than  one.  Then,  GAS  is  guaranteed  via  the  ISS 
small-gain  theorem  of  Teel  et  al.  [39,  93]. 


Figure  1.7:  Closed-loop  system  with  observer  feedback. 

In  Chapter  7  we  illustrate  the  design  on  an  axial  compressor  model  due  to  Moore 
and  Greitzer  [65] ,  which  has  been  the  starting  point  for  numerous  jet  engine  control 
studies.  We  study  the  following  single-mode  approximation  of  the  Moore-Greitzer 
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PDE  model 


4>  —  -1p+-(p+-  —  -((p+l)3—3((p  +  l)R 
Ip  =  -jp(<p  +  l  -u) 

R  =  erR(-2<p  -  <j)2  -  R) ,  R(0)  >  0 , 

where  (p  and  ip  are  the  deviations  of  the  mass  flow  and  the  pressure  rise  from  their  set 
points,  the  control  input  u  is  the  flow  through  the  throttle,  and,  a  and  f3  are  positive 
constants. 

Krstic  et  al.  designed  a  state  feedback  GAS  control  law  in  [53,  Section  2.4], 
and  later  replaced  it  by  a  design  using  <p  and  ip  in  [52].  With  a  high-gain  observer, 
Isidori  [27,  Section  12.7],  and  Maggiore  and  Passino  [58],  obtained  a  semiglobal  result 
using  the  measurement  of  ip  alone.  With  y  =  ip,  our  design  achieves  GAS.  The  exact 
observer  cannot  be  designed  because  of  the  nonlinearities  c pR  and  (p2R.  However, 
the  (cp,  VO-subsystem  (7.31), (7.32)  contains  the  nondecreasing  nonlinearity  (cp  +  l)3, 
and  we  exploit  this  fact  by  treating  the  i?-subsystem  as  unmodeled  dynamics.  The 
small-gain  assignment  is  then  achieved  via  an  observer-backstepping  design. 

1.3  Notation  and  Terminology 

A  function  /  :  lRn  — >  lRq  is  Ck  if  its  partial  derivatives  exist  and  are  continuous 
up  to  order  k,  1  <  k  <  oo.  A  C°  function  is  continuous.  A  C°°  function  is  smooth, 
that  is,  it  has  continuous  partial  derivatives  of  any  order.  The  same  notation  is  used 
for  vector  fields. 

A  function  a  :  1R> o  — >  M>q  is  said  to  be  class -1C  if  it  is  continuous,  increasing, 
and  <t(0)  =  0.  It  is  called  class-ZC^  if,  in  addition,  limt^+oocrpb)  =  +oo. 

A  function  f3  :  1R> o  x  ]R> o  — >  IR>o  is  said  to  be  class-/C£  if  for  each  t  £  IR>o, 
/?(•,£)  is  class -K,  and  for  each  s  £  1R> o,  /3(s,  •)  is  decreasing  and  limt^+oo/3(s,t)  =  0. 

A  locally  Lipschitz  function  (p{-)  :  1R  —>  [—A,  A]  is  called  a  saturation  function 
with  saturation  level  A  >  0,  if  <p{x)  =  x  when  |x|  <  A/2,  and 

A/2  <  sgn(x)<p(x )  <  min{|x|,A}  (1.56) 

when  |a;|  >  A/2.  This  definition  incorporates  the  standard  saturation  function  < p{x)  = 
sgn(x )  min{|x'|,  A}. 

Given  a  vector  field  /  :  Mn  — >  Mn  and  a  differentiable  scalar  function  A  :  Mn  — > 
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Table  1.1:  List  of  acronyms. 


CLF 

control  Lyapunov  function 

LES 

local  exponential  stability 

GAS 

global  asymptotic  stability 

LMI 

linear  matrix  inequality 

GES 

global  exponential  stability 

PR 

positive  real 

IOS 

input-to-output  stability 

SPR 

strictly  positive  real 

ISS 

input-to-state  stability 

]R,  LfX  denotes  the  directional  derivative  of  A  along  /: 

f)\ 

LfKx )  =  ^/O).  (L57) 

A  smooth,  positive  definite  and  radially  unbounded  function  V(x)  is  called  a 
control  Lyapunov  function  (CLF)  for  the  system  x  =  f(x)  +  g{x)u  if,  for  all  r/0, 

LgV(x)  =  0  =$■  LfV(x)  <  0.  (1.58) 

Throughout  the  dissertation,  |  •  |  denotes  the  Euclidean  norm.  We  say  that  a 
measurable  function  u(t)  is  locally  bounded  if,  for  all  T  >  0,  supter0)T]  |u(i)|  <  oo, 
where  sup  stands  for  the  essential  supremum.  For  a  locally  bounded  u(t),  we  define 
the  Coo  and  asymptotic  norms  as 

IMloo  :=  sup  KOI,  ||u||a  :=  limsup  \u(t)\,  (1.59) 

t>  0  t— >oo 

respectively. 

Table  1.1  contains  a  list  of  acronyms  used  throughout  the  dissertation. 
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Part  I 

Robust  Redesigns  Against 
Unmodeled  Dynamics 
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Chapter  2 

Robustification  of  Backstepping 


In  this  chapter  we  redesign  backstepping  schemes  such  as  those  in  Kanellakopou- 


los  et  al. 

[41],  and  Krstic  et  al.  [53],  to  robustify  them  against  input  unmodeled 

dynamics. 

We  consider  systems  of  the  form 

X  =  F(X)  +  G(X)x  i 

xi  =  fi(X,xi)  +  gi(X,xi)x2 

X2  =  f2(X,xi,x2)  +g2(X,xi,  x2)x3 

(2.1) 

Xn  =  fn(X,x)  +gn(X,x)v 

£  =  q(£,u) 

v  =  p(€,u). 

(2.2) 

where  \gi(X, Xi)\  >  go  >  0,  V(X, ...,  Xi)  6  lRr+l,  i  =  1,--*  ,n.  The  ^-subsystem 
(2.2)  with  input  u  6  M  and  output  v  £  M  represents  unmodeled  dynamics,  that  is, 
(2.1)  with  v  =  u  is  the  nominal  system.  When  u  =  0,  the  system  (2.1), (2. 2)  has  an 
equilibrium  at  zero,  whose  stability  properties  are  to  be  analyzed. 

For  the  X-subsystem  with  x\  viewed  as  a  virtual  control  input,  a  CLF  Vq(X) 
and  a  control  law  A°(X),  A°(0)  =  0,  are  known  such  that,  for  all  X  ^  0, 

Lf+ga oV0(X)  :=  -^(F  +  GA°)  =  -U0(X)  <0.  (2.3) 

With  the  knowledge  of  Vq(X)  and  A°(X),  backstepping  can  be  applied  to  guarantee 
GAS  for  the  nominal  system.  However,  a  GAS  control  law  for  the  nominal  system  does 
not  guarantee  GAS  in  the  presence  of  unmodeled  dynamics.  To  robustify  backstepping 
designs,  we  propose  two  redesign  methods:  passivation  and  truncated  passivation.  In 
the  first  redesign  we  use  the  results  of  Harnzi  and  Praly  [24],  and  ensure  GAS  via  the 
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passivity  properties  of  the  closed-loop  system.  In  the  second  redesign  we  passivate 
the  A-subsystem,  and  proceed  with  backstepping.  For  both  redesigns  the  unmodeled 
dynamics  subsystem  (2.2)  is  restricted  to  be  minimum  phase  and  relative  degree  zero. 

In  Section  2.1  we  review  the  two  main  versions  of  backstepping:  cancelation 
backstepping  and  LqV -backstepping.  In  Sections  2.2  and  2.3,  we  present  the  two 
redesigns.  Vo(A),  A°(A),  and  the  system  (2.1),  (2.2)  are  assumed  to  be  sufficiently 
smooth. 

2.1  Cancelation  and  LcF-Backstepping 

Backstepping  design  starts  with  the  virtual  control  law  A0  (A)  designed  for  the  A- 
subsystem  as  in  (2.3).  By  adding  —  i]qLgVo(X)  to  A°(A)  we  obtain  the  virtual  control 
law 

A(A)  =  A0 (A)  -  mLGV{ 3(A),  %  >  0,  (2.4) 

which,  when  r/o  >  0,  increases  the  negativity  of  Lf+G\Vq{X).  With  the  error  variable 
yi  :=  xi  —  A(A),  the  (A, xi)-subsystem  driven  by  X2  is 

A  =  (F  +  GA)(X)  +  G(X)y1  (2.5) 

m  =  h(X,xi)  -  A(A,yi)  +g1(X,xi)x2,  (2.6) 

where  A(A,  yi)  is  explicitly  known  from  (2.4)  and  (2.5). 

Step  1.  To  find  a  virtual  control  law  a^A, xi)  for  the  (A,  xi)-subsystem,  we 
introduce  the  CLF 

V1(X,x1):=V0(X)  +  ^-yl  (2.7) 

2/il 

where  >  0  is  to  be  specified.  Its  time  derivative  along  (2.5)  and  (2.6)  is 
Vi(X,Xl)  <  -Uo(X)  -  V0(LGV0(X))2  +  LGV0(X)yi 

+~y i  xi)  -  A(A,  yi)  +  gi{X,  xi)x2)  •  (2.8) 

Our  aim  is  to  select  x2  =  a5(A, xi)  that  renders  V\  negative  definite.  Using 

-V0(LGV0(X))2  +  LGV0(X)yi  <  -^-yl  (2.9) 

4r/0 

we  obtain 

Vl(X,xi)  <  -C/0(A)  +  ±-y\  +  ±-yi  (/i(A,xi)  -  A(A,  y{)  +  9l(X,  xi)x2)  ,  (2.10) 
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which  suggest  the  virtual  control  law 


a?(X,xi)  = 


gi(X,x i)  L 


-km  -  fi(X,xi)  +  A(X,yi) 


k\  >  0. 


(2.11) 


Then,  X2  =  a\(X,xi)  yields 


Vi(x,Xl)  <  -u0(x)  +  {^---)yl 

V4??0  Mi, 


(2.12) 


Setting  0  <  /.ii  <  4/cir/tb  we  guarantee  that  V\(X,x\)  is  negative  definite. 

We  refer  to  this  design  as  LqV -backstepping  because  the  sign-indefinite  term 
LqVq(X)di  in  (2.8)  is  dominated  by  adding  —tiqLgVq(X)  to  the  previous  virtual 
control  law  A°(X).  In  cancelation  backstepping  we  keep  the  control  law  A(X)  = 
A°(X)  by  setting  rjo  =  0  in  (2.4).  In  this  case,  the  LgVq(X)di  term  in  (2.8)  must  be 
canceled  by  a\(X,x\).  A  virtual  control  law  incorporating  both  types  of  backstepping 
is 

<A[X,x{)  =  — 4 - r[-km  -  Ai LgVo(X)  -  fi(X,Xl)  +  A(X,yi)},  (2.13) 

gi(X,xi) 

where,  for  cancelation  backstepping,  Ai  >  0  and  ji\  =  Ai  in  (2.7).  If  r/o  >  0,  we 
recover  the  LcW-backstepping  control  law  (2.11)  with  Ai  =  0. 

For  the  ( X ,  y\  )-subsystem  (2. 5), (2. 6)  driven  by  X2,  the  input  vector  field  is  given 
by  [0 ,  gi(X,x\)]T .  Differentiating  (2.7)  along  this  vector  field,  the  counterpart  of 
L( jVo(Al)  for  the  A-subsystem  is  gi(X,xi)yi.  Then,  as  in  (2.4),  we  introduce  the 
virtual  control 


ai(X,xi)  =  al(X,x{) -rngi(X,xi)yi,  gi>0,  (2.14) 

which  enables  us  to  avoid  a  cancelation  in  the  next  step,  and  define  the  error  variable 
y2  '■=  x2  -  «i(A,xi). 

Step  i.  (i  =  2,  •  •  •  ,  n)  We  take  the  CLF 

Vi(X,  xi, ...,  Xi)  =  Vo(X)  +  — — y\+  ...  +— — yf,  (2.15) 

Z/il  Z/Jsi 

where  m  >  0  is  to  be  specified.  With  ki  >  0,  calculations  similar  to  Step  1  yield  the 
virtual  control  law 

at}(X,...,Xi)  =  1 - -[-km  -  fi(X,...,Xi) 

gi{^,  Xi) 

+dj_i(A, ... ,yi )  -  \igi-i(X, ..., Xj_i)?/j_i],  (2.16) 
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where  A*  >  0  in  cancelation  backstepping,  while  ry_ i  >  0  and  A,  =  0  in  LqV- 
backstepping.  By  selecting  m  =  Aj//j_i  in  the  first,  and,  0  <  ^  <  Aktyt-i  in  the 
second  case,  V)(X,  x i,  is  rendered  negative  definite.  For  i  <  n,  we  define 

ai(X,  :=  a°(X,  ...,Xj)  -  rngi(X,  ...,Xi)yi,  ry  >  0, 

and  the  next  error  variable  j/j+i  :=  Xj+i  —  «i(X,  ...,Xi). 

This  procedure  results  in  the  CLF 

Vn{X,x)  =  Vo(X)  +  ±-y\  +  ...  +  (2.17) 

Z/il 

and  the  control  law 

U  —  CYn(X,  x)  —  — — :  r [  kn yn  /n(X,  x)  (2.18) 

5n(X,  x) 

+dn_i(X,  y i , . .. ,  yn—  1 )  \ifjn- 1  (X, ...,  xn_i )yn_i] , 

which  renders  Vn(X,  x\,  ...,xn)  negative  definite,  thus  achieving  GAS  for  the  nominal 
system. 

2.2  Passivation  Redesign 


In  this  section  we  consider  input  strictly  passive  unmodeled  dynamics,  charac¬ 
terized  by  a  constant  8  >  0  and  a  positive  definite,  radially  unbounded  function  £(£) 
such  that 


S(£)  <  — 8u 2  +  vu. 

(2.19) 

For  linear  unmodeled  dynamics  v  =  A (s)u,  this  means 

Re{A(ju)}  >  8,  Vcu  G  JR. 

(2.20) 

In  compact  notation,  the  system  (2.1), (2. 2)  is 

x  =  $(x)  +  r(x)^ 

(2.21) 

cTry. 

II 

(2.22) 

u  =  p(£,u), 

where 


X  :=  [ XT ,  xT]T 

$(X)  :=  [F(X)  +  G(X)x\,  -  ■  ■  ,fn(X,x)]T 
T(x)  :=  [0,  ■  ■  ■  ,gn(X,x))]T. 

The  passivation  redesign  makes  use  of  the  following  lemma: 
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Lemma  2.1  Consider  the  system  (2. 21), (2. 22),  and  suppose  that  the  £- subsystem 
satisfies  (2.19)  with  6  >  0;  and  is  GAS  when  u  =  0.  If  there  exists  a  positive  definite 
and  radially  unbounded  function  V(x)  such  that,  for  all  x  £  lRr+n  —  {0}, 

L*V(X)  <  (LrV(X))2,  (2-23) 

then  the  control  law 

u  =  -kLrV(x),  k  >  ^  (2.24) 

guarantees  GAS. 

Proof:  With  the  control  law  (2.24),  we  view  the  closed- loop  system  as  the  feedback 
interconnection  of  the  subsystems  (2.21)  and  (2.22),  and  examine  the  passivity  prop¬ 
erties  of  each  subsystem.  For  (2.21),  we  use  kV(x)  as  a  storage  function,  and  denote 
P(x)  :=  (LpU(x))2  —  L$V(x)  >  0,  Vx  /  0.  Substituting  (2.24),  we  obtain 

kV (x)  =  ^u2  -uv-  k(3(x).  (2.25) 

Adding  (2.19)  and  (2.25),  we  get 

kV{x)  +  5(0  <  -(^  -  l)u2  -  k(3(x)  <  -k(3(X),  (2-26) 

which  establishes  global  stability  of  the  closed-loop  system  (2. 21), (2. 22)  with  (2.24). 
It  also  follows  from  (2.26)  that  x  ~ ^  0.  Since  u  =  0  when  x  =  0,  and  the  ^-subsystem 
is  GAS,  we  conclude  from  LaSalle’s  invariance  principle  that  £  — ►  0,  and  hence,  the 
closed-loop  system  is  GAS.  □ 


The  existence  of  a  function  satisfying  (2.23)  for  general  nonlinear  systems  has 
been  shown  by  Hamzi  and  Praly  [24]  to  be  equivalent  to  the  existence  of  a  CLF  V (x) 
and  a  constant  l  >  0  such  that 


lim  sup 

x->o 


L*V(x) 

(LrV (x))  2 


<  l. 


(2.27) 


Under  this  local  condition,  a  continuous,  positive  definite  function  9  :  M>o  — >  1R> o 
exists  such  that 


0(V(x))  > 


L*V(X) 

(£rW(x))2’ 


Vx  e  Mr+n  -  {0}, 


(2.28) 


Unit 


9(s)d-s  =  Too. 


(2.29) 
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Then,  it  can  be  verified  that 


rV(x) 

V(x)=  /  0(s)ds 


Jo 

is  a  positive  definite  and  radially  unbounded  function  which  satisfies  (2.23). 
therefore  useful  to  determine  when  a  CLF  V(x)  will  satisfy  (2.27). 


(2.30) 
It  is 


Lemma  2.2  Let  P  :  = 


d2V 
dx 2 


j  X=0 


If  kxT Px  is  a  CLF  for  the  Jacobian  linearization 


of  (2.21),  then  V(x)  satisfies  (2.27). 


□ 


We  now  show  when  this  will  be  the  case  for  Vn(x)  constructed  by  backstepping. 

Theorem  2.1  Consider  the  system  (2.21)-(2.22),  and  suppose  that  the  £ -subsystem 
satisfies  (2.19)  with  5  >  0,  and  is  GAS  when  u  =  0.  If 


Pa  = 


d2V0 

dx2 


and  Qo  = 


d2U0 
dX 2 


(2.31) 


J  X=0  L  J  X=0 

in  (2.3)  are  both  positive  definite,  then  there  exists  a  positive  definite  function  V(x) 
such  that  the  control  law 


u  =  -kLrV(x),  &  > 

renders  the  closed-loop  system  GAS. 


(2.32) 


Proof:  With  the  procedure  of  Section  2.1  applied  to  (2.21),  the  quadratic  part  of 
Vn(x )  m  (2.17) 

1 

rT  i 


-XTP0X  +  —  y\  + 

2  2/ii yi 


+ 


1 


2/ir 


Vn 


is  a  CLF  for  the  Jacobian  linearization  of  (2.21).  It  follows  from  Lemma  2.2  and 
the  preceding  discussion  that  there  exists  a  positive  definite  and  radially  unbounded 
V(x)  =  "  6{s)ds  which  satisfies  (2.23).  Then,  by  Lemma  2.1,  (2.32)  achieves 

GAS.  □ 


Substituting  L^V  =  0(Vn)G-gn(X,x)yn  in  (2.32),  the  redesigned  control  law  is 
u  =  —  9(Vn)gn(X,x)yn,  k>\,  (2.33) 

hn  0 

where  6  :  lR>o  — >  IR> o  is  a  continuous,  positive  definite  function  that  satisfies  (2.28) 
and  (2.29).  It  is  important  to  note  that  this  redesign  does  not  require  detailed  in¬ 
formation  about  the  unmodeled  dynamics,  only  a  lower  bound  on  5  is  assumed  to  be 
known. 
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Example  2.1  Consider  the  system  with  linear  unmodeled  dynamics  A(s): 

X  = 
x  = 

v  = 

Since  Re{A(jco)}  >  1,  A(s)  satisfies  (2.19)  with  5  =  1.  Using  the  virtual  control 
A0  (A)  =  —X  —  X 3  for  the  X -subsystem,  we  see  that  (2.3)  is  satisfied  with  Vq(X)  = 
±X2,  U0(X)  =  X2,  that  is,  Pq  =  1 ,  Qo  =  2  in  (2.31).  We  define  the  error  variable 
y  =  x  +  X  +  X 3  and  set  =  1  in  (2.17),  that  is,  a: )  =  \X2  +  \y2 .  In  the 

(X,y)-  coordinates  the  system  is 


X3  +  x 


(2.34) 


12(s  +  35)(s  +  20) 
7(s  +  40)  (s  +  30) 


u  =  A  {s)u, 


X  =  -X  +  y 

y  =  (1  +  3A2)(— X  +  y)  +  v,  (2.35) 


and,  hence, 

L^V2(X,y)  =  -X2  -  3X3y  +  (1  +  3A2)y2 

LvV2{X,y)  =  y.  (2.36) 


We  now  need  to  find  a  continuous  function  0{V2)  such  that 

— 3A  1 


m)>  =  i+3x* +x°- 


(LrV,(X,y))? 


y  y 


Using  the  inequality 


r  1  r 

-  -  —  <  J,  Vy&m, 


y  y 

we  obtain,  upon  completion  of  the  squares, 

ii^wil+3X2+ix‘<3+im{X-y?- 

Choosing  9(s )  =  3  +  18s2,  k  =  \  =  1,  and  substituting  in  (2.33),  we  obtain  the 
redesigned  control  law 


u=-3y-^(X2  +  y2)2y, 


which  achieves  GAS. 


(2.37) 

□ 
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2.3  Truncated  Passivation  Redesign 

For  higher  order  systems  the  task  of  finding  the  function  0{  )  may  be  cumber¬ 
some.  We  now  circumvent  this  difficulty  by  a  ‘truncated’  design  in  which  we  passivate 
the  X-subsystem  only,  and  then  apply  backstepping  to  the  redesigned  virtual  control 
law. 

For  the  Jacobian  linearization  of  the  X-subsystem,  ^ XtPqX  is  a  CLF,  provided 
Po  and  Q o  defined  in  (2.31)  are  positive  definite.  Therefore,  we  conclude  from  Lemma 
2.2  that  there  exists  a  positive  definite  and  radially  unbounded  function  Vq(X)  that 
satisfies 

LfVo(X)  <  LgV0(X)2,  VX  €  Mr  -  {0}.  (2.38) 

Then,  the  virtual  control  law 

A(X)  =  -k0LGV0(X),  k0  >  1,  (2.39) 

renders  Vo(X)  =  LpV q(X)  +  LgVq(X)A(X)  negative  definite. 

To  make  the  main  features  of  the  truncated  passivation  redesign  more  apparent, 
we  present  it  for  the  special  case  of  (2.1)  in  which  the  X-subsystem  is  augmented  by 
a  chain  of  n  integrators,  and  the  unmodeled  dynamics  are  linear: 


X 

=  F(X)  +  G(X)x  i 

Xl 

=  X2 

X2 

=  x3 

(2.40) 

=  A  (s)u. 

For  this  subclass,  the  sequence  of  the  virtual  control  laws  obtained  via  LcW-backstepping 
is 

aq 

OL2 

cxn  —  kn(xn  crn_i)  T  Qn_ i,  (2.41) 

where  ki>  0,  i  =  1, 2,  •  ■  ■  ,  n,  and  A  is  as  in  (2.39). 


=  —k\{x\  —  A)  +  A 
=  -k2(x2  -  a±)  +  oq 
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To  represent  the  closed-loop  system  as  the  feedback  interconnection  of  a  linear 
and  a  nonlinear  subsystem,  we  write  the  xn-equation  of  (2.40)  with  u  =  an  in  the 
Laplace  domain  as 


which  yields 


SXn  —  A(s)[  k'ri  (xn  Otn—  l)  T  SOtn—  l] , 


( S  +  kn)A{s) 
s  +  knA(s) 


a ■ 


i-i  — :  Ai(s)a„_i. 


(2.42) 


(2.43) 


Next,  we  substitute  xn  =  sxn-\  and  an_i  =  —kn- \xn-\  +  (s  +  kn-i)an-2  in  (2.43), 
and  obtain 


(s  +  kn-i)Ai(s)  r 

Xn-l  =  - - - ~  OLn-2  =:  A2(s)an_2. 

S  +  kn-iA^S) 


(2.44) 


Proceeding  recursively,  we  get 


xi  =  An(s)A(X), 


where  A n(s)  is  obtained  from 


A0(s) 


A  i(s) 


A  (s) 

(s  +  kn—i- |_r )  Aj_i(.s) 

s  +  kn—i-\-\ Aj_x  (s) 


Thus,  the  closed-loop  system  (2.40)-(2.41)  is 

X  =  F{X)  +  G{X)x  i 
X!  =  An(s)A(X). 


(2.45) 


(2.46) 


(2.47) 


Because  A(X)  in  (2.39)  is  in  Ley-form,  we  can  use  Lemma  2.1  to  prove  the  following 
stability  margin: 

Lemma  2.3  The  closed-loop  system  (2-4  7)  is  GAS  for  all  stable  An(jcu)  that  satisfy 

Re{An(ju)}  >  Moj  €  JR.  (2.48) 

ko 

□ 


Our  hnal  result  shows  how  the  design  parameters  are  to  be  selected  to  satisfy 


(2.48). 
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Theorem  2.2  Consider  the  system  (2-40),  with  a  controller  designed  according  to 
(2.39)  and  (2.)1).  Suppose  that  A(s)  is  minimum  phase  and  relative  degree  zero, 
with  high-frequency  gain 

h  :  =  lim  A (ju)  >  0. 

id — >00 

If  we  select  fco  >  max{l,  and  ki  >  0,  i  =  1,  ■  ■  ■  ,n  —  1,  then  there  exists  k*  >  0 
such  that  kn  >  k *  guarantees  GAS  for  the  closed-loop  system. 


Proof:  From  (2.46),  An_i(s)  is  minimum  phase,  relative  degree  zero  and  its  high- 
frequency  gain  is  h  >  0.  This  implies  that 


A„(s)  =  <»  +  HA-lW 

<  +  tiViW 


(2.49) 


is  stable  for  sufficiently  large  k\  >  0,  as  verified  from  a  root-locus  argument.  Next,  it 
can  be  shown  from  (2.49)  that,  as  k\  is  increased,  the  Nyquist  plot  of  A n(s)  converges 
to  that  of  .  which  is  a  circle  that  intersects  the  real  axis  at  1  and  h.  Since 

^  <  min{l,  /;,},  (2.48)  is  satisfied  if  k\  is  selected  sufficiently  large.  Thus,  GAS 
follows  from  Lemma  2.3.  □ 


Example  2.2  For  the  system  of  Example  2.1,  we  now  perform  the  passivation  re¬ 
design  for  the  X -subsystem  only,  and  then  apply  backstepping.  A  redesigned  control 
law  for  the  X-subsystem  is  A(A)  =  —2 L0A3.  Then,  the  LqV -backstepping  control  law 

u  =  —k\{x  +  2k$X 3)  —  6koX2(X3  +  x),  (2.50) 

achieves  GAS  for  the  nominal  model.  To  guarantee  GAS  with  A(s),  we  note  that 
h  =  1,  and  choose  h®  >  1.  Then,  by  Theorem  2.2,  the  closed-loop  system  (2. 34),  (2. 50) 
is  GAS  for  sufficiently  large  k\  >  0.  □ 


2.4  Summary 

We  have  presented  two  passivation  redesigns  of  backstepping  which  achieve 
global  asymptotic  stability  for  a  class  of  minimum  phase  unmodeled  dynamics  with 
relative  degree  zero.  The  redesigns  do  not  require  detailed  information  about  the 
unmodeled  dynamics:  the  class  of  admissible  unmodeled  dynamics  are  characterized 
by  their  passivity  properties  in  the  first  redesign,  and  their  high-frequency  gains  in 
the  second  redesign.  Our  analysis  provides  insight  into  the  robustness  properties  of 
backstepping  designs. 
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Chapter  3 

Dynamic  Nonlinear  Damping  Redesign 


The  redesigns  presented  so  far  restrict  the  unmodeled  dynamics  by  small-gain 
or  passivity  conditions.  The  dynamic  nonlinear  damping  redesign  presented  in  this 
chapter  removes  these  restrictions.  Instead,  the  main  restriction  is  that  the  unmod¬ 
eled  dynamics  subsystem  be  relative  degree  zero  and  minimum  phase.  For  nonlinear 
unmodeled  dynamics,  the  minimum  phase  requirement  is  replaced  by  a  robust  stabil¬ 
ity  property  of  the  zero  dynamics.  Our  control  law  employs  a  dynamic  normalization 
signal  to  counteract  the  destabilizing  effect  of  the  unmodeled  dynamics. 

A  practical  feature  of  this  redesign  is  that  it  does  not  require  detailed  information 
about  the  unmodeled  dynamics  subsystem,  only  the  rate  of  exponential  convergence 
and  the  sign  of  the  high-frequency  gain  are  assumed  to  be  known.  The  closed-loop 
solutions  are  bounded,  and  converge  to  a  compact  set  which  can  be  made  arbitrarily 
small  by  increasing  the  controller  gain.  If  the  Jacobian  linearization  of  the  zero  dy¬ 
namics  is  asymptotically  stable,  then  the  redesign  recovers  global  asymptotic  stability 
and  local  exponential  stability  (LES). 

In  Section  3.1  we  introduce  the  class  of  systems  to  be  studied,  and  characterize 
the  admissible  unmodeled  dynamics.  In  Section  3.2,  we  present  our  redesign  and 
illustrate  it  on  an  analytical  example.  The  proofs  are  given  in  Section  3.3. 

3.1  Problem  Statement 

We  consider  the  system 

X 
x 

s 
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F(X,x) 

(3.1) 

f{X,  x)+g(X,  x)v 

(3.2) 

A(0  +  Bu 

(3.3) 

c(£)  +  du, 

(3.4) 

v 


in  which  X  £  IRn,  x  £  1R,  f  £  JRm,  and  |g(X,  x)|  >  go  >  0  for  all  (X,  x)  £  J?n+ 1 .  The 
^-subsystem  with  input  u  £  IR  and  output  v  £  IR  represents  unmodeled  dynamics, 
that  is,  the  (X,  x)-subsystem  with  v  =  u  is  the  nominal  system.  It  is  assumed  that  all 
functions  in  (3.1)-(3.4)  are  Cl ,  and  _F(-,-),  /(-,-),  A(-)  and  c(-)  vanish  at  zero.  The 
stability  properties  analyzed  are  with  respect  to  zero,  which  is  an  equilibrium  for  the 
system  (3.1)-(3.4)  when  u  =  0. 

The  main  restriction  on  the  nominal  system  is  that  the  X-subsystem  be  globally 
stabilizable  with  x  viewed  as  a  virtual  control  input. 

Assumption  3.1  There  exists  a  C 1  function  y(X)  such  that 

X  =  F(X,/x(X))  (3.5) 


is  globally  asymptotically  stable. 

We  will  consider  nominal  control  laws  that  guarantee  an  ISS  property  for  the  X- 
subsystem.  Our  redesign  will  render  such  control  laws  robust  against  the  destabilizing 
effect  of  the  unmodeled  dynamics. 

For  nominal  control  laws  we  employ  a  backstepping  procedure  which  relies  on 
the  following  fact  proved  by  Sontag  [87]  (see  also  Isidori  [27,  Theorem  10.4.3]): 

Proposition  3.1  If  /x(X)  is  as  in  Assumption  3.1,  then  there  exists  a  C 1  function 
0  <  /3(X)  <  1  that  renders 

X  =  F(X,ii(X)+(3(X)y)  (3.6) 

ISS  with  input  y.  □ 

To  proceed  with  the  backstepping  design,  we  note  that  the  variable  y  :=  (x  — 
p(X)) / P(X)  is  governed  by 

V  =  /(X,  x)  +  g(X,  x)v,  (3.7) 

where  \g(X,x)\  =  >  go-  Then,  for  the  nominal  system  (3.1)-(3.2)  with  v  =  u, 

the  control  law 

u  =  a(X,  x)  =  ^  (-/(X,  x)  —  ky'j  ,  k  >  0 ,  (3.8) 

results  in  y  =  — ky ,  and  guarantees  GAS  because  of  the  ISS  property  of  (3.6). 

The  admissible  unmodeled  dynamics  are  characterized  by  the  following  assump¬ 
tions: 
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Assumption  3.2  The  unmodeled  dynamics  subsystem  (3.3)- (8.4 )  has  relative  degree 
zero,  that  is  d  /  0. 


Assumption  3.3  There  exists  a  constant  c  >  0  such  that  |c(£)|  <  c|£|. 
Assumption  3.4  There  exists  a  C 1  Lyapunov  function  V(£)  such  that 


Hei2  <  no  <  Hei2 

<  -2^(0 


dV 

~df 


<  v3\t\, 


where  5,  v\,  V2,  v3  are  positive  constants. 


(3.9) 

(3.10) 

(3.11) 
□ 


The  constant  c  and  the  Lyapunov  function  G(£)  serve  only  to  define  the  class 
of  unmodeled  dynamics.  Their  knowledge  is  not  required  for  the  redesign.  By  As¬ 
sumption  3.4,  if  u  =  0,  then  the  unmodeled  dynamics  (3.3)  are  globally  exponentially 
stable  (GES)  with  the  convergence  rate  5.  The  converse  of  this  statement  is  also  true 
if  A(£)  is  globally  Lipschitz,  because,  if  the  unmodeled  dynamics  subsystem  is  GES, 
then  there  exists1  a  G(£)  as  in  Assumption  3.4. 

The  final  assumption  requires  a  robust  stability  property  for  the  zero  dynamics 
of  the  subsystem  (3.3)-(3.4),  that  is 

i  =  A(z)  -  c(z )  =:  A°(z).  (3.12) 

Assumption  3.5  The  zero  dynamics  subsystem  (3.12)  disturbed  by  d\  and  c?2 

z  =  A°(z  +  d\)  +  g?2  (3.13) 


is  ISS  with  input  (^1,^2)-  LI 

For  linear  unmodeled  dynamics  A(£)  =  A£,  c(£)  =  C(,  Assumptions  3.4  and  3.5 
are  satisfied  if  the  matrices  A  and  A0  :=  A  —  ^BC  are  both  Hurwitz.  This  means 
that  all  relative  degree  zero,  stable,  and  minimum  phase  linear  unmodeled  dynamics 
are  admissible. 

1See  Khalil  [44,  Theorem  3.12],  and  its  recent  extension  by  Corless  and  Glielmo  [12]  which  guar¬ 
antees  the  same  convergence  rate  for  E(£(f))  as  for  £(t). 
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3.2  Dynamic  Nonlinear  Damping  Redesign 

Using  a(X,  x)  as  in  (3.8),  the  control  law  for  system  (3.1)-(3.4)  is 

u  =  sgn(d)[a(X,  x)  —  k  (1  +  \m\  +  \a(X,  .x)|)  g(X,  x)y\ ,  n>  0  (3.14) 

m  =  —  5m  +  \u\,  (3.15) 

which  only  requires  the  knowledge  of  6  and  the  sign  of  d.  Applying  (3.14)  to  (3.1)- 
(3.4),  we  obtain  the  closed-loop  system 

X  =  F(X,n(X)+/3(X)y)  (3.16) 

V  =  f{X,x)  (3.17) 

+3(^,a;){c(0  +  \d\  [a{X,x)  -  re(l  +  \m\  +  \a(X,x)\)  g(X,x)y]} 

£  =  A(0  +  Bu.  (3.18) 

Adding  and  subtracting  (3.8),  we  rewrite  (3.17)  as 

y  =  -ky  +  g(X,x)[c(£)  +  (|d|  -  l)a{X,x)  -  K\d\(l  +  \m\  +  \a(X,x)\)  g(X,x)y\. 

(3.19) 

The  three  robustification  terms  in  —  k(1  +  \m\  + 1 a(X,  x)\)  g(X,  x)y  counteract  the  un¬ 
modeled  terms  c(£)  and  (|d|  —  l)a(X,  x)  in  (3.19).  In  particular,  —  n  \ a(X,  x)|  g(X,  x)y 
counteracts  (|d|  —  1  )a(X,x),  and  can  be  dropped  if  |ci|  =  1.  The  remaining  terms 
— ng(X ,  x)y  and  —  k  \m\  g(X,  x)y  counteract  c(^),  as  we  illustrate  with  the  help  of  the 
following  lemma: 

Lemma  3.1  Consider  equations  (3.15)  and  (3.18),  and  suppose  u(t )  €E  £oo[0,T].  If 
Assumption  3.f  holds,  then  there  exist  constants  0\,  62  >  0  such  that,  for  allt  £  [0,T], 

\m\  <  Oi  (|£(0)|  +  |m(0)|)  e~5t  +  e2\m(t)\.  (3.20) 

□ 

The  proof  is  given  in  Section  3.3.  For  every  time  interval  [0,  T]  in  which  the  closed-loop 
solutions  exist,  Lemma  3.1  implies  that 

|c(£(i))|  <  c|£(t)|  <  o0i(|C(O)|  +  |m.(0)|)e_5f  +  c62\m(t)\.  (3.21) 

Thus,  —ng(X,  x)y  and  —  k  \m\  g(X,  x)y  counteract  c0i(|^(O)|+|m(O)|)e_5t  and  c02\m(t)\, 
respectively. 

Our  redesign  guarantees  boundedness  of  the  states  (X,x,£)  and  the  signal  m 
for  any  k  >  0.  Moreover,  the  states  (X,x,£)  converge  to  a  compact  set  around  the 
origin,  which  can  be  made  arbitrarily  small  by  increasing  k. 
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Theorem  3.1  Consider  the  system  (3.1)-(3-4),  and  suppose  Assumptions  3. 2-3. 5 
hold.  Ifa(X,x )  is  as  in  (3.8),  then  the  control  law  (3.14)-(3.15)  guarantees  that  the 
closed-loop  solutions  are  bounded  for  allt  >  0.  Moreover,  there 

exists  a  class-fCL  function  /?(-,•)  with  the  property  that,  for  any  given  e  >  0,  we  can 
find  n*  such  that  for  all  k  >  n*, 

\(X(t),x(t),f,(t))\  <  max{/?(|(X(0),x(0),£(0),m(0))|,f) ,  e}  Vt  >  0 .  (3.22) 

If,  in  addition,  the  Jacobian  linearizations  of  F(X,  y(X))  in  (3.5)  andA°(z)  in  (3.12) 
are  asymptotically  stable  at  zero,  then  the  closed  loop  system  (3.1)-(3-4),(3.14)-(3.15) 
is  globally  asymptotically  stable  and  locally  exponentially  stable  for  sufficiently  large 

K.  □ 

The  proof  is  given  in  Section  3.3.  Inequality  (3.22)  indicates  that  the  state  variables  of 
the  plant  (3.1)-(3.4)  converge  to  a  ball  of  radius  e,  where  e  can  be  rendered  arbitrarily 
small  by  increasing  k.  It  is  important  to  note  that  increasing  k  does  not  cause  peaking 
in  the  transients  because  /?(•,  •)  is  independent  of  k. 

Example  3.1  The  redesign  for  the  system 

X  = 
x  = 

6  = 
i‘2  = 

V  = 

starts  with  the  design  of  the  nominal  control  law  a(X,x).  The  X -subsystem  satisfies 
Assumption  3.1  with  n(X)  =  —X  —  X3,  and  Proposition  3.1  holds  with  (i{X)  =  1. 
The  y-subsystem  is 

y  =  f{X,  x )  +  g(X,  x)v  =  x2  +  (1  +  3X2)(X3  +  x)  +  (1  +  X2)v  (3.25) 

and,  hence, 

a(X,  x)  =  +x2)  (~x2  ~  (1  +  3X2)(X3  +  x)  -  ky)  ,  k  >  0.  (3.26) 

The  remaining  task  is  to  show  that  the  unmodeled  dynamics  satisfy  Assumptions 
3.2-3. 5.  Assumptions  3.2  and  3.3  hold  because,  from  (3.24), 

d=  1  and  c(g)  =  gi  +  ^  ^  ,  (3.27) 

1  +  ^2 


X3  +  x 

x2  +  (1  +  X2)v 

-  £i  +  6 

— 2^2  +  u 
-6 


6  + 


10£23 


1  +  £9 


+  U 


(3.23) 


(3.24) 
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and  Assumption  3-4  is  verified  with  V(£)  =  ^(£2  +  C|)  an d  $  =  1-5.  To  verify 
Assumption  3.5  for  the  zero  dynamics 

„  I"  —2z\  -  zf  +  z2 

A°(z)  =  Z2-I0z|  ,  (3.28) 

L  "  2~2  +  -tt^t  J 

we  use  the  ISS  Lyapunov  function  V\  =  \zT z  and  note  from  (3.13)  that 

Vi  =  zt[A°(z  +  di)  -  A°(di)]  +  zT[A°(di)  +  d2\,  (3.29) 


which,  in  view  of 

t'1  \r)AQ' 

^4° (z  +  d\)  —  A0 (c?i)  =  /  — —  z  dX  (3.30) 

JO  L  Js=di+Az 

yields 

Vl  =  \{f0  zT(l^  +  )  zdX^+zT[A°(d1)  +  d2].  (3.31) 

v  \  /  s=di~\~Xz  ) 

It  follows  from  (3.28)  that 


'dA°~ 

dz 


T 

<  -21 


Vz  €  IR2 


(3.32) 


and,  hence, 

Vi  —zd'z  +  z^ \AP (d\)  +  d2]  <  —  +  —  | [t4°(c?i)  +  d2]\ ~ ,  (3.33) 

which  is  an  ISS  property  with  input  (d\,d2). 

Having  verified  that  Assumptions  3.2-3. 5  hold,  we  substitute  g(X,x)  =  1  +  X2, 
y  =  x  +  X  +  X3,  sgn(d)  =  1  and  5  =  1.5  in  (3.14)-(3.15),  and  obtain  the  control  law 


u  =  a(X,  x)  —  k  (1  +  \m\  +  \a(X,  x)|)  (1  +  X2)(x  +  X  +  X3) ,  k  >  0  (3.34) 
m  =  —  1.5m +  \u\,  (3.35) 


which  completes  the  redesign  of  the  nominal  control  law  a(X,x)  in  (3.26).  The  closed- 
loop  system  (3.23)-(3.24)  is  GAS  and  LES  for  sufficiently  large  k,  because  the  Jaco¬ 
bian  linearizations  of  A°(z)  and  F(X,  n(X))  =  —X  are  asymptotically  stable.  □ 
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3.3 


Proofs 


3.3.1  Proof  of  Lemma  3.1 


Starting  with  V (£)  as  in  Assumption  3.4,  we  use  the  Lyapunov  function  W (f)  := 
\JV (£)  which  is  not  differentiable  at  the  origin,  but  is  locally  Lipschitz  because  from 
(3.9)  and  (3.11), 


dW 

1  dV 

<  V 3 

di 

2VV  df 

2^/vl 

V£/0. 


(3.36) 


We  need  the  following  result  of  Teel  and  Praly  [97],  which  makes  use  of  the  generalized 
directional  derivative  of  Clarke  [11]. 


Proposition  3.2  Suppose  W(£)  is  locally  Lipschitz,  f(£,u)  and  a(f,u)  are  continu¬ 
ous,  and 

dw 

—/(£,«)  <a(£,«)  V«,V£gn,  (3.37) 

of 

where  ii  is  the  set  in  which  W  is  not  differentiable.  Let  u(t)  be  a  function  defined  on 
[0,  T],  and  let  f(t)  be  an  absolutely  continuous  function  satisfying  £  =  f(f(t),u(t))  on 
[0,  T\.  Then,  for  almost  all  t  €  [0,  T], 


W  <  a(£(t),u(t)). 


(3.38) 

□ 


To  evaluate  (3.37),  we  note  from  (3.10)  and  (3.36)  that,  for  all  £  0, 

dW  1  dV  uq 

—  m]  +  Bn)  =  5^  WO  +  Bn)  <  -SW  +  14  (3.39) 

Denoting  V4  :=  2^jt-||l?||,  we  conclude  from  Proposition  3.2  that,  for  almost  all  t  E 

[0 ,71, 

W  <  -5W  +  v4\u\  (3.40) 


and,  hence, 

W(f(t))  <  W{t(0))e~st  +  v4  [  e-s^\u(T)\dT  Vfe[0,T].  (3.41) 

Jo 

Substituting 

e-S{t-r)  \u(j^\dT  =  m(t)  —  m(0)e~5t,  (3.42) 
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obtained  from  the  solution  of  (3.15),  we  get 


W(£)  <  [W(£(0))  -  v4m(0)\e-st  +  v4m{t)  Vi  €  [0,  T\.  (3.43) 

Using  (3.9),  we  obtain 

VN£(*)I  ^  (\/^2|<f (0)|  +  V4\m(0)\)e~st  +  v4\m(t)\, 
from  which  (3.20)  follows  with  0\  =  -j=  rnaxjy^,  rq}  and  62  =  -^= 

3.3.2  Proof  of  Theorem  3.1 

We  divide  the  proof  into  four  parts.  First,  we  prove  that  in  (3.1)-(3.4)  there 
are  no  finite  escape  times,  and  that  X (t)  and  x(t)  are  bounded.  Next,  we  prove  that 
£(i)  is  also  bounded,  and  derive  (3.22).  In  the  third  part  we  prove  the  boundedness 
of  m[t).  Finally,  we  prove  GAS  and  LES  under  the  additional  assumption  that  the 
Jacobian  linearizations  of  A°(z)  and  F(X,  n(X))  are  asymptotically  stable. 

Part  1:  (X,x)  G  £oo,  (£,m)  G 

The  closed-loop  system  is  locally  Lipschitz  and,  hence,  solutions  exist  and  are 
unique.  To  prove  the  absence  of  finite  escape  times,  we  analyze  the  solutions  on  the 
compact  interval  [0,  T]  where  T  is  in  the  maximal  interval  of  existence,  and  show  that 
they  are  bounded  by  a  continuous  function  of  T  on  [0, 00).  From  (3.19), 

^  (^y2)  -  ~kV2  +  \§(x,x)y\  KOI  +  \g(X,x)y\  |(|d|  -  1)1  \a{X,x)\  (3.45) 
—n\d\  (1  +  \m\  +  |a(A,.x)|)  \g{X,x)y\2. 

Substituting  (3.21)  and  \g(X,x)y\  >  go\y\,  and  rearranging  terms,  we  obtain 

^  <  -ky2  -  K\d\\m\\g(X,x)y\  (^g0\y\- (3.46) 

-n\d\  \g(X,x)y\  (g0\y\  -  ^  (|£(0)|  +  M0)|)e“at^ 

-K\d\  \a(X,x)\  \g(X,x)y\  (g0\y\  -  • 

Defining 

m(0)|)e-a,  ^}, 
(3.47) 


r  :  = 


■  max 


90 


c02 

w 


-1)1 


Ml 


C(i)  :=  max 


c6 1 

n\d\go 


(M(0)|  + 


(3.44) 

□ 
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we  note  from  (3.46)  that 


(3.48) 


\y\  >  |C(*)| 


This  means  that 


d_ 

dt 


<  -ky 2. 


\y(t)\  <  max 


\yme~kt, 


sup  |C(r)| 

0  <T<t 


Vf  £  [0,  T], 


(3.49) 


which  is  an  ISS  property  with  input  £(t).  It  follows  from  (3.47)  that 

\((t)\  <  max  j|C(0)|e-5t,  ,  (3.50) 

that  is  ((t)  is  ISS  with  input  X  Moreover,  the  X-subsystem  (3.6)  is  ISS  with  input 
y.  Recalling  that  the  cascade  interconnection  of  ISS  systems  is  ISS  (see  Sontag  [86]), 
we  conclude  that  there  exists  a  class-/C  function  7o(-)  and  a  class-/C£  function  /3q(-,  •) 
such  that 


|(X(t),x(f))|  <  max  j/30(|(X(0),x(0),C(0))|,f),  70 


Vt  £  [0,T],  (3.51) 


This  gives  an  upper  bound  on  ( X (t),x(t))  which  is  independent  of  T.  Using  this  upper 
bound  in  (3.14)-(3.15),  we  can  find  a  constant  N  such  that  for  almost  all  t  £  [0,T], 

—  |ra|  <  _/V|ra|  +  IV,  (3.52) 

dt 

which  implies  that  \m(t)\  is  bounded  by  a  continuous  function  of  T  on  [0, 00).  Like¬ 
wise,  a  continuous  bound  can  be  derived  for  |£(t)|  using  (3.21).  Since  the  closed-loop 
signals  ( X ,  x,  £,  m)  are  bounded  by  a  continuous  function  of  T  on  [0, 00),  the  maximal 
interval  of  existence  is  infinite,  that  is  (X,x,£,m)  £  Moreover  (3.51)  holds  for 
all  t>  0,  which  proves  that  (X,x)  £  Coo- 


Part  2:  £  £  Coo 

Our  derivations  so  far  have  not  relied  on  Assumption  3.5.  We  now  employ 
Assumption  3.5  to  prove  that  £  £  Coo-  To  eliminate  u  from  (3.3),  we  use  (3.4)  and 
obtain 


i  =  A°{C)  +  B°v, 


(3.53) 


where  A°(£)  :=  A(£) 
that  the  variable 


^Rc(£)  and  B°  :=  Then,  it  follows  from  (3.2)  and  (3.53) 


rx 

z:=£-B°  / 

Jo 


da 

g(x,  a) 


(3.54) 
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is  governed  by 


z  =  A°  I  z  +  B' 


da 


g(X,a) 


-  B[ 


i  f(X,x) 

g{x,x ) 


(3.55) 


which  is  (3.13)  with  d\  =  B{)  g(x  a) 1  and  ^2  =  •  Since  d\  and  d2  are 

continuous  functions  of  (X,x),  and  vanish  at  zero,  it  follows  from  Assumption  3.5  and 
the  ISS  property  of  the  cascade  of  ISS  systems  that  there  exist  a  class-/C£  function 
/?i  (- ,  •)  and  a  class -1C  function  7i(-)  such  that,  for  all  t  >  0. 


\(X(t),x(t),£(t))\  <  max  |/3i(|(A(0),x(0),^(0),C(0))|,t),  71 


(3.56) 


This  proves  that  £  G  To  prove  that  (3.22)  holds,  we  need  to  eliminate  £(0)  from 
(3.56).  To  this  end,  we  let  R  >  0  and  observe  from  (3.47)  that  for  all  k>R , 

|C(0)|  <  max|^^(|C(0)|  +  |m(0)|),  ^  j  .  (3.57) 

Using  this  inequality  in  (3.56),  we  can  find  a  class-/C£  function  /3S(-,  •)  and  a  class -1C 
function  j(-)  such  that  for  all  k>R , 


\(X(t),x(t),£(t))\  <  max  j/?s(|(X(0),z(0),£(0),m(0))|,f),  7  |  Vf  >  0. 

(3.58) 

Then,  (3.22)  follows  by  setting  /?(•,  •)  =  /?«(-,  •)  and  k*  =  max{^_},  . ,  il}. 


Part  3:  m  G  C ^ 

To  prove  boundedness  of  m(t),  we  consider  the  differential  equations  (3.15)  and 
(3.19),  and  analyze  their  solutions  in  the  half-plane  H+  :=  {(y,m)\m  >  0}  which 
is  forward  invariant  from  (3.15).  It  is  sufficient  to  analyze  solutions  in  H+  because 
(3.15)  implies  that  solutions  satisfy  \m(t)\  <  |m(0)|e~5i  as  long  as  they  remain  outside 
H+ .  We  use  the  function 

U(y,  m)  =  +  9m,  (3.59) 

where  the  positive  constant  9  is  to  be  determined.  We  first  note  that  there  exist 
class-/Coo  functions  tti(-)  and  u2{-)  such  that,  V  ( y,m )  G  H+ , 

ui{\{y,m)\)  <  U(y,m )  <  u2(\(y,m)\).  (3.60) 
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Next,  from  (3.15)  and  (3.19),  U(y,m )  in  H+  satisfies 

U  <  -ky2  -  K\d\  \g(X,x)y\2m  +  6(-8m  +  \u\)  +  \g(X,  x)y\ (|c(£)|  +  | (\d\  -  1)|  \a{X,x)\). 

(3.61) 

Substituting 

|u|  <  | a(X,  x)|  +  k(1  +  | a(X,  x)\)\g{X,  x)y\  +  n\g(X,  x)y\  m  (3.62) 

in  (3.61),  and  defining 

R  :=  \g(X,x)y\(\c(£,)\  +  |(M  -  1)|  \a(X,x)\)  +6(\a(X,x)\  +  rc(l  +  \a(X ,  x)\)\g {X ,  x)y\) , 

(3.63) 


we  obtain 


6 


U  <-ky  -  K\d\(\g(X,x)y\  -  0—\g(X,x)y\  +  0—-  m  +  R.  (3.64) 


\d\ 


K\d\ 


We  pick  6  >  0  small  enough  to  guarantee  p2  —  @jj\P  +  >  0,  Vp  £  1R,  so  that  we 

can  find  q\  >  0  such  that 


\g{X,x)y\2  -6-^-\g{X,x)y\  + 9 -f-  )  >  qx. 


Substituting  in  (3.64),  we  get 


\d\ 


U  <  — (72P  +  R , 


n\d\ 


(3.65) 


(3.66) 


where  <72  =  min{2fc,  K^91}.  Since  R(t)  is  bounded,  (3.66)  implies  that  U(t)  is 
bounded.  Then,  m  £  Coo  follows  from  (3.60). 


Part  4:  GAS  and  LES 

Defining  A°L  =  r>Ad[^  [^=o  and  Fl  =  —  [x=(b  the  Jacobian  linearization 

of  the  (A,  y,  z)-subsystem  (3.16),  (3.17)  and  (3.55)  is 

z  =  +  B\X  +  b-py 

X  =  FLX  +  b3y  (3.67) 

y  =  -kci(1  +  \m\)y  +  c2y  +  b±z  +  b5X. 

Since  AQL  and  Fl  are  Hurwitz,  we  let  Pi  =  P{  >  0,  P2  =  P%  >  0  satisfy  A°LT P\  + 
PiA°l  =  —I,  FJP2  +  P2FL  =  —I,  and  use 

S  =  zTP1z  +  1XTP2X  +  y2,  7  >  0,  (3.68) 


43 


as  a  Lyapunov  function  for  the  nonlinear  ( X ,  y,  ^-subsystem  (3. 16), (3. 17), (3. 55).  The 
expression  for  S  shows  that  by  selecting  k  and  7  sufficiently  large,  we  can  find  con¬ 
stants  g  >  0  and  1?  >  0  such  that  S  <  —  gS  is  satisfied  for  all  ( X,y,z )  in  the  set 
11,9  :=  {(X,  y,,z)  :  S'  <  #}.  Thus,  if  a  trajectory  (X(t),y(t),z(t),m(t))  enters  x  1R , 
then  (X(t),y(t),  z(t))  — >  0  exponentially  as  t  — >  00.  This  also  ensures  rri(t)  — >  0  be¬ 
cause  of  (3.66),  where  i?  defined  by  (3.63)  as  a  function  of  ( X,y,z )  vanishes  at  zero. 
Choosing  e  >  0  sufficiently  small  in  (3.22),  we  guarantee  that  the  solutions  (X,  y,  z,  m) 
enter  x  1R  in  finite  time  and,  hence,  (X(t),y(t),  z(t),m(t))  —>  0.  Finally,  LES  is 
proved  with  the  help  of  the  Lyapunov  function  7 S'  +  U,  in  which  U  is  as  in  (3.66)  and 
the  constant  7  >  0  is  sufficiently  large.  □ 

3.4  Summary 

The  redesign  in  this  chapter  renders  a  class  of  nonlinear  control  laws  robust 
against  input  unmodeled  dynamics  which  are  relative  degree  zero,  globally  exponen¬ 
tially  stable,  and  have  ISS  zero  dynamics.  For  linear  unmodeled  dynamics,  the  latter 
two  conditions  are  equivalent  to  stability  and  minimum  phase  properties,  which  are 
less  restrictive  than  the  conditions  required  by  previous  redesigns.  Nonlinear  unmod¬ 
eled  dynamics  considered  here  also  differ  from  those  in  other  redesigns. 
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Chapter  4 

Nested  Saturation  Redesign 


The  redesigns  discussed  so  far  require  that  the  unmodeled  dynamics  be  relative- 
degree  zero  and  minimum-phase.  This  restricts  the  applicability  of  these  redesigns 
because  most  physical  examples  of  unmodeled  dynamics,  like  actuator  models,  have 
a  higher  relative  degree. 

For  systems  in  feedforward  form,  the  nested  saturation  design  of  Teel  [92,  93] 
employs  saturation  elements  to  guarantee  a  small-gain  property.  In  the  absence  of 
unmodeled  dynamics,  Teel’s  design  achieves  GAS  and  LES. 

In  this  chapter  we  redesign  nested  saturation  control  laws  to  guarantee  GAS 
and  LES  in  the  presence  of  input  unmodeled  dynamics.  The  redesign  is  applicable 
to  a  large  class  of  unmodeled  dynamics,  not  restricted  to  be  relative  degree  zero  or 
minimum  phase.  Using  the  asymptotic  small-gain  theorem  of  Teel  [93],  we  prove  not 
only  GAS,  but  also  an  asymptotic  gain  property  from  small  disturbances  to  the  states. 


4.1  Problem  Statement 

We  consider  the  locally  Lipschitz  system 

Xn  —  xn—  1  T  9n{x  1 1  X2 j  '  '  '  >  %n—  1  >  V ,  d ) 

X2  =  xi  +g2(xi,  v,d) 
xi  =  v  +  gi(v,d) 

£  =  q(t,u,d) 
v  =  p(£,u,d). 


(4.1) 


(4.2) 
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where  d  £  Mm  represents  disturbances,  £  £  Mp,  u,  v  £  1R,  the  functions  q(£,u,d)  and 
p(£,u,d)  vanish  at  zero,  and  gf  s  satisfy 

3i(u,0)  =  o(v) 

9i(xi,  ■  ■  ■  ,Xi-i,v,0)  =  o(xi,  ■  ■  ■  ,Xi-i,v),  i  =  2,  ■  ■  ■  ,  n,  (4.3) 


where  the  notation  g(v)  =  o(v)  means  that 


lim 

|v|— »0 


\9(v)  I 
M 


=  0. 


(4.4) 


For  the  nominal  system,  that  is  (4.1)  with  v  =  u  and  d  =  0,  the  nested  saturation 
design  guarantees  GAS  and  LES.  However,  the  following  example  illustrates  the  loss 
of  stability  due  to  unmodeled  dynamics: 


Example  4.1  For  the  system 


x 2  =  xi  +  g2(xi,v) 

xi  =  v  +  gi(v),  (4.5) 

the  nested  saturation  control  law  is 

v  =  u  =  -0i  (?/i  +  02(2/2)),  (4.6) 

where  2/1  =  x\,  2/2  =  x\  +  x2,  and  0«(-)’s  are  saturation  functions  with  saturation 
levels  selected  according  to  Teel  [92]. 

In  the  presence  of  the  input  unmodeled  dynamics 

v{s)  =  '  +  1,  u(s),  (4.7) 

sz  +  s  +  1 

the  resulting  closed-loop  system  is  unstable  as  verified  from  the  Jacobian  linearization. 

□ 

In  this  example  instability  is  established  from  the  linearization,  therefore  it  can¬ 
not  be  prevented  by  adjusting  the  saturation  levels.  To  guarantee  robustness  against 
unmodeled  dynamics,  we  introduce  the  controller  gains  k\,  ■  ■  ■  ,kn  in  the  redesigned 
control  law 


U  =  — 0!  {k\X\  +  02(/C2X2  H - h  <fn{knxn)  •••)),  (4.8) 

where  0j(-)’s  are  saturation  functions  as  defined  in  (1.56).  We  provide  a  procedure 
for  selection  of  the  gains  ki  and  the  saturation  levels  \  to  ensure  robustness. 
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For  unmodeled  dynamics  satisfying  Properties  4.1  and  4.2  below,  our  redesign 
achieves  an  asymptotic  gain  from  the  disturbance  d  to  the  state  (x,£)  when 

||d||a  :=  limsup \d(t)\ 

>  oo 

is  sufficiently  small.  In  particular,  if  ||d||0  =  0,  then  x(t)  and  f(t)  converge  to  the 
origin.  For  d(t)  =  0,  the  origin  (x,£)  =  (0,0)  is  GAS  and  LES. 

Property  4.1  The  vector  field  q(f,u,  0)  and  the  output  function  p(f,u,  0)  are  differ¬ 
entiable  at  the  origin  and  the  Jacobian  linearization 

=  f(9g(g,u,0)l  =  [ dg(£,  it,  0) 

:  L  ^  J  (0,0)  '  l  du 

c  ,=  f9p(^,tt,0)l  D  >=  f  dp(f,  u,  0) ' 

^  -  (0,0)  -  .  (0,0) 

is  such  that 

1.  The  system  matrix  A  is  Hurwitz, 

2.  The  dc  gain  5  =  D  —  CA~XB  >  0  is  positive. 

Property  4.2  For  £  =  q(f,u,d),  there  exist  constants  cu,Cd,Ad  >  0,  Au  >  0  such 
that,  if  u(t),  d(t)  are  locally  bounded,  ||u||0  <  Au,  ||d||0  <  Ad,  then,  for  each  initial 
condition  £(0),  the  solution  f(t)  exists  for  all  t  >  0,  and 

Hello  <  max{cu||u||a,crf||(i||a}.  (4.10) 

Inequality  (4.10)  expresses  an  asymptotic  gain  from  the  input  (u,d)  to  the  state  £. 
If  the  ^-subsystem  has  a  linear  input-to- state  stability  gain  as  in  Sontag  [86],  then 
Property  4.2  is  satisfied  with  Au  =  Ad  =  oo. 

4.2  Nested  Saturation  Redesign 

Substituting  v  =  p(f,u,  d)  and  using  (4.9),  we  rewrite  the  system  (4.1), (4. 2)  as 


Xn 

=  xn-i  +  Gn(xi,x2,  ■  ■  ■  ,xn-i ,f,u,d) 

Xl 

=  xi  +  G2{xi,f,u,d) 

(4.11) 

X\ 

=  Cf  +  Du  +  G\ (£,  u,  d) 

i 

=  q(€,u,d), 
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where  G i(£,  u,  d )  :=  p(£,  u,  d )  —  C £  —  Du  +  g\ (p(£,  u,  d),  d)  and,  for  i  =  2,  •  •  •  ,  n, 

Gi(x i,---  ,Xi-i,£,u,d)  :=  gi(x i,---  ,Xi-i,p(£,u,d),d).  (4.12) 

Because  of  (4.3),  GV s  satisfy 

Gi(^,u,0)  =  o(£,it)  /413v 

Gi(xi,  ■  ■  ■  ,  Xi~i,£,  u,  0)  =  o(xi,  ■  ■  ■  ,Xi-i,£,u),  i  =  2,  ■  ■  ■  ,n. 

The  saturation  levels  A*  >  0  and  the  gains  ki  >  0  in  (4.8)  are  designed  recur¬ 
sively:  In  the  first  step  we  set 


u  =  -  ui), 


(4.14) 


where  u\  is  to  be  designed  in  the  next  step.  Selecting  Ai  >  0  and  k±  >  0  according 
to  Lemma  4.1  below,  we  guarantee  that  the  composite  system 


6  := 
vi  \-- 


Xj 

£ 

X!  =  [1  0 


=  qi(Z.i,ui,d)  :  = 

0 ]Ci  =:  C& 


C£  +  Du  +  G\  (£,  u,  d) 

q(£,u,d) 


u=—(f)i(kixi—ui) 


(4.15) 


with  input  (u\,d)  and  output  v\  possesses  the  same  properties  as  the  unmodeled 
dynamics  subsystem,  that  is,  Property  4.2  is  satisfied  with  q,  £,  u  replaced  by  qi,  £i,  u\, 
and  Property  4.1  is  satisfied  with  the  Jacobian  linearization  A i,  B±  of  the  vector  field 
gi(£i,  ui,  0),  D i  =  0  and  C\  as  in  (4.15). 

In  the  second  step  we  treat  the  =  [xi,£t]t  subsystem  of  (4.11)  with  input 
(ui ,  d)  and  output  v\  =  C\^\  as  the  “virtual”  unmodeled  dynamics  subsystem  for  the 
(x'2,  •  •  •  ,xn)  subsystem  and  rewrite  (4.11)  as 


Xn  —  Xn—\  +  Gn  2  ■  ,  Xn—i,  ,  U\ ,  d) 


X‘2  =  Cit,i  +  (*2,2  (£i,  ui)  d)  (4.16) 

£l  =  Ql(£l:Ul,d), 

where,  for  j  =  2,  ,n,  the  function  Gjt  2  is  obtained  by  substituting  (4.14)  in  Gj. 

We  note  that 

G2, 2(^,1, u\,0)  =  o(£i,  ui)  (4  17) 

Gj,2(x2,-”  =  o(x2,---  ,Xj-l,£l,Ul),  j  =  3,---,ra, 
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and  proceed  to  design 


til  =  —4>2(k2X2  -  U2) 


(4.18) 


such  that  the  composite  system 


£2  ;— 


®2 

£1 

V2  ■=  X2  =  [1  0  0 


=  q2(&,u2,d)  := 


C&  +  G2,2(£i,ui,d) 

qi(£i,ui,d) 

0]^2  =:  C2& 


Ul=— (f)2{k‘2X2  —  U2) 


(4.19) 


with  input  (u2,d)  and  output  V2  satisfies  Properties  1  and  2. 

This  recursive  design  relies  on  the  assumption  that  at  the  ith  step  A*  and  ki 
can  be  selected  such  that  the  ^-subsystem  with  input  (v,i,d)  and  output  v,L  satisfies 
Properties  1  and  2.  We  now  prove  that  this  is  guaranteed  by  selecting  both  Aj  >  0 
and  ki  >  0  sufficiently  small. 


Lemma  4.1  Consider  the  system  (4-11)  and  let  fi  he  defined  by  £0  :=  £>  Ci  := 
[xi,£j_i\T ,  i  =  1,  •  •  •  ,n.  Let  the  control  u  be  constructed  recursively  by 

Ui-i  =  -fifikiXi  -  m),  i  =  l,  (4.20) 

where  is  a  saturation  function  with  saturation  level  A,;,  u$  =  u  and  un  :=  0. 
Let  q0  '■=  q,  Gpi  :=  Gi,  C0  :=  C,  D0  :=  D,  and  Cl  :=  [1  0ix(p+i_i)],  A  :=  0, 
i  =  1,  •  •  •  ,n,  and  write  the  £j- subsystem  of  (4-11)  as 

±i  =  Vi-i  +  GM(^_i,Uj_i,d) 

£i—  1  =  qi—l{fii—l-iui—hd\  (4.21) 

where  G^(&-i> «i-i,  0)  =  o(^_i,Uj_i). 

If  the  fi- 1  -sub system  with  input  (uj_i,d)  and  output 

Vi— 1  =  Ci-i^i-i  +  Di—iUi—i  (4.22) 


possesses  Properties  1  and  2,  then  there  exist  A*  >  0  and  k*  >  0  such  that,  with 
A i  6  (0,  AT]  and  ki  6  (0,  A;*]  in  (4-%0),  the  subsystem 


ti 


Xi 

ti- 1 


qi(€i,ui,d) 

Ci—l^i—l  +  Di-iUi-i  T  Gi  i (£j_  1 ,  Ui—\ ,  d) 

qi-i(fi-i,Ui-i,d) 


Uj, —  1 —  Ui) 

(4.23) 
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with  input  ( Ui ,  d)  and  output 


Vi  :=  Xi  =  C&  +  DiUi  (4.24) 

satisfies  Properties  1  and  2.  □ 

The  proof  is  given  in  Section  4.3.  Our  control  law  is  constructed  by  recursively 
applying  Lemma  4.1  for  i  =  1,  •  •  •  ,  n  to  select  A j  and  k{.  At  the  nth  step,  the 
saturation  level  An  and  the  gain  kn  in 


'U'n—l  • —  4>n(knXn') 


(4.25) 


guarantee  that  the  closed-loop  system 


£n—  1 


qn(in,d) 


%n—  1  Gn^n  (^n—  1 5  1 5  ^0 

Qn—  1  (£n—  1 5  1 5  ^0 


1 —  tpniknXn) 

(4.26) 


possesses  Properties  1  and  2,  which  yield  the  following  result: 


Theorem  4.1  Consider  the  system  (4-1)  with  the  unmodeled  dynamics  (4-2)  satisfy¬ 
ing  Properties  1  and  2.  If  the  control  law  (4-8)  is  constructed  by  recursively  applying 
Lemma  4-1  for  i  =  1,  •  •  •  ,  n,  then  there  exist  A  >  0  and  c  >  0  such  that, 

1.  If  d(t)  is  locally  bounded  and  ||d||a  <  A,  then  (x(t),f(t))  exist  for  all  t  >  0 

and 


ll(z>Olla  <  4d\\a-  (4.27) 

2.  If  d(t)  =  0,  then  the  origin  is  globally  asymptotically  stable  and  locally  expo¬ 
nentially  stable. 

Proof:  Inequality  (4.27)  follows  directly  from  Property  4.2.  If  ||d||a  =  0,  then  (4.27) 
guarantees  global  attractivity  of  the  origin  (x,f)  =  (0,  0).  If,  in  addition,  d(t)  =  0,  we 
have  LES  because  the  Jacobian  linearization  An  of  qn(fm  0)  is  Hurwitz  by  Property 
4.1.  Since  the  resulting  closed-loop  system  is  time-invariant,  global  attractivity  and 
LES  together  imply  GAS.  □ 

It  is  important  to  note  that  this  redesign  does  not  require  detailed  knowledge  of 
the  unmodeled  dynamics  (4.2).  Indeed,  from  the  proof  of  Lemma  4.1,  the  selection 
of  A i  and  ki  is  based  on  upper  bounds  on  5  =  D  —  CA~1B,  |CA4_1|,  cu  and  Au  in 
Properties  1  and  2,  and,  on  functions  that  form  upper  bounds  on  |p(£,  u,  0)|,  |g(£,  u,  0)| 
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and  \p{f,  u,  0)|,  where  q  and  p  represent  the  deviations  of  q(f,  u,  0)  and  p(f,  u,  0)  from 
their  linearizations  +  Bu  and  +  Du,  respectively.  Therefore,  the  control  law 
constructed  as  in  Theorem  4.1  guarantees  robustness  against  a  class  of  unmodeled 
dynamics  characterized  by  these  upper  bounds. 

Example  4.2  We  now  apply  the  redesign  to  the  system  of  Example  4-  1,  with  <72(24,  v)  = 
x\  +  v2  and  g\(v)  =  v2.  The  linear  unmodeled  dynamics  subsystem  (4-7)  is  Hurwitz 
and  has  a  positive  dc  gain  5  =  1,  thus  it  satisfies  Properties  1  and  2.  The  redesigned 
nested  saturation  control  law 


u  =  ~(pi(kixi  +  (fo(k2x2 )),  (4.28) 

with  parameters  Ai  =  0.5,  k\  =  0.4,  A2  =  0.05,  k2  =  0.02  selected  according  to 
Lemma  f.l  renders  the  closed-loop  system  (4-5),  (4-7),  (4-28)  GAS,  as  illustrated  by 
numerical  simulation  in  Figure  4-1.  The  plots  of  x2  and  u,  along  with  the  small  values 
of  X2  and  k2  illustrate  the  low-gain  character  of  this  design. 
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Figure  4.1:  Closed-loop  solutions  with  the  redesigned  control  law  (4.28). 
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4.3  Proof  of  Lemma  4.1 


We  give  the  proof  for  %  =  1,  that  is,  we  show  that  if  Ai  >  0  and  k\  >  0  in  (4.14) 
are  selected  sufficiently  small,  then  the  £ i -subsystem  (4.15)  with  input  u\  and  output 
v\  satisfies  Properties  1  and  2. 

We  prove  Property  4.2  by  using  the  asymptotic  small-gain  theorem  of  Teel  [93]. 
To  put  the  system  (4.15)  in  a  form  suitable  for  the  small-gain  formulation  we  write 


<?(£,  u,  d)  =  A£  +  Bu  +  q(£,  u,  d),  (4.29) 

where  q(£,u,0)  =  o(£,u).  From  (4.15)  and  (4.29),  the  variable 

z  :=  x i  —  Ch4_1£  (4.30) 

is  governed  by 

z  =  (D  -  CA~lB)u  +  G{i,u,d),  (4.31) 


where  G(£,u,d)  :=  Gi(£,u,d)  —  CA  lq(^,u,d),  thus  G(£,u,  0)  =  o(£,u).  We  substi¬ 
tute  5  =  D  —  CA~lB  >  0  in  (4.31),  and  write  (4.15)  as 


From  (4.14),  we  have 


z  =  5u  +  G{^,u,d) 

£  =  q(£,u,d). 

5u  =  -5(f>i  (ki(xi  -  J-) 


(4.32) 


(4.33) 


It  follows  from  the  definition  of  saturation  functions  (1.56)  that  4>(s)  :=  6(f>i(s/5)  is 
also  a  saturation  function  with  level  A  :=  4Ai .  Defining  k  :=  Ski,  we  write  (4.33)  as 


Su 


-*  {Hzi  -  £>)  - 


Substituting  (4.34)  in  (4.32),  and  using  (4.30),  we  obtain 


Defining 


i 

£ 


-4>  [k(z  +  CA~li 
q(Z,u,d). 


+  G(£,u,  d) 


Vi 


(2)  CA-i 


e,  vr-= 


(2)  G(£,u,d) 


k 


(4.34) 


(4.35) 


52 


^:=d,  2  +  j,<2>  + 

K\ 

we  represent  (4.35)  as  in  Figure  4.2,  where  the  subsystems  £(1*  and  E(2^  are 

:  i  =  —0  ( k{z  +  y +  i/^)^  +  icy ^  (4.36) 

S(2):  £  =  g(C,u,i/(2)),  u  =  -<&i(fcU/(1))-  (4-37) 


Figure  4.2:  The  feedback  interconnection  (4.36)-(4.37). 


We  now  prove  that  if  \\  and  Aq  are  sufficiently  small,  then  the  -subsystem 
(4.15)  with  input  (u\,d),  represented  as  in  Figure  4.2,  satisfies  Property  4.2.  We 
will  equivalently  show  that  if  A,  k  are  sufficiently  small,  then  the  composite  output 
(y(’),  y(2))  of  the  feedback  interconnection  in  Figure  4.2  with  input  satisfies 

the  asymptotic  gain  property 

||(y(1),J/(2))IU  <  maxfc^i  ||r/(1)||a,  ^2  ||i/(2) ||a } ,  (4.38) 

Ci/,i,  c,/^  >  0,  for  small  ||i/^||0  and  To  prove  (4.38),  we  rely  on  the  following 

result  adapted  from  Teel  [93,  Theorem  2]. 

Proposition  4.1  For  the  feedback  interconnection  in  Figure  f.2,  suppose 

Al.  For  all  locally  bounded  i/W,  i/(2),  the  solutions  are  unique  and  there  are  no 
finite  escape  times, 

(2) 

A 2.  There  exist  constants  Ay  >  0,  ci,  C2,  C3  >  0  such  that,  if  \\y2  q|a  <  Ay,  then 
||2/(1)IU  <  max{c!  ||yS2)||a,  c2  ||y£2)||a,  c3  ||^(1)||a},  (4.39) 

A3.  There  exist  continuous,  nondecreasing  functions  7i(-),72()  satisfying 
7i(0)  =  72(0)  —  0,  71(00)  <  00,  72(00)  <  00  and  72(00)  <  Ay,  (4.40) 
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and  there  exist  constants  Au,  04,05  >  0  such  that,  if  ||z^2)||a  <  Au,  then 

Wy^Wa  <  max{7i(||y(1)||a),  c4  ||^(2)||a}  (4.41) 

h^Wa  <  max{72(||y(1)||a),  c5  ||;/(2)||a},  (4.42) 

Af.  For  all  s  /  0, 

ci  71  (s)  <  sj.  c272(s)<s.  (4.43) 

Then,  there  exist  constants  cV)\,  cv^,  A,y.2  >  0,  A^i  >  0  such  that  (4-38)  holds  for 
||^(1)||a  <  and  ||^(2)||a  <  AVj2.  □ 


We  now  prove  that  if  Ai  and  k\  are  sufficiently  small,  then  the  conditions  A1 
through  A4  of  Proposition  4.1  are  satisfied  for  X/1)  and  X/2)  defined  in  (4.36)  and 
(4.37). 

Uniqueness  of  solutions  follows  from  the  Lipschitz  continuity  of  the  right  hand 
sides  of  (4.36)  and  (4.37).  The  absence  of  finite  escape  times  can  be  argued  from  the 
boundedness  of  u,  as  in  Teel  [93,  Lemma  3.5]. 

To  prove  A2,  we  introduce  the  Lyapunov  function  V(z )  =  z2  for  (4.36),  and 
note  that  <  A/2 k  and  \z\  >  +  v^\  +  [y^l  together  imply  V  <  0.  Thus,  if 

||y22^||a  <  A/2 k  =  Ai/2Aq,  then 

\\z\\a  <  II y[2)  +  ^(1)  |a  +  ||y22)  lla-  (4.44) 

It  follows  from  (4.44)  and  y^>  :=  z  +  y ^  that  if  ||y22,||a  <  Ai/2fci,  then 

||y(1)||a  <  2||y(J)  ||a  +  ||y22)||a  +  2||i'(1)||a  <  max{6||yf)  ||0, 6|| y[22 )  1 1 a,  6|| u(1)  ||0} .  (4.45) 


We  conclude  from  (4.45)  that  A2  holds  with  <7  =  c2  =  C3  =  6  and  Ay  <  Xi/2k±,  say, 

Ai 


Av  =  ,  . 
y  3k\ 


(4.46) 


To  prove  A3,  we  recall  from  Property  4.2  that,  if  ||u||a  <  Au,  that  is,  Ai  <  At 
and,  if  ||z/®||a  <  A^,  then 


IlClla  <  max{Cu  1 1  U 1 1 a  ,  Cd  ||^(2)  ||a} 


(4.47) 


From  (1.56),  we  have  |«|  =  \4>i(kiy^)\  <  min  {fci|y(1)|,Ai}.  Substituting  in  (4.47), 
:=  CA _1£  satisfies  (4.41)  with 

7i(s)  =  |Cb4_1|  cumm{kis,  Ai},  c4  =  | CA~l\cd,  A„  <  Ad.  (4.48) 
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It  follows  from  (4.48)  that  71(00)  =  \CA  x|  cuX\  <  00  and,  for  all  s  >  0, 

7i(s)  <  facu\CA~l\s.  (4.49) 

To  show  that  (4.42)  holds  for  tpp  :=  we  use  Lipschitz  continuity  of  G(£,u,d) 

and  note  that,  if  £,  u,  d  belong  to  a  compact  set,  then  3L  >  0  such  that 

\G(fau,d)\<\G(fau,0)\+L\d\.  (4.50) 

By  selecting  Av  to  be  finite,  we  restrict  ||d||0  <  Au  to  a  compact  interval.  ||it||a  and 
||^||0  also  belong  to  a  compact  set  because  |u|  <  Ai.  Then,  from  G(£,u,  0)  =  o(£,u), 
and  the  gain  property  (4.10),  there  exists  Lq  >  0  and  a  nondecreasing,  continuous 
function  70  (s)  =  o(s)  such  that 

\\G(tu,d)\\a  <  max{7o(||u||a),L0  ||d||a},  (4.51) 

thus  (4.42)  is  satisfied  with 

72 («)  =  ^7o(min{A;1s,  Ai}) ,  c5  =  (4.52) 

It  is  clear  that  72(00)  <  00  as  in  (4.40).  Comparing  72(00)  =  ^770(^1)  to  Ay  in 
(4.46),  and  recalling  7o(s)  =  o(s),  we  conclude  that  72(00)  <  Ay  is  satisfied  if  Ai  is 
selected  sufficiently  small.  From  (4.52),  a  linear  bound  on  72  (s)  is 

72  (s)  <  res,  (4.53) 

where  k  :=  sup^^  ^  =  sup^^  Defining  s  :=  fas,  we  get 

k  =  sup  ^0}-,  (4.54) 

«e[o,Ai]  ° s 

which  is  independent  of  fa.  Moreover,  since  70 (5)  =  o(s),  k  can  be  rendered  arbitrarily 
small  by  making  Ai  sufficiently  small. 

We  prove  A4  using  the  upper  bounds  on  71  (s)  and  72(5)  given  in  (4.49)  and 
(4.53),  respectively.  From  (4.49),  the  first  inequality  in  (4.43)  is  guaranteed  by  select¬ 
ing  fa  such  that 

ci£:icu|Ck4-1|  <  1.  (4.55) 

From  (4.53),  the  second  inequality  in  (4.43)  is  guaranteed  by  selecting  Ai  and,  hence, 
k  sufficiently  small. 
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To  show  that  Property  4.1  holds  for  (4.15),  we  note  that  the  Jacobian  lineariza¬ 
tion  of  <7i(£i,«i,0)  is 


Ai  = 


-fci  D  C 
-kiB  A 


B  i 


D 

B 


Using  C\  =  [1  Oixp],  =  0,  and 


ki 


—A~1B  A-1 


-j^CAr1 

+  A~1BCA~ 1  ’ 


(4.56) 


(4.57) 


we  obtain 

Si  =  D1  -  CiA^Bi  =  - -(D-  CA~1B)  >  0  .  (4.58) 

With  k\  selected  to  satisfy  the  small-gain  condition  (4.55),  Ai  is  also  guaranteed 
to  be  Hurwitz.  This  follows  because  the  small  gain  analysis  for  (4.35)  also  holds  with 
G  =  0,  ui  =  0,  d  =  0,  and  with  q(f,u )  :=  q{f,u,  0)  replaced  by  its  linearization 
H£  +  Bu , 


i  =  -k{z  +  CA^f)  (4.59) 

£  =  H£  +  Bu,  u  =  —  k\{z  +  Cb4_1£), 


which  is  £i  =  A^i  expressed  in  the  (2,  £)  coordinates.  In  Lemma  4.2  below,  we  show 
that  the  linearization 


£  =  H£  +  Bu,  (4.60) 

has  the  same  asymptotic  input-to-state  gain  as  the  nonlinear  system  £  =  q(f,u). 
Then,  since  the  small-gain  condition  (4.55)  holds,  the  asymptotic  small-gain  theorem 
guarantees  attractivity  of  the  origin  for  £i  =  Hi£i,  which  is  equivalent  to  Hurwitz 
stability  of  Hi.  □ 

Lemma  4.2  Let  £  =  q(£,u)  be  such  that  its  Jacobian  linearization  (4-60)  is  Hurwitz, 
and  for  each  locally  bounded  u  satisfying  ||u||a  <  Au,  Au  >  0,  and  for  each  initial 
condition  £(0),  the  solution  £(t)  exists  for  all  t  >  0  and  satisfies 

Il£||a  <  CulMlo-  (4-61) 

Then,  for  each  locally  bounded  u,  the  solutions  of  (4-60)  also  satisfy  (4-61). 
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Proof:  Suppose,  on  the  contrary,  there  exists  a  locally  bounded  u*  such  that  the 
solution  of  the  linearization  (4.60)  satisfies 

||£l|U  =  c||u*||a,  c  >  cu.  (4.62) 

Since  (4.60)  is  a  linear  system,  multiplying  u *  by  a  constant  does  not  change  (4.62), 
thus  we  can  set  ||u*||0  =  e  <  Au.  We  rewrite  £  =  </(£,  u)  as 

£  =  4£  +  Bu  +  q(£,u),  (4.63) 


where  q(£,u)  =  o(£,  u),  and  note  that  its  solution  with  u  =  u*  can  be  expresses  as 
£(i)  =  £i(t)  +  £g(i),  where  £g  is  due  to  g(£,tt*),  and, 


||£||a  >  |c||lt*||o  -  ||£q||a|  = 


C  — 


lie 


q\\a 


U 


(4.64) 


Because  of  (4.61),  ||g(£, u*)||a  =  o(e),  and,  hence,  ||£g||a  =  o(e).  Thus,  for  sufficiently 
small  e,  \c  —  |  >  cu,  which  contradicts  (4.61).  □ 


4.4  Summary 

We  have  presented  a  redesign  of  nested  saturation  control  laws  that  makes  them 
robust  against  input  unmodeled  dynamics,  without  any  restrictions  on  their  zero 
dynamics  or  relative  degree.  The  achieved  robustness  property  is  due  to  the  low-gain 
design,  which  is  sufficient  for  stabilization  of  systems  in  feedforward  form.  Without 
the  feedforward  structure,  global  stabilization  may  not  be  possible  in  the  presence  of 
unmodeled  dynamics,  because  low-gain  designs  may  not  be  applicable. 
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Part  II 

Output-Feedback  Designs 
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Chapter  5 

A  New  Nonlinear  Observer 


The  control  designs  we  have  presented  in  Part  I  assume  that  the  full  state  of  the 
plant  is  measured.  We  now  remove  this  assumption  and  consider  the  problem  where 
only  a  part  of  the  state  is  available  for  measurement.  In  the  rest  of  the  dissertation 
we  introduce  new  tools  for  nonlinear  observer  and  observer-based  control  designs,  and 
use  them  to  design  robust  output-feedback  control  laws  for  systems  with  unmodeled 
dynamics. 

In  this  chapter  we  introduce  a  new  global  observer  design.  The  class  of  systems 
for  which  our  observer  is  applicable  are  characterized  by  two  restrictions  which  allow 
the  observer  error  system  to  satisfy  the  multivariable  circle  criterion.  First,  a  linear 
matrix  inequality  (LMI)  is  to  be  feasible,  which  implies  a  positive  real  property  for 
the  linear  part  of  the  observer  error  system.  The  second  restriction  is  that  the  non- 
linearities  be  nondecreasing  functions  of  linear  combinations  of  unmeasured  states. 
This  restriction  ensures  that  the  vector  time-varying  nonlinearity  in  the  observer  er¬ 
ror  system  satisfies  the  sector  condition  of  the  circle  criterion.  The  observer  design  in 
Section  5.1  is  accompanied  by  its  reduced-order  variant  in  Section  5.2. 

The  proposed  observer  design  is  constructive  in  the  sense  that  the  issues  of  exis¬ 
tence  and  the  evaluation  of  observer  matrices  satisfying  the  circle  criterion  are  resolved 
by  efficient  LMI  computations.  A  further  advantage  of  our  design  is  its  robustness 
against  inexact  modeling  of  nonlinearities.  This  robustness  property  is  analyzed  in 
Section  5.3,  and  bounds  are  given  within  which  the  observer  error  gradually  increases 
with  an  increase  in  the  modeling  error. 
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5.1  The  New  Observer  Design 


We  now  present  the  observer  design  for  the  control  system 

x  =  Ax  +  G^(Hx)  +  g(y,u)  (5.1) 

y  =  Cx, 

where  x  £  Mn  is  the  state,  y  6  Mp  is  the  measured  output,  u  £  Mm  is  the  control 
input,  the  pair  ( A,C )  is  detectable,  and,  q(-)  and  g(-,-)  are  locally  Lipschitz.  The 
state-dependent  nonlinearity  7 (Hx)  is  an  r-dimensional  vector  where  each  entry  is  a 
function  of  a  linear  combination  of  the  states 

n 

7 i  =  liC^2Hijxj),  *  =  1,  ■■■  ,r.  (5.2) 

3= 1 

The  main  restriction  is  that  each  7 j(-)  is  nondecreasing,  that  is,  for  all  a,  b  £  M,  it 
satisfies 

(a  -  &)[7i(a)  -  7i(6)]  >  0.  (5.3) 

If  7 j(-)  is  continuously  differentiable,  then  (5.3)  means  that  d^i(v)/dv  £  [0,  00)  for 
all  v  £  IR.  If,  instead,  7j(-)  satisfies  d-)i(y)/dv  £  [^,00)  with  ^  /  0,  we  can  still 
represent  the  system  as  in  (5.1)-(5.3)  by  defining  a  new  function  7 j(u)  :=  7i(n)  —  g{V 
which  satishes  dr/i (v)/ dv  £  [0, 00),  and  absorbing  giV  in  the  linear  part  of  the  system. 
The  observer  for  system  (5.1)-(5.3)  is 

x  =  Ax  +  L(Cx  —  y)  +  G^(Hx  +  K(Cx  —  y))  +  g(y,  u ) ,  (5.4) 

and  the  design  task  is  to  determine  the  observer  matrices  K  £  Mrxp  and  L  £  Mnxp. 
Note  that  the  nonlinear  injection  7 (Hx  +  K(Cx  —  y ))  is  analogous  to  tc3  in  the  van 
der  Pol  example  in  Section  1.2.4.  The  same  observer  design  can  be  applied  to  the 
system  (5.1)  when  the  nonlinearity  7 (Hx)  also  depends  on  y  and  u.  In  this  case  we 
require  that  the  nondecreasing  property  (5.3)  hold  for  each  y  £  1RP  and  u  £  Mm. 

For  the  observer  equation  to  be  defined,  the  uniqueness  of  the  solutions  x(t)  of 
(5.1)  is  guaranteed  by  restricting  the  control  law  u  =  a(y,x,t)  to  be  locally  Lipschitz 
in  (y,x),  uniformly  in  t.  As  will  be  further  clarified  in  Chapter  7,  u  =  a(y,x,t )  is 
also  assumed  to  ensure  that  x(t)  does  not  escape  to  infinity  in  finite  time,  that  is 
x(t)  £  £%. 

From  (5.1)  and  (5.4),  the  dynamics  of  the  observer  error  e  =  x  —  x  are  governed 
by 

e  =  (A  + LC)e  +  G[y(v) (5.5) 
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where 


v  :=  Hx  and  w  :=  Hx  +  K(Cx  —  y) .  (5.6) 

We  begin  the  observer  design  by  representing  the  observer  error  system  (5.5)  as  the 
feedback  interconnection  of  a  linear  system  and  a  multivariable  sector  nonlinearity. 
To  this  end,  we  introduce  a  new  variable 

z:=v-w  =  {H  +  KC)e,  (5.7) 

and  view  y(n)  —  7 (w)  as  a  function  of  t  and  z 

p(t,z)  :=  j(v)  -  y(w) ,  (5.8) 

where  the  time  dependence  of  <p(t,z )  is  due  to  v(t)  =  Hx(t).  Substituting  (5.8),  we 
rewrite  the  observer  error  system  (5.5)  as 

e  =  {A  +  LC)e  +  Gcp(t,  z)  (5-9) 

z  =  {H  +  KC)e, 

and  note  from  (5.3)  that  each  component  of  <p(t,  z)  satisfies 

Zi<pi(t,  Zi)  >  0  ,  \/zi£lR.  (5.10) 

Thanks  to  this  sector  property,  we  can  employ  the  multivariable  circle  criterion  and 
derive  a  condition  that  guarantees  the  exponential  convergence  of  the  observer  error 
e(t)  to  zero. 

Theorem  5.1  Consider  the  plant  (5.1),  observer  (5-4),  and  suppose  x(t)  exists  for 
all  t  >  0.  If  there  exists  a  matrix  P  =  PT  >  0 ,  a  constant  v  > 
matrix  A  >  0  such  that 

'  (A  +  LC)tP  +  P(A  +  LC)  +  ul  PG  +  {H  +  I<C)TA  ' 

GtP  +  A  {H  +  KG)  0 

then  the  observer  error  eft)  satisfies  for  all  t  >  0 

|e(i)|  <  ac| e(0) |  exp(— fit) ,  (5.12) 

where  K  =  \/fcS h 13  =  2 


0,  and  a  diagonal 

<  0,  (5.11) 
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Proof:  From  (5.9),  the  time  derivative  of  V(e)  =  e1  Pe  is 

V  =  eT[(A  +  LC)tP  +  P(A  +  LC)]e  +  2 eTPGp(t,  z) ,  (5.13) 

and,  in  view  of  (5.11), 

V  <  -v\e\2  -  2 eT(H  +  KC)T  A(p(t ,  z) .  (5.14) 

Substituting  ( PI  +  KC)e  =  z,  we  rewrite  (5.14)  as 

r 

V  =  -v\e\2  -  2  ^  \iZi<pi(t,  Zi) ,  (5.15) 

1=1 

where  A;  >  0,  i  =  1,  •  •  •  ,  r,  are  the  entries  of  the  diagonal  matrix  A.  Then,  the  sector 
property  (5.10)  yields 

V  <  —  u\e\2  ,  (5.16) 

from  which  (5.12)  follows.  □ 

With  Theorem  5.1,  the  observer  design  for  system  (5.1)  is  reduced  to  the  problem 
of  finding  observer  matrices  K  and  L  such  that  (5.11)  is  satisfied  with  some  P  =  P 1  > 
0,  A  >  0,  and  v  >  0.  The  existence  of  such  K  and  L  depends  on  A,  C,  G  and  H ,  and 
cannot  be  ascertained  a  priori.  However,  (5.11)  is  a  LMI  in  P,  PL,  A,  A K  and  u. 
Therefore,  we  can  use  the  efficient  numerical  tools  available  for  LMI’s  to  determine 
whether  the  problem  is  feasible  and,  if  so,  to  compute  K  and  L. 

We  have  gained  additional  design  freedom  in  the  LMI  (5.11)  by  introducing 
A  =  diag(X\,  ■  ■  ■  ,  Xr)  as  a  parameter  instead  of  A  =  I.  Due  to  the  special  struc¬ 
ture  of  <p(t,z ),  in  which  every  individual  entry  <pi(t,Zi)  satisfies  the  sector  condition 
Zi<fi(t,Zi)  >  0,  the  sum  Yli=i  XiZppiit,  zl)  =  ip(t,z)TAz  is  nonnegative  for  any  choice 
of  Aj  >  0,  i  =  1,  •  •  •  ,r.  This  means  that  the  framed  block  in  the  feedback  loop 
in  Figure  5.1  below  is  a  multivariable  sector  nonlinearity  for  any  diagonal  A  >  0. 
Thus,  from  the  multivariable  circle  criterion,  asymptotic  stability  of  the  closed-loop 
is  guaranteed  if  the  linear  system  with  input  i?  and  output  A z  is  SPR.  Indeed,  the 
LMI  (5.11)  constitutes  the  required  SPR  condition. 

If  the  LMI  (5.11)  holds  with  v  =  0,  that  is,  if  the  linear  system  in  Figure  5.1  with 
input  iJ  and  output  A z  is  PR  but  not  SPR,  then  (5.16)  guarantees  |e(f)|  <  «;|e(0)|, 
but  not  necessarily  e(t)  — ►  0.  Under  the  additional  assumption  that  x(t )  £  Coo,  we 
now  prove  that  e(t)  still  converges  to  zero  exponentially. 
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Figure  5.1:  Observer  error  system. 


Theorem  5.2  Consider  the  plant  (5.1)  and  observer  (5.4).  If  P  =  PT  >  0  satisfies 
(5.11 )  with  v  =  0,  and  A  +  LC  is  Hurwitz,  then  for  any  compact  set  C  C  M2n,  with 
( x(t),x(t ))  6  C  for  all  t  >  0,  there  exist  positive  constants  nc  and  fie  such  that 

|e(f)|  <  «c|e(0)|exp(— /?cf),  Vf>0.  (5.17) 


Proof:  We  will  use  the  following  lemma  proved  by  Teel  [94]: 

Lemma  5.1  If  there  exist  positive  constants  9\  and  62  such  that,  for  each  to  >  0, 
and  for  all  t  >  to, 


<  0i|e(to)| 

<  O2\e(t0)\2, 


then  e(t )  satisfies  (5.17)  with  constants  k c,fic  >0  that  depend  on  9\,  02. 


(5.18) 

(5.19) 


The  first  inequality  (5.18)  holds  because  V  <  0  by  (5.16).  To  prove  (5.19),  we  first 
note  that  the  observer  error  system  (5.9)  driven  by  ip(t,  z )  is  testable,  because  A+LC 
is  Hurwitz.  This  means  that  there  exist  >  0  and  c  >  0  such  that 


[  \e(r)\2dT  <bi\e(t0)\2  +c  f  |y?(r, z(r))|2dr . 

Jto  Jt0 

To  prove  that  (5.19)  holds,  we  show  that 


/ 

Jto 


|^(r,2(r))|2dr  <  &2|e(t0)|2. 


(5.20) 


(5.21) 
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Since  (x(t),x(t))  6  C  for  all  t>  0,  the  components  Vi,Wi  of  v  and  w,  defined  in  (5.6), 
remain  for  all  t  >  0  in  a  compact  interval  Ct,  in  which 


di 


sup 

Vi,Wi&Ci,Vi^Wi 


7 i(vj)  - 

Vi  -  Wi 


(5.22) 


is  finite,  because  7 j(-)  is  locally  Lipschitz.  It  follows  from  (5.3),  (5.22),  and  Zi  =  Vi~Wi 
that  whenever  z*  0,  the  inequality 

0  <  <  ck  (5.23) 

holds  for  all  t  >  0.  Multiplying  the  right  inequality  in  (5.23)  by  ^Zi(p(t,  Zi),  which  is 
nonnegative  by  the  left  inequality,  we  obtain 

-j-¥>i(i,  zi )2  ^  zm(t,  Zi) ,  (5.24) 

di 

and,  in  view  of  (5.15), 


A 

iznpi(t,zi)  <  -^^i(t,Zif  <  —b  \ip(t,  z)\2  , 

i  ''' 


(5.25) 


where  b  :=  min, { ^ }  >  0.  Integrating  both  sides  of  (5.25)  from  to  to  t,  we  verify  that 
(5.21)  holds  with  b-2  =  Xma^p^ .  Substituting  (5.21)  in  (5.20),  we  see  that  (5.19)  holds 
with  6-2  =  b\  +  062  •  hi 


We  now  illustrate  the  use  of  Theorem  5.2  in  a  situation  where  the  LMI  (5.11) 
is  not  feasible  for  v  >  0. 


Example  5.1  The  system 

±1  =  X2  +  x\ 

x 2  =  X2  +  X3  -  exp(x2)  +  u 

x3  =  2  u 

V  =  x  i 

is  of  the  form  (5.1)  with  g(y,u)  =  [y2  u  2 u]T , 


"  0 

1 

0  " 

0 

A  = 

0 

1 

1 

,  C=  [1  0  0],  G  = 

-1 

0 

0 

0 

0 

(5.26) 


H  =  [0  1  0]  . 
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Since  the  nonlinearity  7 (Hx)  =  71(2:2)  =  exp(,T2)  is  nondecreasing  as  in  (5.3),  we 
proceed  with  the  observer  design.  With  A,  C,  G  and  H  as  above,  the  LMI  (5.11)  is 
not  feasible  with  v  >  0.  However,  with  v  =  0,  it  is  feasible  and  a  solution  is 

3  -1  -1 1  r-2- 

P=  -1  1  0  ,  L  =  -4  ,  K  =  - 1,  A  =  1 .  (5.27) 

-1  0  1  J  |_  -1 

Because  A  +  LC  is  Hurwitz,  we  conclude  from  Theorem  5.2  that  the  resulting  observer 
x\  =  x2-  2{xi  -  y)  +  y2 

x2  =  x2  +  £3  ~  4(27  -  y)  -  exp(x2  -  (£i  -  y))  +  u  (5.28) 

£3  =  -{x1-y)  +  2u 

guarantees  e(t)  —*■  0  as  t  — >■  00,  if  a  control  law  can  be  designed  to  ensure  x(t )  G  £oo. 

5.2  Reduced-Order  Observer 

In  applications  it  may  be  more  convenient  to  employ  a  reduced-order  observer, 
which  generates  estimates  only  for  the  unmeasured  states.  The  design  of  such  an  ob¬ 
server  starts  with  a  preliminary  change  of  coordinates  such  that  the  output  y  consists 
of  the  first  p  entries  of  the  state  vector  x  =  [yT  x((]T .  In  the  new  coordinates,  the 
system  (5.1)  is 

ij  =  A\x0  +  Gi'y(H\y  +  H2x0)  +  Qi{y,  u ) 

x0  =  A2x0  +  G2'y(Hiy  + H-2x0)  + g2(y,u) ,  (5.29) 

where  the  linear  terms  in  y  are  incorporated  in  g\(y,u)  and  g2(y,u).  An  estimate  of 
x0  will  be  obtained  via  x  :=  x0  +  Ny ,  where  N  G  jR(n_P)xP  is  to  be  designed.  From 

(5.29) ,  the  derivative  of  x  is: 

X  =  (A2  +  NAJx  +  (G2  +  NGMH2X  +  -  H2N)y )  +  g(y,  u ) ,  (5.30) 

where  g{y,u)  :=  Ng\(y,u)  +  g2{y,u)  —  (A2  +  NA\)Ny.  In  this  ^-subsystem,  the 

output  injection  matrix  N  has  altered  the  A2  and  G2  matrices  of  the  x0-subsystem 

(5.29) . 

To  obtain  the  estimate 

Xo  =  X~Ny,  (5.31) 


67 


we  employ  the  observer 

{  =  (A2  +  N A\)x  +  (G2  +  NGMH2x  +  (ffi  -  H2N)y)  +  g(y,  u ) .  (5.32) 

From  (5.30)  and  (5.32),  the  dynamics  of  ea  :=  xQ  —  xQ  =  x  ~  X  are  governed  by 


ea  =  (A2  +  NAi)e0  +  (G2  +  NGi)[y(v0)  -  j{w0)\ , 


(5.33) 


where  vQ  :=  H2x+(Hi~H2N)y  and  wQ  :=  H2x+(Hi -H2N)y.  We  let  2  :=  v0-w0  = 
H2ea,  and  denote  <p(t,z)  =  7 (y0)  —  7 (w0).  Then,  the  nondecreasing  property  (5.3) 
implies  that  Zi<pi(t,  zi)  >  0  for  all  i  =  1,  •  •  •  ,  r.  Derivations  similar  to  those  in  Section 
5.1  yield  the  following  LMI  in  PQ,  P0N ,  v  and  A: 


(A2  +  N A\)T PQ  +  P0{A2  +  N A\)  +  vl  Pq(G2  +  NG\)  +  H ^  A 
(G2  +  A^Gi)P0  +  A^2  0 


<0.  (5.34) 


If  this  LMI  is  satisfied  with  a  matrix  Pa  =  Pj  >  0,  a  constant  u  >  0,  and  a  diagonal 
matrix  A  >  0,  then  it  is  not  difficult  to  show  that,  with  appropriate  modifications, 
Theorems  5.1  and  5.2  hold  for  the  observer  error  e0(t). 

To  illustrate  the  analog  of  Theorem  5.2,  we  design  a  reduced-order  observer  for 
the  system  (5.26). 

Example  5.2  The  system  (5.26)  is  of  the  form  (5.29)  with  xa  =  [x2  xf\T ,  6i{y,u)  = 
V2,  62{y,u)  =  [u  2 u]T, 


Ai  =  [l  0] 


A2  = 


1  1 
0  0 


G\  =0,  G2  = 


-1 

0 


Hi  =  0,  H2  =  [  1  0], 
(5.35) 


The  LMI  (5.34)  is  feasible  with  v  =  0,  and  a  solution  is 


Po  = 


1  0 
0  1 


N  = 


-2 

-1 


A  =  1. 


(5.36) 


Since  A2  +  N A\  is  Hurwitz,  the  resulting  observer 

X2  =  ~X2  +  X3  -  exp(y2  +  2 y)  +  (-2 y2  -  y  +  u) 

X3  =  -X2  +  ( ~y 2  -2 y  +  2 u) 

£2  =  X/2  +  2  y 

X3  =  X3  +  y 


(5.37) 


guarantees  (e2(t),  e3(t))  — ^  0  as  t  —>  00,  if  a  control  law  can  be  designed  to  ensure 
x(t)  e  Coo.  □ 
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5.3  Robustness  Against  Inexact  Modeling  of  Nonlinear¬ 
ities 

Thus  far  we  have  assumed  exact  knowledge  of  the  system  nonlinearities.  To 
analyze  the  effects  of  modeling  errors,  we  suppose  that  instead  of  (5.1),  the  system  is 

x  =  Ax  +  G['y(Hx)  +  A (Hx)fi(t)\  +  g(y,  u ) ,  (5.38) 

where  y(t)  is  a  bounded  disturbance.  Then,  the  nominal  observer  (5.4)  yields  the 
observer  error  system 

e  =  (A  +  LC)e  +  G  [y(u)  -  7 (w)  +  A (v)n(t)] ,  (5.39) 


where  v  =  Hx,  w  =  Hx  +  K(Cx  —  y ). 

In  Theorem  5.3  below,  we  characterize  nonlinearities  A(-)  for  which  the  observer 
(5.4)  guarantees  an  ISS  property  from  the  disturbance  fi(t)  to  the  observer  error  eft). 


Theorem  5.3  Consider  the  plant  (5.38)  and  the  observer  (5-4).  Suppose  xft)  exists 
for  all  t  >  0,  and  that  the  LMI  (5.11)  holds  with  a  matrix  P  =  PT  >  0,  a  constant 
v  >  0,  and  a  diagonal  matrix  A  >  0.  If,  for  each  i  =  1,  •  •  •  ,  r,  there  exists  a  class-K, 
function  cij(-)  such  that 


{a  -  b)(ji(a)  -  7 i(b )  +  A ,(a)/x]  >  -<Jj(|/x|)  Va,6,/x  G  1R , 
then  the  observer  error  e(t)  satisfies,  for  all  t  >  0, 


|e(t)|  <  «|e(0)|  exp(-/3t)  +  p\  sup  |/u(r)| 

\o  <T<t 

where  k  =  \J  X\ia(L<(j)^ ,  P  =  2A  —  ^ ,  and  the  ISS-gain  from  y(t)  to  eft)  is 


(5.40) 


(5.41) 


P(-)  =  K 


\ 


Vii- 

7/  ^ 


(5.42) 


1=1 


Proof:  We  use  V  =  eTPe  as  an  ISS-Lyapunov  function,  and  evaluate  its  derivative 
for  (5.39): 

r 

V  <  — ^|e|2  -  2  ^  A i  (vi  -  Wi)  [7 i(vi)  -  jiiwi)  +  A i(vi)n\ .  (5.43) 

i=  1 
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Substituting  (5.40),  we  obtain 


V<-v |e|2  +  2^AlCj.8(|/u|)  <  -2(3V  +  2^\iai(\n\) 


i= 1 


i= 1 


from  which  it  follows  that 


(5.44) 


V(t)  <  V (0)  exp(— 2/?i)  +  -  (  VAi  sup  <Tj(|//(r)|) 


This  yields 


(5.45) 


A i  sup  o-j(|/i(r)|) 

0 <T<t  ) 

(5.46) 

1  _  2k2 

P^min{P')  V 

□ 

The  ISS  property  established  by  Theorem  5.3  shows  that  e(t)  degrades  grace¬ 
fully  with  the  increase  in  the  magnitude  of  the  disturbance  /x(f).  As  fi{t)  vanishes, 
we  recover  the  convergence  result  of  Theorem  5.1.  The  dependence  of  admissible 
nonlinearities  A(-)  on  y(-)  is  characterized  by  (5.40).  For  example,  if  y(-)  is  cubic, 
then  A(-)  is  allowed  to  be  linear.  In  this  case,  (5.40)  is  satisfied  because 

(a  —  6)  [a3  —  b 3  +  a/j]  >  — /r2  (5.47) 

3 

holds  for  all  a,  b,  n  £  M,  due  to  the  identity  a3  —  b3  =  (a  —  b) (a2  +  ab  +  b2).  On  the 
other  hand,  (5.40)  does  not  hold  for  cubic  y(-)  and  quadratic  A(-).  To  see  this,  we 
evaluate  (a  —  6)  [a3  —  b3  +  a2/j]  with  b  =  a+  and  note  that,  for  any  fixed  [i  >  0,  the 
resulting  function  3  +  -A  +  \  —  afi  tends  to  —  oo  as  a  — ►  +oo,  thus  violating  (5.40). 

It  is  not  difficult  to  prove  the  analog  of  Theorem  5.3  for  the  reduced-order 
observer  (5.32),  in  which  case  e(f)  is  to  be  replaced  by  eQ(t). 


5.4  Summary 

An  observer  design  is  presented  for  systems  with  monotone  nonlinearities  in 
the  unmeasured  states.  The  new  design  represents  the  observer  error  system  as  the 
feedback  interconnection  of  a  linear  system  and  a  multivariable  sector  nonlinearity. 
The  issues  of  existence  and  the  evaluation  of  the  observer  matrices  K  and  L  satisfying 
the  circle  criterion  are  resolved  by  efficient  LMI  computations.  The  robustness  of  the 
new  observer  to  inexact  modeling  of  nonlinearities  is  analyzed,  and  ISS  bounds  are 
derived  within  which  the  observer  error  increases  with  an  increase  in  the  modeling 
error. 
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Chapter  6 

Feasibility  Conditions  for  the  Observer 
Design 


Thus  far,  the  feasibility  of  our  observer  design  was  left  to  be  resolved  by  iterative 
LMI  computations.  In  this  chapter  we  derive  structural  conditions  that  characterize 
feasibility.  We  first  show  that  rendering  a  linear  system  SPR  by  output  injection 
is  equivalent  to  rendering  its  dual  system  SPR  by  control.  Then,  we  characterize 
the  feasibility  of  this  dual  problem,  and  use  it  to  establish  necessary  and  sufficient 
feasibility  conditions  for  the  observer  design. 

The  duality  property  is  discussed  in  Section  6.1.  The  feasibility  of  the  dual 
problem  is  studied  in  Section  6.2.  In  Section  6.3,  we  give  the  feasibility  conditions 
for  the  observer  design,  and  illustrate  them  on  examples.  Finally,  in  Section  6.4,  we 
prove  that  the  feasibility  conditions  for  the  reduced-order  observer  are  the  same  as 
those  for  the  full-order  observer.  The  lengthy  proofs  are  given  in  Section  6.5. 

6.1  The  Dual  Problem 

The  observer  design  in  the  previous  chapter  relied  on  solving  the  LMI  (5.11), 
rewritten  here  as 


(A  +  LC)TP  +  P(A  +  LC)  <  0  (6.1) 

PG  +  (H  +  KC)t  A  =  0.  (6.2) 

Multiplying  the  inequality  (6.1)  from  both  sides  by  V  :=  P -1,  and  multiplying  the 
equality  (6.2)  from  the  left  by  V,  and  from  the  right  by  II  :=  A-1,  we  obtain 

V{A+BP)  +  {A  +  BP)TV  <  0  (6.3) 

p{g  +  B£)  +  nTn  =  o,  (6.4) 
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where 


A:=At,  B  :=  CT ,  £  :=  KT ,  T:=Lt ,  Q  :=  HT,  H:=GT.  (6.5) 

In  view  of  the  PR  Lemma,  (6.3)  and  (6.4)  mean  that  the  control  system 

x  =  Ax  +  Bu  +  Qw  (6.6) 

z  =  —Hx  (6.7) 

is  rendered  SPR  from  the  disturbance  w  to  the  output  f Iz  by  the  control  law 

u  =  Tx  +  £w.  (6.8) 

Thus,  our  observer  design  is  dual  to  a  control  problem  in  which  the  system  (6.6)-(6.7) 
is  to  be  rendered  SPR  by  state  feedback  and  disturbance  feedforward.  The  state 
feedback  matrix  T  is  the  dual  of  the  observer  gain  matrix  L,  and  the  feedforward 
matrix  £  is  the  dual  of  the  nonlinear  injection  matrix  K. 

This  dual  problem  encompasses  the  circle  criterion  control  design  of  Jankovic 
et  al.  [30].  For  the  system  in  Figure  6.1,  where  the  nonlinearity  w  =  —(j)(t,  z)  satisfies 
the  sector  property  z(j)(t,  z)  >0,  the  design  in  [30]  renders  the  linear  system  SPR  with 
respect  to  the  input  w  by  using  a  control  law  of  the  form  (6.8),  thus  achieving  global 
uniform  asymptotic  stability  via  the  circle  criterion.  This  means  that  T  and  £  are  to 
be  designed  such  that  (6.3)  and  (6.4)  hold  for  some  V  =  VT  >  0. 


Figure  6.1:  Circle  criterion  control  design. 


In  [30],  the  feasibility  issue,  that  is,  the  existence  conditions  of  a  control  law 
of  the  form  (6.8)  that  renders  the  resulting  closed-loop  system  SPR  from  w  to  z,  is 
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left  to  be  resolved  by  iterative  LMI  computations.  None  of  the  feedback  passivation 
conditions  in  the  literature  are  applicable  to  the  feasibility  problem  for  the  design  in 
Figure  6.1.  The  results  such  as  Fradkov  [16],  and  Kokotovic  and  Sussmann  [47],  deal 
with  ‘direct’  passivation  for  the  control  input  u,  and  not  with  ‘indirect’  passivation 
for  the  disturbance  input  w.  Recent  indirect  passivation  conditions  derived  in  the 
Tfoo-framework  by  Safonov  et  al.  [80],  Haddad  and  Bernstein  [22],  Sun  et  al.  [90], 
and  Turan  et  al.  [100]  assume  that  the  relative  degree  from  w  to  z  is  zero,  which  does 
not  hold  for  (6.6)-(6.7). 

6.2  Feasibility  of  the  Circle  Criterion  Design 

In  this  section  we  derive  structural  conditions  that  completely  characterize  the 
feasibility  of  the  circle  criterion  design.  This  result  will  be  used  in  the  next  section 
to  establish  feasibility  conditions  for  the  observer  design.  To  simplify  the  derivations, 
we  restrict  our  analysis  to  a  single  nonlinearity,  that  is,  Q  and  £  in  (6.3)-(6.4)  are 
scalars.  We  let  £1  =  1,  £  =  p.  and  analyze  the  existence  of  J~ .  p  and  V  =  VT  >  0 
satisfying 

(A  +  BT)tV  +  V{A  +  B£F)  <  0  (6.9) 

V{Q  +  PB)  +  HT  =  0.  (6.10) 

The  case  where  p  is  constrained  to  be  zero  is  of  separate  interest,  because  then  the 
control  law  (6.8)  can  be  implemented  without  the  exact  knowledge  of  the  nonlinearity 
w  =  —(j)(t,  z). 

After  a  change  of  coordinates  and  a  preliminary  state  feedback,  the  system 
(6.6)-(6.7)  is  represented  as 


£ 

=  A0£  +  E0yi  +  G0w 

(6.11) 

in 

=  2/2  +  9iw 

in 

=  2/3  +  92W 

(6.12) 

Vr 

=  u  +  grw 

z 

=  2/1, 

(6.13) 

where  £  £  Mn~r,  and  r  is  the  relative  degree  from  the  output  z  to  the  control  input  u. 
As  will  be  shown  later,  a  crucial  ingredient  of  the  feasibility  problem  is  the  existence 
of  a  matrix  Y  =  YT  >  0  satisfying 

A0Y  +  Y Aq  +  2 EQGl  +  2GqEq  <  0,  (6.14) 
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or,  equivalently, 


A0T  +  Y Aq  <  ~(E0  +  Gq)(Eq  +  Gof  +  (E0  -  G0)(E0  -  G0)T . 


When  Aq  is  Hurwitz,  the  existence  of  Y  is  immediate.  If  Aq  is  not  Hurwitz,  we 
decompose  the  ^-subsystem  (6.11)  into  three  subsystems, 


ki  —  Aq^  +  EqX\  +  GqW,  i  —  1,2,  3,  (6.15) 

such  that  ct(Aq)  C  C+,  <t(Aq)  C  C°,  ct(Aq)  c  C“  and  consider  U  =  UT ,  V  =  VT 
defined  by 

AlU  +  UAf  =  (El-GbiEl-Glf  (6.16) 

AlV  +  VAf  =  (E^+GM  +  Glf.  (6.17) 


Theorem  6.1  (p  =  0)  A  state  feedback  control  lawu  =  Ex  that  renders  (6.11)-(6.13) 
SPR  from  w  to  z  =  y\  exists  if  and  only  if 

9i  >  0,  <72  <0,  U-V>  -G^Gf,  (6.18) 

g  i 


and 


w*(Eq  -  Gl){El  -  Gl)Tw  >  iv* (Eq  +  Gl){El  +  Gl)Tw  (6.19) 


for  every  eigenvector  w  of  Aq 


□ 


The  proof  is  given  in  Section  6.5. 

As  a  corollary,  we  give  the  feasibility  conditions  for  the  case  p  /  0,  when  the 
control  law  u  =  Ex  +  pw  can  arbitrarily  assign  gr  =  gr  +  p  in  (6.12).  If  r  >  3,  the 
feasibility  conditions  for  p  0  are  the  same  as  in  Theorem  6.1,  because  they  do  not 
depend  on  gr.  However,  if  r  =  2,  then  <72  <  0  is  not  required  because  we  can  use  p  to 
satisfy  <72  <  0.  Likewise,  if  r  =  1  then  <71  >  0  is  no  longer  required.  Moreover,  since 
gi  >  0  can  be  arbitrarily  large,  U  —  V  >  is  replaced  by  the  less  restrictive 

condition  U  >  V . 


Corollary  6.1  ( p  /  0)  When  r  =  1,  a  control  law  u  =  Ex  +  pw  that  renders  (6.11)- 
(6.13)  SPR  from  w  to  z  =  y\  exists  if  and  only  if  U  >  V  and  (6.19)  holds  for  every 
eigenvector  w  of  AqT .  When  r  =  2,  g\  >  0  and  U  —  V  >  ^-GqGqT  are  required  in 
addition.  When  r  >  3,  all  the  conditions  of  Theorem  6.1  are  required.  □ 
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6.3  Feasibility  of  the  Observer  Design 

We  now  derive  necessary  and  sufficient  conditions  for  the  feasibility  of  the  ob¬ 
server  design.  For  the  system 

x  =  Ax  +  G^(Hx)  +  g(y,u)  (6.20) 

y  =  Cx, 

with  a  single  output  y  €  M,  and  a  single  nondecreasing  nonlinearity  yf),  we  determine 
when  the  LMI 

(A  +  LC)TP  +  P(A  +  LC)  +  vl  PG  +  {H  +  KC)T 
GtP  +  (H  +  KC )  0 

is  feasible.  When  u  =  0  is  allowed  as  in  Theorem  5.2,  the  feasibility  conditions  become 
cumbersome,  therefore  we  restrict  our  discussion  to  feasibility  with  v  >  0. 

We  omit  g(y,u)  and  the  linear  terms  in  y  from  the  right-hand  side  of  (5.1)  be¬ 
cause  they  do  not  affect  feasibility,  and  represent  the  resulting  system  in  the  following 
canonical  form 

V  =  !J\ 

in  =  <j2 

m  =  m 

(6.22) 

yr  =  nf  -  7(££  +  CTiyi  H - b  aryr) 

£  =  st, 

where  r  is  the  relative  degree  from  the  output  y  to  the  nonlinearity  7(-)-  Next,  we 
decompose  S,  II,  and  S  as 

'  Si  0  o' 

5=  o  s2  o  ,  n  =  [n1n2n3],  s  =  [s1s2s3],  (6.23) 

_  0  0  s3  _ 

where  cr(Si)  CC+,  cr (52)  CC°,  cr(53)  CCh  and  dehne  U  =  UT ,  V  =  VT  by 

SfU  +  USi  =  (Ili  -Si)T(n!  -Si)  (6.24) 

SfV  +  VS!  =  (Ili  +  Si)T(n!  +  Si).  (6.25) 

Using  Theorem  6.1,  Corollary  6.1,  and  the  duality  argument  in  Section  6.1,  it  is  not 
difficult  to  prove  the  following  result: 


<  0,  (6.21) 
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(6.26) 


Theorem  6.2  The  LMI  (6.21)  is  feasible  with  v  >  0  and  K 


0  if  and  only  if 


ar  >  0,  ar~ i  <  0, 


2  rp 

U  -V  >  — SfSi, 

Gy 


and 


™*(n2  -  s2)T(n2  -  t2)w  >  W*(U2  +  s2)T(n2  +  s2)™  (6.27) 

for  every  (possibly  complex)  eigenvector  w  of  S2. 

If  the  restriction  K  =  0  is  removed,  then  these  conditions  are  relaxed  as  follows: 
When  r  =  1,  the  LMI  (6.21)  is  feasible  with  v  >  0  if  and  only  if  U  >  V  and  (6.27) 
holds  for  every  eigenvector  w  of  S2.  When  r  =  2,  ar  >  0  and  U  —  V  >  are 

also  required.  When  r  >  3,  all  the  conditions  for  feasibility  with  K  =  0  above  are 

required.  □ 

Example  6.1  With  y  =  x\,  and  the  nondecreasing  nonlinearity  7(0:3)  =  0:3,  the 
system 

x\  =  xi  +  x2  +  xl 

x2  =  -x2  +  x3  (6.28) 

x3  =  u, 

is  of  the  form  (6.20).  Omitting  y  =  x\  and  u  from  the  right-hand  side  of  (6.28),  and 
using  the  change  of  variables  f2  =  —x3,  £3  =  x2  —  x3,  we  obtain 

y  =  -£2  +  £3  -  7(6) 

6  =  0  (6.29) 

£3  =  -£3, 

which  is  of  the  form  (6.22).  With  the  restriction  K  =  0,  the  observer  design  is  not 
feasible  because  ar  =  0.  When  the  restriction  K  =  0  is  removed,  the  design  is  feasible 
because  r  =  1,  the  eigenvalues  of  S  are  X2  =  0,  A3  =  —1,  and  (6.27)  holds  for 

the  £2- subsystem  with  S2  =  0,  II2  =  —  1,  S2  =  1.  A  solution  of  the  LMI  (6.21)  is 

L  =  [—0.6929  —  0.5621  —  0.4814]r,  K  =  —3.2417  and,  hence,  the  resulting  observer 
is 

x\  =  — 0.6929(.xi  -  y)  +  y  +  x2  +  (x3  -  3.2417(£i  -  y))3 

x2  =  — 0.5621(.ti  -  y)  -  x2  +  x3  (6.30) 

x3  =  — 0.4814(.ti  —  y)  +  u. 
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Example  6.2  Consider  again  the  system  (6.28),  but  with  —X2  in  the  second  equation 
replaced  with  X2  ■' 

x\  =  xi+x2  +  xl 

±2  =  X2  +  X3  (6.31) 

x3  =  u. 

With  y  =  x\,  the  omission  of  x\  and  u  from  the  right-hand  side,  and  the  change  of 
variables  =  X2  +  x3,  £2  =  —  x3,  we  get 

V  =  £1  +  6  -7(6) 

£1  =  6  (6-32) 

£2  =  0, 

where  the  eigenvalues  of  S  are  X\  =  1  and  X2  =  0.  The  observer  design  is  not  feasible 
even  with  K  /  0  because  U  >  V  does  not  hold  for  the  £1  -subsystem  with  S 1  =  1. 
IIi  =  1,  Si  =  0.  Thus,  the  sign  in  front  of  x 2  in  the  second  equation  of  (6.28)  is 
crucial.  □ 


6.4  Feasibility  of  the  Reduced-Order  Observer 

We  now  show  that  the  feasibility  conditions  for  the  reduced-order  observer  design 
are  the  same  as  those  for  the  full-order  observer. 

Theorem  6.3  Let  the  constant  u  >  0  and  the  diagonal  matrix  A  >  0  be  given.  Then, 
the  following  two  statements  are  equivalent: 

1.  There  exist  matrices  P  =  PT  >  0,  K  and  L  satisfying  the  full-order  observer 
LMI  (5.11). 

2.  There  exist  matrices  P0  =  Pj  >  0  and  N  satisfying  the  reduced-order  observer 
LMI  (5.34). 


Proof: 


(1  =>  2)  Suppose  the  full-order  observer  LMI  (5.11)  is  feasible  for  the  system  (5.29). 
Partitioning  P  and  L  as 


P  = 


Pi  Pi 
Pi  P3 


L  = 


(6.33) 
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and  substituting 


A  = 


0  Ai 

0  A2 


,  G  = 


G\ 

G2 


H  =  [H1H2],  and  C=[I  0]  (6.34) 


in  (5.11),  we  obtain 

(A  +  LC)tP  +  P(A  +  LG)  +  vl  = 


~k  ~k 

*  P3(A2  +  P^PjA,)  +  (A2  +  P-1P?A1)TP3  +  vl 


PG  +  {H  +  KC)tA  = 


=  0. 


p3(g2  +  p^1p^g1)  +  h^a 

Defining  PQ  :=  P3,  and  N  :=  P^lP2  ,  we  note  that  (6.35)  and  (6.36)  imply 

Po{A2  +  N  A\ )  +  (A2  +  NA\)tP0  +  z/J  <  0 
Po(G2  +  NGx)  +  Hi  A  =  0, 

which  is  the  reduced-order  observer  LMI  (5.34). 


(6.35) 

<  0 

(6.36) 


(6.37) 

(6.38) 


(2^1)  Suppose  (6.37)  and  (6.38)  hold  for  the  system  (5.29).  To  prove  the  existence 
of  P  =  PT  >  0,  K  and  L  satisfying  (5.11),  we  rewrite  system  (5.29)  in  y  and 
X  =  x0  +  Ny  coordinates,  so  that 


’  u 

A\ 

,  G  = 

G\ 

7/2 

A2  +  NAi 

G2  +  NG\ 

A  +  LC  = 


H  +  KG  =[H-i-  H2N  +  K  H2], 
where  L\  :=  L\  —  A\N  and  L2  =  L2  —  ( A2  +  N A\)N .  We  let 


(6.39) 


P  = 


I  0 

0  Pn 


(6.40) 


and  obtain 

(A  +  LC)TP  +  P(A  +  LG)  = 


L\  +  L\  A\  +  l\P0 


PG  +  (H  +  KC)t  A  = 


PoL2  +  AJ  P0(A2+NA1)  +  {A2+NA1)tP0 

G\  +  {H\  -  H2N  +  K)T A 
Po(G2  +  NG\)  +  H2  A 


(6.41) 
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In  view  of  (6.37)  and  (6.38),  the  choice 

Li  =  ~I,  L2  =  -Pq1AJ,  K  =  -H1+H2N-A-1Gj  (6.42) 

results  in 

(A  +  LC)TP  +  P(A  +  LC)  +  uI<  0 
PG  +  {H  +  KC)t  A  =  0, 

which  is  the  full-order  observer  LMI  (5.11).  □ 


6.5  Proof  of  Theorem  6.1 


We  prove  the  theorem  in  two  steps.  In  Step  1  we  prove  that  a  state  feedback 
control  law  u  =  Px.  rendering  (6.11)-(6.13)  SPR  from  w  to  z,  exists  if  and  only 
if  g\  >  0,  g2  <  0,  and  Y  =  YT  >  j^GqGq  can  be  found  such  that  (6.14)  holds. 
In  Step  2,  we  show  that  Y  =  Yr  >  satisfying  (6.14)  exists  if  and  only  if 

U  —  V  >  -^GqGqT  and  (6.19)  holds  for  every  eigenvector  of  A^ . 


Step  1  -  Necessity:  Suppose  u  =  Px  renders  (6.11)-(6.13)  SPR  from  w  to  z  =  y i, 
let  X  :=  V _1,  and  rewrite  (6.9)-(6.10)  as 


(A  +  BP)X  +  X(A  +  BPf  <  0 
XUT  +  Q  =  0 

which,  for  system  (6.11)-(6.13),  yields 


^0 

Go 

k 

k 

^  a 

II  11 

.  Olx  (n— r) 

[Go  gi  ■  ■ 

-  1  0  •• 

•  9r] 

■  =►  X  = 

rT 

Ctq 

k 

9i 

k 

9r 

k 

k 

9r 

k 

k 

(6.43) 

(6.44) 


(6.45) 


Then,  g\  >  0  because  X  >  0,  and  g2  <  0  because,  from  (6.11)  and  (6.45), 


(A  +  BP)X  +  X(A  +  BP)t 


+  X0Al  +  EqGq  +  GqEq  -k  -k  -k 
~k  2g2  *  'k 

'k  'k  it 


<0. 


-k  k  k  k  J 

(6.46) 
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This  also  shows  that  Y  =  2Xq  satisfies  (6.14).  Then,  Y  >  Aq0qT  follows  from  (6.45) 
because 


Go 

G'o  9i 


>0, 


and  the  Schur  complement  of  g\  is  Xq  — -GqGq  >  0. 


(6.47) 


Step  1  -  Sufficiency:  For  r  =  1,  the  equations  (6.11)-(6.13)  are 


£  =  A£  +  E0yi  +  G0vj 
Vi  =  u  +  giw. 


For  this  system  7i  =  [0lx(n_r)  —  1],  and  (6.44)  is  satisfied  by 


X  = 


where  Y  is  as  in  (6.14).  Substituting 


T  = 


Go 

GT, 

9i 

y-1 

K 

>  0 


2  J 


K  >  0 


in  (6.43),  we  obtain 


(6.48) 

(6.49) 


(6.50) 


(6.51) 


(A  +  BT)X  +  X(^l  +  BT)t  = 

\AqY  +  iy  Aq  +  EqGq  +  GqEq 
(A0Gq  +  giE0)T 


(6.52) 

A)  Go  +  giEo 
-k(9i-2G^Y-1G0) 


From  (6.50),  g\  —  2G^Y  1Gq  >  0  and,  hence,  the  right-hand  side  of  (6.52)  is  rendered 
negative  definite  by  selecting  k  >  0  sufficiently  large. 

For  r  >  2,  we  need  the  following  lemma  proved  at  the  end  of  the  section: 


Lemma  6.1  With  input  w  £  M  and  output  z  G  1R,  the  system 

r)  =  My  +  Nz 
z  =  Ly  —  az  +  gw 

is  SPR  if  and  only  if  a  >  0,  g  >  0,  and  the  C2-gain  of  the  system 

1 


r )  =  (  M  +  - — NL  )  g  +  Nv 

2a 


y  =  Lr) 

from  input  v  to  output  y  is  7  <  2a. 


(6.53) 

(6.54) 


(6.55) 


□ 
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*  =  2,  ■  ■  ■ 


With  the  change  of  coordinates 


Z  =  yi,  £  =  £ - G0yi,  yi  =  yi-—yi, 

9i  9i 

the  system  (6.11)-(6.13)  is  rewritten  as 


,r,  (6.56) 


£  —  A)£ - Goy2  +  f — AqGq  +  Eq  —  \  z 

9i  \9i  9\  ) 

-  -  9i  -  9i92  .  0  . 

Vi  =  Vi+1 - 92 - o“2,  i  =  2,- ■ ■  ,  r  —  1, 

9i  91 

'i.  9r  _ 

yr  =  U - y2 - —  Z 

91  91 

-  9 9 

2  =  2/2  +  —  z  +  giw, 

9i 


(6.57) 

(6.58) 

(6.59) 

(6.60) 


which  is  (6.53)-(6.54)  with  g  =  {£r,y2,---  ,9r)T ,  Lr]  =  y2,  a  =  ,  and  g  =  g1. 

Since  <71  >  0  and  g2  <  0,  we  conclude  that  a  >  0  and  g  >  0  as  in  Lemma  6.1.  The 
system  (6.55)  in  Lemma  6.1  is 


where 
Ro  = 


£ 


9i 


9r 

y 


^4o£  +  Rq92  +  Q  ov 

9i92 


m+\ 


9i  - 
^—92 
2gi 


9i 


i  =  2,  •  •  •  ,r  —  1, 


9r  _ 
U  -  —y2 

2<?i 

92, 


9r9  2 
2  r> 
91 


(6.61) 

(6.62) 

(6.63) 

(6.64) 


( ^~A0Go  +  i^E0  +  ^-Go]  and  Q0  =  ( -A0G0  +  E0-  % G0 )  . 

\2g2  2g2  2g{  )  \gx  g{  ) 

(6.65) 


Next,  we  show  that  there  exists  a  state  feedback  control  law  for  u  that  assigns  an 
£2-gain  7  <  2 a  from  input  v  to  output  y  =  y2.  It  was  shown  by  Isidori  [27,  Section 
13.2]  and  Chen  [10,  Chapter  5]  that  a  control  law  ensuring  7  <  2 a  exists  if  and  only 
if  there  exists  a  matrix  Z  =  ZT  >  0  satisfying 

AqZ  +  ZAq  +  —  9  QqQq  —  RqRo  <  0-  (6.66) 

\^a) 

Substituting  a  =  —  Qo  and  Rq  from  (6.65),  it  is  not  difficult  to  verify  that  (6.66) 
holds  for 
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(6.67) 


where  Y  satisfies  (6.14).  Moreover,  Z  >  0  because  g 2  <  0  and  Y  >  AGqGq  .  This 
means  that  there  exists  u  =  Ex  such  that  7  <  2a,  thus  rendering  (6.11)-(6.13)  SPR 
from  w  to  z,  by  Lemma  6.1. 


Step  2  -  Necessity:  Suppose  that 

'  Yu  Y12  Y13  I 

Y  =  Yt=  Y21  Y22  Y23  >  —G0Gq  (6.68) 

_  Y31  y3  2  y33  \  91 

satisfies  (6.14).  Then,  Yu  >  ■^■Gi0GqT ,  i  =  1,2,3,  and 

A*0Y,  +  YnAf  <  -(E*  +  G*)(i$  +  Gq)t  +  (£*  -  G^Ei  -  Gt0)T.  (6.69) 

Because  —  Aq  is  Hurwitz,  it  follows  that  there  exists  Y\\  >  Yu  >  ■^-GqG'q^  satisfying 

AhYu  +  YuAf  =  -{El  +  Gq){Eq  +  Gl)T  +  (^  -  Gj)^1  -  G(,)t.  (6.70) 

From  (6.16)  and  (6.17),  the  solution  of  (6.70)  is  Y\  \  =  U  —  V,  and  U  —  V  > 
follows  from  Yu  >  ^-GqGqT. 

For  (6.69)  with  i  =  2,  the  following  lemma  due  to  Scherer  [82,  Theorem  4]  shows 
that  (6.19)  holds  for  every  eigenvector  w  of  A^T: 

Lemma  6.2  Let  all  the  eigenvalues  of  A  be  on  the  imaginary  axis,  and  let  Q  be  an 


arbitrary  matrix.  Then,  there  exists  X  =  XT  >  0  such  that 

AX  +  XAt  <  Q  (6.71) 

if  and  only  if 

w*Qw  >  0  (6.72) 

for  every  eigenvector  w  of  AT .  Moreover,  if  (6.72)  holds,  then  for  any  Xq,  there 
exists  X  >  Xq  satisfying  (6.71).  □ 

Step  2  -  Sufficiency:  Let  Y33  =  Y^  >  0  be  the  solution  of  A'lY33  +  Y33AqT  =  —kl, 
and  denote 


Al  0 
0  Al 


,  G0  = 


(6.73) 


If  there  exists  Y  =  YT  >  ^-GqG^  such  that 
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(6.74) 


AqY  +  Y Aq  <  —(Eq  +  Gq)(Eo  +  Gq)t  +  (Eq  —  G0)(E0  —  Gq)t , 


then 


Y  = 


Y  0 
0  Y>;> 


(6.75) 


satisfies  (6.14),  and  Y  >  A-GoGq  for  large  enough  k  >  0.  To  prove  (6.74),  we  use 

AlYn  +  YuAlT  <  -(El  +  GbXEl+Glf  +  iEl-GM-Glf  (6.76) 
A20Y22  +  Y22AlT  <  -(E2  +  G2)(E2  +  G2)T  +  (E2-G2)(E2-G2)T,  (6.77) 

where  Aq^qA  <  yn  <  Y\  \  =  JJ  —  V  because  —  is  Hurwitz,  and  Y22  =  L22  >  0  in 
view  of  Lemma  6.2  and  (6.19).  Then 


satisfies 


Y  = 


iii  il2 
y£  y22 


(6.78) 


+  Ti0  +  (E0  +  Go)(^o  +  Gof  ~  (E0  -  G0)(E0  -  G0)T  = 

Qi  =  AqYh  +  YuAf  +  (El  +  Gq)(Eq  +  Gl)T  -  (E^  - 

Q2  =  A2Y22  +  Y22Af  +  (E2  +  G2)(E2  +  G2)T  -  (E2  -  G2)(E2 


Q 1  0 
0  Q2 

Gl)T  <  0 

G20)T  <  0, 


and  Y\2  satisfies  the  corresponding  Sylvester  equation.  Finally,  Y  >  Aq^qT  follows 
because  Y\\  >  Aq^qA  and,  from  Lemma  6.2,  Y22  can  be  selected  sufficiently  large. 
□ 


Proof  of  Lemma  6.1 


Let  the  system  (6.53)-(6.54)  be  SPR,  that  is,  let  P  =  PT  >  0  exist  such  that 


where 


AtP  +  PA  <  0 
PB  =  CT , 


A  = 

M  N 

,  B  = 

0 

,  c  = 

0 

L  —a 

9 

1 

(6.79) 

(6.80) 


(6.81) 
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6.6  Summary 

We  have  derived  necessary  and  sufficient  conditions  for  the  feasibility  of  the 
observer  design.  The  result  is  established  by  showing  that  the  observer  design  is 
dual  to  a  control  design  in  which  a  linear  system  is  rendered  SPR  with  respect  to  a 
disturbance  input.  We  have  proved  that  the  feasibility  of  the  reduced-order  observer 
design  coincides  with  that  of  the  full-order  observer. 
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Chapter  7 

Robust  Output-Feedback  Design 


As  with  any  other  observer,  the  main  purpose  of  the  new  observer  is  stabiliza¬ 
tion  and  tracking  in  conjunction  with  a  control  law  u  =  a(y,x,t)  which  employs  the 
available  measurement  y,  the  state  estimate  x  and,  possibly,  exogenous  signals.  Be¬ 
cause  the  separation  principle  does  not  hold  for  nonlinear  systems,  our  observer  is 
to  be  used  with  control  laws  that  guarantee  boundedness  of  the  states  for  bounded 
observer  errors.  This  design  strategy  is  detailed  in  Section  7.1,  and  illustrated  on 
an  analytical  example.  In  Section  7.2  we  study  the  effect  of  unmodeled  dynamics, 
and  propose  a  small-gain  design  for  robust  output-feedback  control.  This  small-gain 
design  is  illustrated  on  the  jet  engine  compressor  example  in  Section  7.3. 

7.1  Observer-Based  Control  Design 

We  view  the  state  observer  error  e  =  x  —  x  as  a  disturbance  acting  on  the 

system  (5.1)  through  the  control  law  u  =  a(y,x  —  e,t)  and  require  that  the  control 

law  guarantee  one  of  the  following  properties: 

No  Finite  Escape  (NFE)  Property: 

e(t)  €  Coo  =►  x(t)  G  Clc  ,  Vs(0)  G  Mn.  (7.1) 

Bounded  Error  -  Bounded  State  (BEBS)  Property: 

e(t)  G  £oo  =>  x(t)  G  £oo ,  Vx(0)  G  Mn.  (7.2) 

If  u  =  a(y,x,t )  is  a  NFE  control  law  for  the  system  (5.1),  then  the  assumption  that 
x(t)  G  C ^  made  in  Theorem  5.1  holds.  To  see  this,  suppose  that  the  maximal  interval 
of  existence  is  finite.  On  this  interval,  e{t)  is  bounded  because  of  (5.16).  If  x{t)  were 
unbounded,  this  would  contradict  the  NFE  property  (7.1),  thus  x(t)  exists  for  all 
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t  >  0.  Likewise,  if  u  =  a(y,x,t )  is  a  BEBS  control  law,  then  x(t)  £  Coo,  as  required 
in  Theorem  5.2.  To  guarantee  x(t)  £  C%q  for  the  system  (5.38),  the  NFE  property  is 
to  hold  for  all  bounded  disturbances  pit). 

The  BEBS  property  incorporates  the  NFE  property.  Requiring  these  properties 
is  meaningful  because  BEBS  control  laws  have  already  been  designed  for  classes  of 
nonlinear  systems  by  Freeman  and  Kokotovic  [18],  Krstic  et  al.  [53],  and  Marino  and 
Tornei  [63].  We  now  illustrate  the  use  of  our  observer  for  output-feedback  design  in 
conjunction  with  the  observer  backstepping  design  of  Krstic  et  al.  [53]. 

Example  7.1  For  the  system 


XI 

=  X2+  x\ 

X2 

=  X2  +  X3  -  exp(x2)  +  u 

(7.3) 

X3 

=  2  u , 

studied  in  Example  5.1,  suppose  that  the  output  y  =  x\  is  required  to  track  yd(t), 
where  yd(t),  i/d(t),  and  yd{t)  are  known,  continuous  and  bounded  on  [0,  oo).  Using  the 
observer  (5.28),  we  design  a  BEBS  control  law  u  =  a(y,x,t)  that  ensures  asymptotic 
tracking  and  boundedness  of  the  states.  Since  the  system  (7.3)  with  output  y  has 
relative  degree  two,  we  apply  two  steps  of  observer  backstepping  as  in  Krstic  et  al. 
[53,  Section  7.1]. 

Step  1.  We  let  Ci  ■=  y  —  yu,  substitute  X2  =  x2  +  e2  and  x\  =  y  in  the  first  equation 
of  (7.3),  and  obtain 

Cl  =  x2  +  y2  -  yd  +  ■  (7.4) 

For  X2,  we  design  the  ‘virtual’  control  law 

ao(y,yd,yd)  =  -ci(i-y2  +  i/d,  >  0,  (7.5) 

which  results  in 

Ci  =  -ci Ci  +  C2  +  e2  ,  (7.6) 


where  (2  :=  x2  -  a0(y,  yd,  yd) . 

Step  2.  Using  the  observer  equation  (5.28),  we  get 


(2  =  11  +  g(y,  yd,  yd,  tid,  x)  -  ^e2 

dy 


(7.7) 
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where 


g(y,yd,yd,yd,z )  =  x2  +  x3  -  exp(x 2  -  Ob  -  y))  -  4(x  i  -  y) 

dao  .  dao ..  dao 


-~^—yd  -  jr^Vd  -  -jr-  [x2  +  y  )  •  (7.8) 

oyd  oyd  dy 


We  let  c2  >  0,  d  >  0,  and 

u  =  afy,  x,  t )  =  -Ci  -  C2C2  -  d 
which  yields 

(2  =  ~Ci  —  C2C2  -  d 


dao 

dy 


(2  -  g(y,yd,Vd,yd,x) ,  (7.9) 


da0\  dao 

-5—  C2  -  -s— ' e2 

dy  /  dy 


(7.10) 


We  yrcwe  that  (7.9)  is  a  BEBS  control  law  with  the  Lyapunov  function  W (Ci,  C2)  = 
^Ci  +  5C2  -  Differentiation  along  (7.6),  (7.10),  and  completion  of  squares  yield 


w  ^  ~ "^"Ci  —  c2Cl  +  —  ~c°w  ^0  e2  >  (7.ii) 

where  cq  =  min{ci,2c2}  and  do  =  ^  From  (7.11),  boundedness  of  eft)  guar¬ 

antees  that  Ci,  C2  and,  hence,  x\  and  X2  are  bounded.  To  see  that  X3  is  also  bounded, 
we  define  the  new  variable  i]  :=  x 3  —  2x2,  governed  by 


r]  =  —  2rj  —  6x2  +  2exp(x2) .  (7.12) 

Since  x 2  is  bounded,  (7.12)  guarantees  boundedness  ofr),  therefore  X3  is  also  bounded. 

The  conditions  of  Theorem  5.2  being  satisfied,  the  observer  (5.28)  guarantees 
eft )  — ►  0  as  t  — »  00.  From  (7.11),  we  see  that  eft)  — >  0  ensures  ( (t)  — >  0  as  t  — >  00. 

Since  (1  =  y  —  yd,  we  conclude  that  the  observer-based  control  law  (7.9)  guarantees 
yft)  — >  ydft)  as  t  — >  00.  □ 


7.2  Robust  Output-Feedback  Stabilization 

We  now  consider  the  problem  of  output-feedback  stabilization  for  the  locally 
Lipschitz  system 

x  =  Ax  +  Gfy(Hx)  +  A(Hx)fj]  +  g(y,  u) ,  y  =  Cx,  (7-13) 

£  =  q(£,hfx))  (7.14) 

g  =  p{£,Kx)), 
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where  7(-)  and  A(-)  are  as  in  (5.40),  and  the  ([-subsystem  (7.14)  represents  unmodeled 
dynamics.  This  formulation  extends  the  applicability  of  our  observer  because,  if  there 
is  a  subsystem  to  which  the  observer  is  not  applicable,  it  can  be  treated  as  unmodeled 
dynamics. 

The  unmodeled  dynamics  (7.14)  are  assumed  to  possess  the  following  input-to- 


output  stability  (IOS),  and  ISS  properties: 

IM0I  <  max|^(^(0)|,£),pM/l 

(  SUP  IMzCO)!")  } 
\0<r<t  /  J 

(7.15) 

|£(0I  <  max  |/%(|£(0)|,£),pc/l  J 

sup  |/i(z(t))|')  ]  , 
,0<T<t  /  J 

(7.16) 

where  /3M(- ,  •),  /%(•,  ■)  are  class-/C£  functions,  and  p^h(-),  P^h(-)  are  class- /C  functions. 

The  ISS  property  of  the  observer  error  (5.41)  is  now  rewritten  in  the  ‘max’  form 
of  Teel  [93],  which  is  suitable  for  the  small-gain  analysis  we  will  employ.  Using  the 
fact  that  for  each  p(-)  E  /C0 0  ,  Va,  b  >  0,  a  +  b  <  max{(7  +  p)(a),  ( I  +  p_1)(6)},  where 
/(•)  =  (•)  represents  the  identity  function,  we  get 

|e(£)|  <  max  j/3e(|e(0)|,£),peM  ^supJ/i(r)|H  .  (7.17) 


h(x) 


Figure  7.1:  Closed-loop  system  with  observer  feedback. 

To  prepare  for  a  small-gain  design  of  u  =  a(y,x),  we  represent  the  closed-loop 
system  with  the  block-diagram  in  Figure  7.1.  The  gains  of  the  unmodeled  dynamics 
block  and  the  observer  error  block  are  pAt/l(-)  and  pe/i(-),  respectively.  Let  phe(-)  and 
Phfi  ( ■ )  be  the  plant  gains  from  e  and  p  to  h(x),  respectively.  The  task  for  the  control 
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law  u  =  a(y,x )  is  to  render  phe{' )  and  Ph.^(-)  small  enough  for  the  inner  loop  gain 
Phn  °  Pfihi-),  and  the  outer  loop  gain  phe  o  pefl  o  p^hfi)  to  satisfy,  for  all  s  >  0, 

Phn  0  Pnh(s)  <  s  (7.18) 

Phe0  Pen0  P/xh(s)  <  s  •  (7.19) 


Then,  GAS  of  the  closed-loop  system  will  be  guaranteed  as  in  the  nonlinear  small-gain 
theorem  of  Teel  et  al.  [39,  93]. 

Theorem  7.1  Consider  the  system  (7.13)-(7.14),  in  which  y(-)  and  A(-)  satisfy 
(5.40),  and  the  £- subsystem  satisfies  (7.15)  and  (7.16).  Suppose  that  the  observer 

x  =  Ax  +  L(Cx  —  y)  +  G^(Hx  +  K(Cx  —  y))  +  g(y,  u)  (7.20) 

is  such  that  the  LMI  (5.11)  holds  with  a  matrix  P  =  PT  >  0 ,  a  constant  v  >  0,  and 
a  diagonal  matrix  A  >  0.  If  the  control  law  u  =  a(y,x)  guarantees 

|  h(x(t))\  <  max  <  (3h(\x(0)\,t),  php  (  sup  \p{r)\ 

l  Vo <T<t 

\x(t)\  <  max  <  (3x(\x(0)\,t),  pXfM  (  sup  \p{r)\ 

l  Vo <T<t 

where  PhnO  and  Ph.e()  satisfy  (7.18)  and  (7.19),  respectively,  then  the  origin  of  the 
closed-loop  system  (7.13),  (7.14),  (7.20)  is  globally  asymptotically  stable. 


phe  f  sup  \e(T)\]  [>(7.21) 

Vo <T<t  J 


,  pxe  (  sup  |e(r)|  )  [>(7.22) 

\0<T<t 


Proof:  Using  the  notation  \\h\\t  :=  sup0<r<t  |h(x(r))|  for  each  t  in  the  maximal 
interval  of  existence  [0,  tf),  and  substituting  (7.15)  in  (7.17),  we  get 

||e||t  <  max  {/3e(|e(0)|,  0),  peix  o  /?M(|£(0)|,  0)|  pefl  o  p^h(\\h\\t)  }  .  (7.23) 

We  then  substitute  (7.23)  and  (7.17)  in  (7.21),  and  obtain 

\\h\\t  <  max{4(|x(0)|,  |£(0)|,  |e(0)|),  php  o  p^h{\\h\\t),  phe  0  Pen  0  Pnh(\\h\\t)}  , 

(7.24) 

where 

*fc(|z(0)U(0)|,|e(0)|)  :=  (7.25) 

max{/^(|x(0)|,0),  phli  o  ^(|C(0)|,  0),  phe  o  /3e(|e(0)|,  0),  phe  o  pefl  o  /?M(|£(0)|,  0)}  . 

The  substitution  of  the  small-gain  conditions  (7.18)  and  (7.19)  in  (7.24)  results  in 

||%<4(|x(0)|,|e(0)|,|e(0)|),  (7.26) 
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and,  from  (7.15)  and  (7.23),  it  is  not  difficult  to  derive  functions  6 M  and  5e  of  the 
initial  conditions  such  that  \\p>\\t  <  <5/t  and  ||e||t  <  5e.  Then,  in  view  of  (7.16)  and 
(7.22),  we  can  find  a  class -1C  function  fi(-)  such  that,  for  each  t  £  [0,  t /) , 

\(x(t),£(t),e(t))\  <  fi(|(x(0),f(0),e(0))|).  (7.27) 

Since  the  right  hand  side  of  (7.27)  is  independent  of  t,  we  conclude  that  tf  =  oo,  and 
the  origin  is  Lyapunov  stable. 

To  prove  GAS,  the  remaining  task  is  to  prove  convergence  of  the  solutions  to 
the  origin.  To  this  end,  we  denote  ||/i||0  :=  lirnsup^^  \h(x(t))\  as  in  Teel  [93],  and 
note  from  (7.27)  that  \\h\\af  \\/j,\\ a  and  ||e||a  are  finite.  Then,  from  (7.15),  (7.23)  and 
(7.21), 

Ma  <  P„h(\\h\\a )  (7.28) 

||e||a  <  Pen  0  Pnh(\\h\\a)  (7.29) 

IHIa  <  max-{Phn(\\P-\\a)i  Phe(\\e\\a)}  •  (7.30) 

The  substitution  of  (7.28)  and  (7.29)  in  (7.30),  and  the  use  of  (7.18)  and  (7.19)  yield 
1 1 || a  =  || e || a  =  \\h\\a  =  0,  thus  proving  the  convergence  of  p(t),  e(t)  and  h(x(t ))  to 
zero.  Then,  the  ISS  conditions  (7.16)  and  (7.22)  imply  (x(t),£(t),e(t))  — >  0.  □ 

7.3  Design  Example 

An  axial  compressor  model,  which  has  been  the  starting  point  for  jet  engine 
control  studies,  is  the  following  single-mode  approximation  of  a  PDE  model  due  to 
Moore  and  Greitzer  [65], 


<p  = 

-^P+\cP+l--l-{cP  +  lf -3(<P  +  1)R 

(7.31) 

Ip  = 

jp((P  +  l-u) 

(7.32) 

R  = 

crR{-2(p  -  (p2  -  R) ,  i?(0)>0, 

(7.33) 

where  4>  and  ip  are  the  deviations  of  the  mass  flow  and  the  pressure  rise  from  their  set 
points,  the  control  input  u  is  the  flow  through  the  throttle,  and,  a  and  f3  are  positive 
constants.  This  model  captures  the  main  surge  instability  between  the  mass  flow  and 
the  pressure  rise.  It  also  incorporates  the  nonnegative  magnitude  R  of  the  first  stall 
mode. 

Krstic  et  al.  designed  a  state  feedback  GAS  control  law  in  [53,  Section  2.4],  and 
later  replaced  it  by  a  design  using  ( j>  and  ip  in  [52],  With  a  high-gain  observer,  Isidori 
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[27,  Section  12.7],  and  Maggiore  and  Passino  [58],  obtained  a  semiglobal  result  using 
the  measurement  of  ip  alone.  With  y  =  ip,  we  will  now  achieve  GAS  for  (7.31)-(7.33). 
The  exact  observer  cannot  be  designed  because  of  the  nonlinearities  cpR  and  <p2R. 
However,  the  {<p,  ^-subsystem  (7. 31), (7. 32)  contains  the  nondecreasing  nonlinearity 
(4>  +  l)3,  and  is  of  the  form  (7.13)  with  disturbance  y  =  R.  This  suggests  that  we 
treat  the  i?-subsystem  (7.33)  as  unmodeled  dynamics  and  apply  the  design  of  Section 
7.2. 

First,  we  prove  that  y  =  R  satisfies  the  IOS  property  (7.15)  with  h(x)  =  cp  as  the 
input.  With  V  =  R2  as  an  ISS-Lyapunov  function,  R  >  2.1\<p\  implies  V  <  — 0.09cri?3, 
because  R(t)  >  0  for  all  t  >  0.  This  means  that  (7.15)  holds  with  the  linear  gain 

Pm(')  =  2-!(-) ,  (7-34) 

and,  since  y  =  £  =  R,  the  ISS  property  (7.16)  is  also  satisfied. 

To  design  the  reduced-order  observer  of  Section  5.2  for  the  (cp,  ^-subsystem,  we 
let  x  =  <p  +  Nip,  and  obtain 

X  =  Q  +  |Q  X  -  \(x  -  Nip  +  i)3  -  3(x  -Nip  +  1  )R  +  g{ip,  u) ,  (7.35) 

where 

g(ip,  u)  :=  -  Q  +  |Q  Nip  -  ip  +  1  +  |^(1  -  u) .  (7.36) 

The  resulting  observer  is  the  scalar  equation 

x  =  + -p)  x- Nip  +  if  +  e(ip,u) 

t  =  X-Nip.  (7.37) 

For  its  implementation,  the  LMI  (5.34)  is  satisfied  by  selecting  N  such  that 

k  :=  -  ^  +  /jQ  >  O'  (7-38) 

To  prove  the  ISS  property  (7.17)  for  the  observer  error  =  <p  —  <p,  we  employ  the 
ISS-Lyapunov  function  Ve  =  e|,  and  evaluate  its  derivative  for 

e<f,  =  -ke^  -  -(a3  -  b3  +  6ai?) ,  (7.39) 

where  a  :=  x  ~  Nip  +  1  and  b  =  \  —  Nip  +  1.  Employing  the  inequality  (5.47),  and 
substituting  a  —  b  =  e^,  we  get 

Ve  <  -2 ke\  +  12 R2  ,  (7.40) 
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from  which  |e^|  >  y  ^-\R\  implies  Ve  <  —0.03 ke^  and,  hence,  the  ISS  property  (7.17) 
holds  with  the  linear  gain 

PeM  =  ]/*£(  ■)■  (7.41) 

Clearly,  the  gain  pefi{-)  can  be  rendered  as  small  as  desired  by  making  k  sufficiently 
large  with  a  choice  of  N  in  (7.38). 

We  are  now  ready  to  design  a  control  law  as  in  Theorem  7.1.  Noting  that  (7.31)- 
(7.32)  is  in  strict  feedback  form,  we  apply  one  step  of  observer  backstepping.  For  -0, 
we  design  the  virtual  control  law  =  ci4>-  Denoting 

uj  ■=  ip  -  C\<j)  =  —  ci</>  +  cie«£,  (7.42) 


we  rewrite  (7.31)  as 

3  1 

4>  =  —ci 4>  —  -(j)2  —  - 4 i3  —  3 4>R  —  uj  —  3 R  +  cie^  .  (7.43) 

The  substitution  of  (7.37)  in  (7.42)  yields  uj  =  (1  +  Nci)iJj  —  ciy,  and,  from  (7.32) 
and  (7.37), 

^  =  1  +fpCl  <t>  +  ^2  (X  ~  u)  +  r(^  i  (7-44) 

where  T( (j),ip)  :=  Clip  +  c\k<j>  +  ^-(^  +  l)3  —  -  Then,  the  control  law 

u  =  1  +  (1  +  Nc\)4>  +  f32(c2UJ  +  T(0,  ijj))  (7.45) 


is  implementable  using  the  signals  i/j  and  <f>,  and  results  in 

1  +  iVci 

UJ  =  —C2UJ  H - ^2 - e</>- 


(7.46) 


The  remaining  task  is  to  select  the  design  parameters  ci  and  C2  such  that  (7.21) 
and  (7.22)  are  satisfied.  For  the  ISS-Lyapunov  function  W((j),uj)  :=  \4>2  +  \uj2 ,  the 
inequalities  — 1<^>3  <  |<^>2  +  —(j)uo  <  \ 4> 2  +  o2,  —307?  <  |c j)2  +  R2,  —3 4>2R  <  0 
(because  R(t)  >  0),  cife#  <  ^4>2  +  and  1+^C1  u2  +  yield 


IF  <  - 


£i 

2 


31 


C2 


(1  +  JVC,)J  V  +  tf  +  ceJ.  (7.47) 


2/34ci 


We  let  c  >  0,  and  select  ci  and  C2  to  satisfy 


ci 


31 


°2  ~  2 


1  (1  +  Nc\Y 


2/34ci 


>  c, 


(7.48) 
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so  that 


31 

W  ^  — c(02  +  u}“)  +  i?2  +  (2c  H — — )e;}  )  (7.49) 

from  which  (7.22)  follows  for  x  =  ((j),ip).  For  h(x)  =  4>  and  p  =  R,  we  now  compute 
the  gains  Ph.fi(-)  and  ph.e{')  in  (7.21).  Using  the  fact  that  for  each  constant  0  >  0, 
a  +  b  <  max  {(1  +  0_1)a,  (1  +  0)6}  for  all  a,  b  >  0,  we  obtain 

31 

W  <  — c((/>2  +  cj2)  +  i?2  +  (2c  H — — ) e 0  (7.50) 

<  — 2clF  +  max  |(1  +  0-1).R2,  (1  +  0)(2c  +  “^)e0|  > 

from  which  it  follows  that 

TU>max  ii±^l(2c+^)e2  |  =t>  W  < -O.lcW .  (7.51) 

Then,  (7.21)  follows  because  \(f>\  <  \/2 W,  and  the  gains  are 

+  <7-52) 

Using  (7.34),  (7.41)  and  (7.52),  the  inner  and  outer  loop  small-gain  conditions,  (7.18) 
and  (7.19)  are,  respectively, 


2.1 


(1  +  g-1) 

0.95c 


2.1 


'(1  +  0) 


2  31 

+ 


0.95  3.8c 


<  1 


<  1 


(7.53) 

(7.54) 


Selecting  c  >  0  and  k  >  0  sufficiently  large  ensures  that  (7.53)  and  (7.54)  hold. 
Additional  freedom  for  the  selection  of  c  and  k  is  obtained  from  0  >  0,  which  allocates 
the  inner  and  outer  loop  gains. 


7.4  Summary 

We  have  discussed  how  the  new  observer  developed  in  preceding  chapters  can 
be  incorporated  in  output-feedback  control  design  with  control  laws  that  guarantee 
boundedness  in  the  presence  of  bounded  observer  errors.  The  combined  use  of  the  new 
observer  design  and  small-gain  control  design  tools  has  led  to  an  output-feedback  de¬ 
sign  procedure,  illustrated  on  the  jet  engine  compressor  example.  Such  combined  use 
of  observer  and  controller  design  tools  is  a  promising  research  direction  for  nonlinear 
output-feedback  control. 
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